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PREFACE 


The  High-Speed  Research  Program  and  NASA  Langley  Research  Center  sponsored  the  NASA 
High-Speed  Research  Program  Aerodynamic  Performance  Workshop  on  February  25-28, 
1997.  The  workshop  was  designed  to  bring  together  NASA  and  industry  High-Speed  Civil 
Transport  (HSCT)  Aerodynamic  Performance  technology  development  participants  in  areas  of 
Configuration  Aerodynamics  (transonic  and  supersonic  cruise  drag  prediction  and  minimiza- 
tion), High-Lift,  Flight  Controls,  Supersonic  Laminar  Flow  Control,  and  Sonic  Boom  Predic- 
tion. The  workshop  objectives  were  to  (1)  report  the  progress  and  status  of  HSCT  aerodynamic 
performance  technology  development;  (2)  disseminate  this  technology  within  the  appropriate 
technical  communities;  and  (3)  promote  synergy  among  the  scientist  and  engineers  working 
HSCT  aerodynamics.  In  particular,  single-  and  multi-point  optimized  HSCT  configurations  and 
HSCT  high-lift  system  performance  predictions  were  presented  along  with  executive  summa- 
rizes for  all  the  Aerodynamic  Performance  technology  areas. 

The  workshop  was  organized  in  three  sessions  as  follows: 

Session  I Plenary  Session 

Session  II  Independent  Session 

Session  III  Executive  Summaries 

The  proceedings  are  published  in  two  volumes: 

Volume  I,  Parts  1 and  2 Configuration  Aerodynamics 

Volume  II  High  Lift 

Conference  Chairmen:  Daniel  G.  Baize  and  Robert  L.  Calloway 

NASA  Langley  Research  Center 
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The  Configuration  Aerodynamics  technology  development  activities 
consist  of  four  primary  subtasks.  These  include: 

1 . Non-linear  rigid  and  aeroelastic  methods  adaption  and  validation. 

2.  Developing  and  validation  non-linear  aerodynamic  design 
optimization  capability. 

3.  Nacelle  design  and  airframe  integration  studies. 

4.  Assessments  of  the  baseline  TCA  configuration  and  determining 
the  benefits  of  the  Configuration  Aerodynamics  technology 
development  activities. 
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The  Non-Linear  Rigid  and  Aeroelastic  Methods  subtask  activities  consist  for 
four  primary  task  elements.  These  Include: 

1 . Prediction  of  complete  configuration  forces  and  moments. 

2.  Inviscid  aeroelastic  analysis  methods  development. 

3.  Viscous  aeroelastic  analysis  methods  development. 

4.  Prediction  of  power  induced  effects  on  the  complete  rigid  configuration. 


The  overall  objective  is  to  adapt,  apply  and  validate  non-linear  aerodynamic 
analysis  methods  for  prediction  of  the  forces,  moments  and  pressure 
distributions  on  complete  HSCT  type  configurations  including  elastic  and 
propulsion  system  power  induced  effects.  The  CFD  methods  are  used  as 
critical  elements  in  then  other  three  subtasks. 

Also  included  in  the  subtask  are  experimental  activities  conducted  to  validate 
the  analysis  tools. 


Boeing  did  not  have  any  activity  in  this  subtask  in  1996 
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The  Aerodynamic  Design  Optimization  Capability  subtask  consists  of  four 
subtask  elements.  These  are  the  four  essential  sequential  steps  necessary  to 
develop  the  ultimate  capability  of  non-linear  aerodynamic  design  optimization 
including  both  viscous  and  aeroelastic  effects. 


The  Boeing  Non-linear  cruise  point  activities  were  focused  on  including 
typical  realistic  design  constraints  in  the  TRAN  AIR  non-linear  optimization 
process  and  assessing  the  potential  cruise  performance  benefits  when  applied 
to  the  TCA  linear  design  baseline  configuration.  Robyn  Wittenburg  and  Rusty 
Conner  will  discuss  the  Boeing  point  design  optimization  work. 


The  Boeing  rigid  multi-point  design  optimization  activity  was  focused  on 
transonic  flap  optimization  of  the  TCA  configuration.  An  important  element  is 
to  understand  the  nature  of  the  flow  over  the  flaps  for  the  optimum  condition. 
This  is  a prelude  for  developing  multi -point  design  optimization  capability. 
Max  Kandula  will  summarize  the  work  performed  by  Dynacs  for  Boeing  in 
this  area. 
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4,  “Nacelle  and  Diverter  Integration  Studies  Mike  Malone  ( Northrup/G ruimnan ) 

5.  “Bifurcated  Inlet  Nacelle  Integration  and  Baseline  TCA  PAI  Assessment”  Steve  Chaney 


The  Nacelle  Design  and  Airframe  Integration  subtask  consists  of  four  three 
subtasks  elements: 

• Inlet  design  and  flow  assessment 

• Nozzle  design  and  boattail  drag  assessment 

• Nacelle  / Diverter  design  integration 


Mike  Malone  will  summarize  the  nacelle  and  diverter  integration  studies 
conducted  by  Northrup/Grumann  for  Boeing. 


Steve  Chaney  will  discuss  the  bifurcated  inlet  nacelle  integration  studies  and 
the  assessments  of  the  nacelle  / diverter  installation  on  the  baseline  TCA 
configuration. 
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The  Technology  Concept  Assessment  subtask  consists  of: 

• Aerodynamic  performance  assessments  of  the  baseline  TCA  configuration 

• Stability  and  control  assessments  of  the  TCA 

• Propulsion  induced  effects  test  programs  and  related  prediction  methods  validation 
activities. 

• Developing  methods  for  predicting  full  scale  aerodynamic  predictions. 

• Fundamental  aerodynamic  design  studies 

• Development  of  wind  tunnel  database  for  the  TCA 

• Aerodynamic  technology  metrics  developments  and  tracking  the  technology  process. 


Doug  Wilson  will  summarize  results  of  computing  aeroelastic  S&C  characteristics  using  the 
B oeing  AEOLAS  program. 


Bob  Kulfan  will  review  the  metrics  development  and  tracking  activities  during  the  plenary 
session  of  the  workshop. 
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The  Boeing  technology  development  areas  in  1 997  will  be  expanded  relative  to 
1996. 

Rigid  Configuration  Forces  and  Moments:  The  focus  will  be  on  prediction 
on  complete  configuration  forces  and  moments  including  aftbody  / empennage 
and  trim  drag.  This  will  support  the  planned  aft-body  closure  wind  tunnel  test 
program  planned  for  later  this  year. 

Rigid  Power  Effects:  The  objective  of  this  activity  is  to  develop  CFD 
modeling  and  analysis  capability  to  represent  Propulsion  Induced  Effects  (PIE) 
on  the  TCA  configuration.  Inlet  spillage,  bleed,  bypass  and  nozzle  flow  effects 
will  be  included.  The  initial  focus  will  be  on  identifying  the  magnitude  of 
these  effects  in  support  of  the  planned  Configuration  Aerodynamics  PIE  wind 
tunnel  test  program.  This  work  will  be  conducted  by  Northrup/  Grumann 

Non-Linear  Cruise  Point  Design:  .The  TRANAIR  design  optimization 
method  will  be  enhanced  to  include  optimization  with  viscous  effects. 

Complete  configuration  including  trim  drag  effects,  and  inlet  flow  quality 
constraints.  The  enhanced  methodology  will  be  applied  to  the  TCA 
configuration. 
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Rigid  Multi-Point  Design:  Techniques  will  be  explored  to  determine  if 
TRANAIR  can  be  used  to  conduct  meaningful  off-design  flap  optimization.  A 
strategy  for  conducting  multi-point  design  optimization  with  TRANAIR  will 
be  formulated. 

Inlet  Design  and  Flow  Assessment:  Studies  are  underway  to  develop 
techniques  to  model  inlet  spillage  effects  with  the  OVERFLOW  Navier  Stokes 
CFD  code.  The  objective  is  to  develop  and  validate  the  capability  to  model 
actual  flight  spillage  rates  so  that  the  axi-symmetric  and  bifurcated  nacelles  can 
be  accurately  compared  at  the  transonic  acceleration  flight  conditions.  This 
work  is  in  support  of  inlet  downselect  decisions. 

Nacelle/Diverter  Design  Integration:  This  will  be  an  extension  of  the  detailed 
nacelle  design  integration  studies  currently  underway.  The  planned  wind 
tunnel  test  programs  of  various  nacelle  integration  concepts  on  the  TC  A will 
be  supported. 

Aerodynamic  Performance  Assessment:  Wind  tunnel  performance  data 
obtained  with  the  TCA  configuration  will  be  evaluated  especially  at  the  off- 
design  conditions  in  which  the  outboard  flaps  are  deflected. 

Stability  and  Control  Assessments:  Stability  and  control  assessments  of  the 
TCA  configuration  will  be  continued. 

Propulsion  Induced  Effects:  Efforts  to  support  the  planned  PIE  wind  tunnel 
test  program  are  currently  under  way. 

Develop  Configuration  Aero  Wind  Tunnel  Database:  Boeing  will 
participate  in  supporting  the  development  of  the  wind  tunnel  database  for  the 
TCA  configuration. 

Aerodynamic  Technology  Metrics:  The  L/D  projection  processes  will  continue 
to  be  refined.  These  metrics  will  be  used  through  the  year  to  track  the 
technology  development  progress. 
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Boeing  Presentation  Schedule 

8:00  - 8: 10  * ‘Overview  of  CA  Activities  at  Boeing 

Bob  Kuifan 

8:10  - 8:50  “TCA  Cruise  Point  Design  Optimization'*  Robyn  Wittenberg 

8:50  - 9:30  “Observations  on  the  Process  and  Results  of  Optimization** 

Rusty  Conner 

9:30  - 10:00  “Optimum  Transonic  Flap  Studies** 

Max  Kandula 

10:00-10:15  Break 

10:15  - 10:45  “Nacelle  and  Diverter  Integration  Studies 

Mike  Malone 

10:45  - 11:45  “Bifurcated  Inlet  Nacdle  Integration  and  Baseline  TCA  PAI  Assessment** 

Steve  Chaney 

11:45-12:00  “Computation  of  Aeroelastic  SAC  Characteristics  Using  AEOLAS” 

Doug  Wilson 

This  is  the  schedule  for  the  Boeing  Presentations. 
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The  resulting  TRANAIR  non-linearly  optimized  configuration  is  taken,  as  the  final  step,  into  the  OVERFLOW 
code  for  Navier— Stokes  analysis  of  the  inviscid  design.  During  the  optimization,  checks  are  made  of  parameters 
which  we  have  found  can  indicate  potential  viscous  problems. 
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Years  of  experience,  on  multiple  configurations,  have  shown  that  except  for  viscous  effects,  TRANAIR’s  full 
potential  formulation  accurately  predicts  the  flow  about  and  forces  on  HSCT  geometries  at  supersonic  cruise. 
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process  was  based  on  subjective  opinion  as  to  what  looked  right.  The  wing  and  body  lofts  were  smoothed  by 
different  individuals  but  were  coordinated  for  consistency  of  approach  and  result. 
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Upwash  Plane  Layout : Closeup  of  Plane  Family 
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Surface  grids  for  Ref.  H & TCA 

» body,  inboard  wing,  collar,  box,  ellipsoid 
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Second  NASA/industry  High  Speed  Research  Configuration 

Aerodynamic  Workshop 


Bryan  W.  Westra,  Michael  B.  Malone,  and  Charles  C.  Peavey 
Northrop  Grumman  Military  Aircraft  Systems  Division 
Pico  Rivera,  California  90660 

The  1996  Nacelle/Diverter  Design  and  Airframe  Integration  Study  (WBS  4.3.13)  was 
initiated  to  develop  an  understanding  of  how  the  nacelle  and  diverter  integration  affect  the 
performance  of  the  High  Speed  Civil  Transport  Technology  Concept  Airplane  (TCA). 
Boeing  subcontracted  Northrop  Grumman  Corporation  (NGC)  to  assist  in  the  Assessment  of 
Nacelle  Integration  (WBS  4.3.13.1)  and  the  Nacelle/Diverter  Integration  (WBS  4.3.13.3). 
NGC’s  contribution  to  WBS  43.13.1  and  WBS  4.3.133  consisted  of  two  tasks:  a Nacelle 
Orientation  Study  and  a Diverter  Shape  Study.  The  computations  for  this  study  were 
performed  on  the  NAS  Cray  C-90  at  the  cost  of  approximately  300  CPU  hours. 

Both  tasks  involved  a parametric,  viscous  CFD  analysis  and  were  performed  on  the  full 
scale  TCA  at  Mach  2.4,  altitude  56,500  feet.  The  objective  of  the  Nacelle  Orientation  task 
was  to  provide  guidance  for  determining  the  minimum  aircraft  drag  and  minimum  inlet 
distortion  configurations.  For  this  task,  the  orientation  of  the  nacelle  was  changed  by  varying 
the  toe-in  angle,  pitch  angle,  and  diverter  height.  For  the  Diverter  Shape  Study,  the  objective 
was  to  provide  design  criteria  for  the  diverter  shape  and  wing  trailing  edge  to  nacelle  spacing 
to  minimize  the  installed,  2D  nozzle/nacelle  drag. 

The  initial  phase  of  the  Diverter  Shape  Study  was  to  reconcile  a viscous  drag  discrepancy 
observed  in  the  1995  Ref.  H Transonic  Nozzle  Boattail  Drag  Study.  It  was  suspected  that  the 
discrepancy  was  due  to  flow  solver  difference  between  OVERFLOW  and  GCNSfv.  It  was 
determined  that  the  solvers  agree  within  0.3%  (0.4  drag  counts)  and  that  the  discrepancy  was 
due  to  different  amounts  of  off-body  grid  stretching. 

For  the  Nacelle  Orientation  task,  eight  variants  were  developed  by  repositioning  the 
nacelle/diverters.  Two  variants  were  created  for  each  geometric  variable,  i.e.  nacelle  pitch, 
inboard  nacelle  toe  angle,  outboard  nacelle  toe  angle,  and  diverter  height.  Based  on  the 
variants,  it  was  predicted  that  a 0.68  count  reduction  could  be  achieved  by  pitching  both 
nacelles  down  and  toeing-out  the  inboard  nacelle.  This  optimum  configuration  was 
generated  and  resulted  in  a 0.66  count  reduction.  The  minimum  inlet  distortion  configuration 
was  achieved  by  primarily  two  adjustments:  an  inboard  nacelle  pitch  up  and  an  outboard 
nacelle  toe-in. 

Another  eight  variants  were  created  for  the  Diverter  Shape  task.  In  addition  to  the  baseline 
diverter  geometry,  a thin  diverter  (similar  to  the  Ref.  H diverter)  and  a NGC  developed  hybrid 
diverter  shape  were  examined  at  different  wing  trailing  edge  to  nacelle  spacings.  It  was 
concluded  the  best  diverter  shape  for  drag  reduction  was  the  hybrid  diverter.  The  best  trailing 
edge  to  nacelle  spacing  was  when  the  aft  end  of  the  nacelle  was  pitched  down  until  it  no 
longer  protruded  through  the  upper  surface  of  the  wing,  giving  a 0.65  drag  count  reduction. 

Finally,  the  results  of  both  studies  were  used  to  define  a configuration  that  reduced  drag  by 
a total  of  1.10  counts. 


1072 


Second  NASA/Industry  High  Speed  Research  Configuration 

Aerodynamic  Workshop 

Objectives 

Nacelle/Diverter  Design  and  Airframe  Integration  Task 
(WBS  4.3.1 .3)  

Assessment  of  Nacelle  Integration  (WBS  4.3.1 .3.1) 

Nacelle  Orientation 

• Perform  a Parametric,  Viscous  CFD  Analysis 

• Provide  Guidance  for  Determining  Minimum  Drag  and  Minimum  Inlet 

Distortion  Configuration  M , L1 

• Use  Nacelle  Toe-in,  Installation  Pitch  Angle  and  Diverter  Height  as 

Variable  Parameters 

Nacelle/Diverter  Integration  (WBS  4.3.1. 3.3) 

Diverter  Shape 

• Perform  a Parametric,  Viscous  CFD  Analysis 

• Provide  Guidance  for  Integration  of  Diverters,  Nacelles  and  Wing  for 
2D  Nozzle-Nacelles 

• Use  Diverter  Shape  and  Wing  T.E.  to  Nacelle  Spacing  (Pitch)  as 
Variables 

• Verify  OVERFLOW/GCNSfv  agreement  and  Determine  Boundary 
Layer  Grid  Stretching  Schemes  to  Obtain  the  Most 

Drag  Levels  — — 


In  1996  Boeing  subcontracted  Northrop  Grumman  Corporation  (NGC)  to  contribute 
to  the  Nacelle/Diverter  Design  and  Airframe  Integration  Task  (WBS  4.3. 1.3).  NGC 
assisted  in  the  Assessment  of  Nacelle  Integration  (WBS  4.3. 1.3.1)  and  Nacelle/Diverter 
Integration  (WBS  4.3.1.3.3)  tasks  by  performing  two  studies:  a Nacelle  Orientation  and 
a Diverter  Shape  Study.  Both  studies  were  performed  at  the  cruise  condition. 

The  objective  of  the  Nacelle  Orientation  study  was  to  perform  a parametric,  viscous 
CFD  analysis  to  provide  guidance  for  determining  the  minimum  aircraft  drag  and 
minimum  inlet  distortion  configurations.  For  this  task,  the  orientation  of  the  nacelle 
was  changed  by  varying  the  toe-in  and  pitch  angles.  Diverter  height  was  also  varied 
by  a vertical  displacement  of  the  nacelles. 

The  objective  of  the  Diverter  Shape  study  was  to  provide  design  criteria  for  the 
diverter  shape  and  wing  trailing  edge  to  nacelle  spacing  (varied  by  pitching  the  nacelle) 
to  minimize  the  installed,  2D  nozzle/nacelle  drag. 

The  initial  phase  of  the  Diverter  Shape  Study  was  to  reconcile  a viscous  drag 
discrepancy  observed  in  the  1995  Ref.  H Transonic  Nozzle  Boattail  Drag  Study.  This 
was  done  to  verify  that  the  NGC  flow  solver  (GCNSfv)  agrees  with  OVERFLOW  and 
to  determine  the  drag  effect  of  boundary  layer  grid  stretching  schemes. 
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Generalized  Compressible  Navier-Stokes  Code 

• NASA  Ames  ARC  Thin-Layer  Navier-Stokes  Algorithm 

• Implicit,  Node-Based  Finite-Volume  Scheme 

• Multi-Block  Structured  Grids  for  Complex  Geometries 

• Class  1 , 2 ,3,  & 4 Patched  Block  Interface  Mappings 

• Chimera  Overlapping  Grid  Block  Option 

• Grid  Sequencing  & Multigrid  Convergence  Acceleration 

• 2-Equation  Menter’s  SST,  Girimaji  ARSM, 

Spalart-Allmaras,  & Baldwin-Barth  Turbulence  Models 

• Extensive  Boundary  Condition  Menu 


MORTHROR  GRUMMAN 


The  CFD  code  used  was  GCNS,  developed  by  Northrop  Grumman.  It  is  based  on  the 
ARC3D  thin-layer  Navier-Stokes  algorithm  created  at  NASA  Ames.  The  convergence 
method  is  an  implicit,  node-based  finite-volume  scheme.  Complex  geometries  are 
analyzed  by  using  multi-block  structured  grids.  The  boundary  conditions  between 
blocks  can  be  specified  as  patched,  class  1 through  4,  where  class  1 is  point— to— point 
matching,  class  2 is  incremental  point-to-point  matching,  class  3 is  arbitrary  face 
matching,  and  class  4 is  arbitrary  sub-face  matching.  A Chimera  overlapping  grid  block 
option  is  also  available.  To  speed  steady— state  convergence,  grid  sequencing  and 
multigrid  convergence  schemes  can  be  used.  GCNS  provides  four  turbulence  models  to 
the  user:  the  Menter’s  k-to  SST  2-equation  model  (used  in  this  study),  the  Girimaji 
Algebraic  Reynolds  Stress  model,  the  Spalart-Allmaras  model,  and  the  Baldwin-Barth 
model.  GCNSfv  offers  many  boundary  conditions,  including  propulsion  specific 
conditions  such  as  characteristic  inflow  (mass  flow  ratio,  corrected  mass  flow,  and  inlet 
bleed)  and  outflow  (total  pressure  and  temperature)  conditions.  The  code  runs  at 
approximately  12  ps/iteration/grid-point  on  the  NAS  Cray  C-90,  and  parallelization 
allows  the  code  to  utilize  six  of  the  available  sixteen  processors,  allowing  effective  use 
of  the  multi-task  batch  queue. 


1074 


Second  NASA/Industry  High  Speed  Research  Configuration 

Aerodynamic  Workshop 


Solver  Comparison  & Viscous  Grid  Stretching 

Results  - Ref.  H 

• GCNS  & OVERFLOW  Agree  within  0.3%  / 0.4  cts.  CD 

• Grid  Stretching  (Using  GCNS)...  Notes: 


Notes: 

M = 2.4,  a = 4.4  deg. 

Results  do  not  include  forces 

on  body  aft  of  station  2,904' 

Turbulence  models:  Grid 

1-Menter  SSt,  Grids  2, 3, and 

4- Baldwin  Barth 

Grid  1 data  is  corrected  to 

Re=7,000,000 


Grid  1 


Right  Re 
Ay  = 0.0005* 
stretching  at 
surface 


Grid  2 


[ Boeing  Grid  j 
Ay  = 0.011' 

4 cells  const,  then 
minimal  stretching 


Grid  3 


Grid  4 


Ay  = 0.011' 

4 cells  const,  then 
significant  stretching 


0.011' 


stretching  at  surface 

NorrmKOf’  Grumman 


It  was  suspected  that  the  viscous  drag  discrepancy  observed  in  the  1995  Ref.  H 
Study  was  due  to  flow  solver  difference  between  OVERFLOW  and  GCNS.  After 
running  the  same  grid  in  both  solvers,  the  drag  results  were  within  0.3%  (0.4  drag 
counts).  The  boundary  layer  grid  stretching  schemes  were  then  investigated. 

Four  grid  schemes  were  used  to  determine  the  stretching  effects.  Grid  1 was 
taken  from  the  1995  Study  at  a flight  Reynolds  number  of  40  million.  The  second, 
third,  and  fourth  grids  shown  were  generated  for  a wind  tunnel  Reynolds  number 
of  7 million.  Grid  2 was  developed  by  Boeing. 

Grid  4 stretched  in  a similar  manner,  and  produced  results  similar  to  grid  1 after 
the  grid  1 results  were  corrected  for  Reynolds  number  differences.  The  small  drag 
discrepancy  is  likely  due  to  the  different  turbulence  model  used.  The  grids  with  a 
region  of  constant  off-body  spacing,  grids  2 and  3,  agree  well  with  the  test  results. 
This  determined  that  the  grid  stretching  near  the  wall  should  be  minimal  for  the 
first  several  off-body  cells.  It  is  important  to  note  that  stretching  after  the  region 
of  constant  cells  could  be  increased  without  a significant  change  in  drag  results. 
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Nacelle  Orientation  Variant  Matrix 


Variation 

1 

MBOARD  NACELLE 

OUTBOARD  NACELLE 

PITCH 

TOE-IN 

HEIGHT 

PITCH 

TOE-IN 

HEIGHT 

Baseline 

0 

0 

0 

0 

0 

0 

1 

tangent 

(+0.58*) 

0 

0 

tangent 

(+0.40°) 

0 

0 

2 

2xVar  1 
(+1.1 6} 

0 

0 

2xVar  1 
(+0.8Cf) 

0 

0 

3 

0 

0 

0 

0 

+1.0° 

0 

4 

0 

0 

0 

-1.0° 

0 

5 

0 

-1.0° 

0 

0 

Min  Cd 

0 

6 

0 

+1.0° 

0 

0 

Min  Cd 

0 

7 

0 

Min  Cd 

+0.25H 

0 

Min  Cd 

+0.25H 

8 

Klntae-  Inm 

0 

Min  Cq 

— 1 •!  _ 

+0.50H 

0 

Min  Cd 

+0.50H 

H is  the  inboard  diverter  height 


MO  ft  THRO/*  CftUMMAM 


Before  the  study  began,  a variant  matrix  was  developed  to  isolate  the  effects  of 
each  reposition  variable.  This  table  shows  the  nacelle  positions  for  the  eight 
variants  generated. 

The  first  two  variations  have  the  aft  end  of  the  nacelles  pitched  down.  Each 
nacelle  on  the  first  variation  was  pitched  until  the  aft  end  of  the  nacelle  no  longer 
penetrated  the  upper  surface  of  the  wing.  At  this  point  the  nacelle  upper  surface 
was  tangent  to  the  wing  trailing  edge.  The  second  variant  was  pitched  twice  the 
amount  used  for  variant  1.  The  outboard  toe-in  angle  was  varied  for 
configurations  3 and  4.  Variation  5 and  6 have  the  inboard  nacelle  toe-in  angle 
varied  and  the  outboard  nacelle  toe-in  fixed  at  the  minimum  drag  position 
determined  by  variant  3,  4,  or  baseline.  Finally,  the  diverter  height  was 
increased  by  dropping  the  nacelles  down  while  using  the  optimum  toe-in  angle 
for  each  nacelle. 
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An  axis  system  was  developed  for  each  nacelle.  The  origin  of  each  nacelle 
axes  system  was  located  at  the  leading  edge  of  the  diverter  where  it  intersects  the 
nacelle.  Two  rotational  axes,  pitch  and  toe  (yaw),  were  used  to  define  the 
nacelle  orientation. 

The  pitch  axis  direction  was  identical  to  the  isolated  nacelle  "y"  axis  (the 
nacelle  axis  before  it  is  repositioned  or  "rigged"  onto  the  wing).  A pitch  up 
rotation  was  defined  as  pitching  the  front  end  of  the  nacelle  up  and  aft  end  down. 
The  toe  axis  was  used  to  yaw  the  nacelle.  It  was  defined  as  a vector  normal  to 
the  plane  created  by  three  points:  the  coordinate  system  origin  and  the  inboard 
and  outboard  trailing  edge/diverter  intersection  points.  This  maintained  the  aft 
nacelle  elevation  at  the  trailing  edge  during  changes  in  toe  angle.  The  diverter 
was  rotated  with  the  nacelle  for  the  toe  cases  to  maintain  the  centerline 
relationship  between  the  components.  A positive  toe  rotation,  toe-in,  brings  the 
aft  end  of  the  nacelle  away  from  the  fuselage.  The  diverter  height  (at  diverter 
leading  edge)  was  changed  by  translating  the  nacelle  down  in  the  z direction. 
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The  nacelles  for  the  baseline  TCA  protrude  through  the  upper  wing  surface 
near  the  wing  trailing  edge.  When  the  aft  end  of  the  nacelle  was  pitched  down  to 
the  "tangent  to  wing  trailing  edge"  position,  the  nacelles  no  longer  intersected 
the  upper  wing  surface.  At  twice  this  pitch  amount,  a fairing  was  required  to 
close  the  gap  between  the  trailing  edge  and  the  nacelle.  A ramp  was  created 
which  extended  from  the  wing  trailing  edge  to  the  nacelle  upper  surface. 
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The  grid  consisted  of  5.3  million  grid  points  distributed  over  15  blocks.  Both 
patched  and  Chimera  (overset)  grid  methods  were  applied.  Flow  through 
nacelles  were  modeled  after  the  wind  tunnel  model. 

At  the  body  surface,  an  off-body  spacing  of  0.0005"  was  used.  After  post 
processing  the  baseline  TCA  results,  it  was  observed  that  the  values  of  y+  were 
less  than  3.0  everywhere  on  the  wing  surface.  This  should  be  adequate  for  the 
turbulence  model  (Menter’s  SST).  As  previously  mentioned,  the  off-body 
stretching  was  near  zero  for  the  first  4 cells  to  ensure  viscous  force  accuracy. 

Each  nacelle  was  defined  by  approximately  1.4  million  points  (2.8  million  for 
both  nacelles)  and  consisted  of  5 blocks.  The  nacelles  were  modeled  to  simulate 
the  wind  tunnel  model  by  incorporating  a flow  through  duct,  eliminating  the 
nozzle  ramp,  and  excluding  the  inlet  spike.  The  remaining  2.5  million  points 
were  used  to  model  the  aircraft  fuselage,  collar,  wing,  wing  tip  and  far  field. 
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A freestream  Mach  number  of  2.4,  angle  of  attack  of  3.0  degrees,  and  flight 
Reynolds  number  of  171.58  x 103  /in  (Re^c  = 195.50  x 106)  were  used  for  all 
runs.  This  Reynolds  number  was  obtained  from  the  full  scale  configuration  at  an 
altitude  of  approximately  56,500  ft  and  a free  stream  temperature  of  390.0°R. 
The  baseline  TCA  configuration  was  run  at  two  additional  angles  of  attack,  2.5 
and  3.5  degrees.  Forces  obtained  were  referenced  by  the  mean  aerodynamic 
chord  (MAC)  of  94.952  ft  and  wing  reference  area  of  8,500  ft2. 

A lift  coefficent  of  83.56  counts  and  a drag  coefficent  of  117.00  counts  was 
calculated  for  the  baseline  configuration. 
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Residual  History 


0.0 


-5.0 


ACDy 

acdp 


< 0.1  cnt.  / 50  iterations 

< 0.001  cnt.  / 50  iterations 


Coars 

e Gric 

^4 

-Fine  Grid 


200 


400 


600  800 
iterations 


1000  1200 


NAS  Cray  C-90 

17  CPU  hrs./Variant 
-300  CPU  hrs.  Total 


1400 

NORTHROP  GRUMMAN 


All  grids  for  this  study  were  run  in  GCNSfv  v3.6  for  an  identical  number  of 
iterations.  A sequenced  grid  method  was  used  to  reduce  the  computational  time. 
Eight  hundred  coarse  grid  iterations  were  initially  performed.  During  the  coarse 
run,  the  exhaust  blocks  failed  to  converge  due  to  the  lack  of  points  at  this 
sequence  level.  All  blocks  were  at  an  adequate  level  of  convergence  after 
running  600  fine  grid  iterations.  After  the  fine  grid  iterations,  the  viscous  drag 
changed  less  than  0.1  counts  per  50  iterations,  and  the  pressure  drag  changed 
less  than  0.001  counts  per  50  iterations.  All  variants  in  the  nacelle  orientation 
study  were  run  the  same  number  of  iterations.  The  total  CPU  time  required  for 
convergence  on  the  NAS  Cray  C-90  was  approximately  17  CPU  hours  per 
configuration.  Approximately  300  CPU  hours  were  used  for  both  studies. 
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Cross  Flow  vs.  Nacelle  Orientation 

-Obtained  from  Baseline  TCA  Flow  Field 


M=2.4,  a=3.0°,  ReMAC=1 95.5x1 06,  Tinf=390°R  A"9le  (d69) 


2.0  2.5 

MOftTHftOP  GRUMMAN  I 


The  predicted  optimum  orientation  to  obtain  minimum  inlet  crossflow  for  the 
inboard  nacelle  is  0.26  degrees  toe-in  and  1.76  degrees  pitch-up.  For  the 
outboard  nacelle,  the  predicted  optimum  orientation  to  minimize  crossflow  is 
1.24  degrees  toe-in  and  0.23  degrees  pitch-up. 

The  cross  flow  was  determined  by  moving  a disk  positioned  at  the  spike  tip 
through  the  baseline  TCA’s  flow  field.  It  was  not  necessary  to  reposition  the 
nacelle  grid  and  run  the  flow  solver  because  the  supersonic  flow  field  upstream 
of  the  nacelle  is  unaffected  by  the  nacelle  orientation.  The  cross  flow  (flow  in 
the  plane  of  the  disk)  for  the  inboard  and  outboard  nacelles  was  reduced  by  87% 
and  82%,  respectively.  Since  the  minimum  cross  flow  configuration  was  not 
modeled,  it  is  uncertain  what  amount  of  drag  penalty  must  be  paid  for  this 
minimum  cross  flow  configuration.  It  is  likely  that  the  drag  would  increase  by 
more  than  0.6  counts,  based  on  the  drag  increase  observed  for  the  variants  of  this 
study. 
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The  effect  of  each  parameter  on  drag  is  shown  in  the  above  cross  plots.  Pitch 
angles  76%  of  the  variant  1 amount  (i.e.,  0.44  deg.  inbd./0.31  deg.  otbd.)  are 
predicted  to  reduce  the  drag  by  0.16  counts  from  the  baseline  drag  polar.  The 
outboard  nacelle  should  be  toed-out  by  0.11  degrees  to  decrease  drag  an 
additional  0.04  count.  When  the  inboard  nacelle  is  toed-out  by  at  least  1.0 
degree,  a further  drag  reduction  of  0.54  counts  is  predicted.  The  diverter  height 
phase  of  the  study  showed  that  any  increase  in  diverter  height  increases  the  drag. 

If  we  assume  that  the  effects  of  pitch,  toe,  and  diverter  height  are  linear,  then 
the  principal  of  superposition  would  suggest  that  the  independent  effects  are 
additive.  Thus,  a configuration  that  uses  the  minimum  drag  value  of  each 
variable  would  be  the  minimum  drag  configuration,  and  reduce  drag  by  0.74 
counts. 

Note  that  the  minimum  drag  configuration  is  not  the  configuration  with 
minimum  cross  flow. 
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Reduced  Drag  Configuration 
Variant  9 Baseline 

Configuration 


Cp 


-0.2  -0.1  0.0  0.1  0.2  03  0.4  03 


Pitch  ->  ACD  = 0.14  coiin 

Inbd.  Nac.  Toe-out  ->  ACD  = 0.54  counts 
Superposition- > ACd  = 0.68  counts 

CFD  Model->  ACD  = 0.66  counts 

M=2.4,  a=3.0°,  ReMAC=1 95.5x1 06.  Tinf=390°R 


To  investigate  if  the  nacelle  variables  were  orthogonal  (no  cross  coupling)  with 
additive  effects,  a variant  9 was  modeled  where  both  nacelles  were  pitched  to  the 
"tangent  to  wing  trailing  edge"  (0.58  degree  inboard  / 0.40  degree  outboard) 
position  and  the  inboard  nacelle  was  toed-out  1 degree.  The  diverter  height  for 
both  nacelles  was  identical  to  baseline,  as  was  the  outboard  nacelle  toe  angle. 

From  earlier  configurations,  the  drag  decreased  by  0.14  counts  due  to  pitch  and 
0.54  counts  due  to  inboard  nacelle  toe-out.  Based  on  superposition,  a drag 
decrease  of  0.68  counts  was  predicted  for  the  new  geometry.  CFD  results 
revealed  that  the  variant  9 configuration  decreased  drag  0.66  counts.  Thus,  it 
can  be  concluded  that  the  superposition/additive  method  is  valid  for  the  range  of 
angle  variations  examined. 
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Diverter  Shape  Variant  Matrix 


Variation 

Diverter  Shape 

Nacelle  Pitch 

Baseline 

Wide  (TCA  geometry) 

0 

1 

Wide  (TCA  geometry) 

tangent  at  T.E. 

2 

Wide  (TCA  geometry) 

2.0  x variant  1 

10 

Thin  (Ref.  H.  geometry) 

0 

11 

Thin  (Ref.  H.  geometry) 

tangent  at  T.E. 

12 

Thin  (Ref.  H.  geometry) 

2.0  x variant  1 

13 

Hybrid 

0 

14 

Hybrid 

tangent  at  T.E. 

15 

Hybrid 

2.0  x variant  1 

Notes:  Increments  are  from  the  baseline  geometry. 


NOR  THROR  GRUMMAN 


For  the  Diverter  Shape  Study,  8 variants  were  generated.  Note  that  the 
numbering  sequence  continues  from  the  variants  developed  in  the  Nacelle 
Orientation  section  of  this  report.  Three  diverter  shapes  were  used  in  this  study. 

The  wide  diverter  is  70.6  inches  wide  and  is  the  baseline  TCA  diverter.  It  has  a 
leading  edge  angle  of  22  degrees. 

The  thin  diverter  was  based  on  the  Ref.  H diverter.  At  32  inches  wide  with  a 
17  degree  leading  edge,  it  represents  the  minimum  volume  configuration, 
bounded  by  structural  requirements. 

The  hybrid  diverter  begins  with  a 17  degree  leading  edge,  but  expands  to  the 
wide  geometry  at  the  trailing  edge  by  following  a parabolic  curve  . 

The  wing  trailing  edge  to  nacelle  spacing  was  changed  by  pitching  the  nacelle 
as  in  the  Nacelle  Orientation  study.  As  previously  discussed.  The  "tangent  at 
wing  trailing  edge"  condition  occurs  when  the  nacelle  is  pitched  until  the  aft  end 
no  longer  penetrates  the  upper  surface  of  the  wing  . 
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Geometry  Development 


Thin  Diverter  Wing  Lower 
Boundary  Surface 


Nacelle  Upper 
Surface 
View  A-A 
Aft,  looking  Fwd. 
(rotated  90  deg) 


Wide- 
Thin  - 

Hvbrid . 


22° 


Fairing 

Wide  Diverter 
Boundary 

Hybrid  Diverter 
Boundary 


Wing  Lower  Surface  / 
Nacelle  Intersection  at 
"tangent"  position 


NORTMROR  GRUMMAN 


A plan  view  shows  the  diverter  geometry  for  the  condition  where  the  nacelle  is 
pitched  tangent  to  the  wing  trailing  edge.  For  the  thin  diverter  geometry,  the 
nacelles  intersected  the  wing  lower  surface  near  the  trailing  edge  outside  the 
diverter  bounds.  In  many  cases,  the  resulting  geometry  had  a thin  "V''-channel 
created  where  the  nacelle  upper  surface  and  wing  lower  surface  met  at  sharp 
angles.  This  channel  was  removed  by  creating  a small  fairing,  less  than  1 inch 
tall,  which  spanned  from  the  wing  lower  surface  to  the  nacelle. 

The  hybrid  diverter  concept  arose  from  observing  the  flow  near  the  thin 
diverter  fairing.  Specifically,  there  were  regions  of  high  viscous  stress  where 
the  fairing  met  the  thin  diverter.  These  regions  were  eliminated  by  constructing 
the  hybrid  diverter. 
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Diverter  Shape  Results 


Relative  Pitch  Amount 

MOftTMftOf*  GRUMMAN 

M=2.4,  a=3.0°,  ReMAC=1 95.5x1 06,  Tinf=390oR 

The  thin  diverter  did  not  improve  from  its  wide  diverter  equivalent  until  the 
nacelle  was  pitched  to  the  "2x  tangent”  position.  At  this  pitch  angle,  the  thin 
diverter  was  more  beneficial  than  the  wide  diverter  because  the  inboard  diverter 
becomes  a true  pylon,  supporting  the  nacelle  completely  off  of  the  wing,  and  the 
flow  does  not  "pinch-off'  between  the  wing  and  nacelle.  For  this  configuration, 
the  drag  decreased  by  0.11  counts  from  the  baseline  wide  diverter  at  the  "Ox 
tangent"  position.  Note  that  the  drag  is  0.36  counts  lower  than  the  "2x  tangent" 
position  wide  diverter. 

The  hybrid  diverter  performed  best  at  reducing  drag.  It  decreased  the  drag  by 
0.40  counts  from  the  baseline  TCA  configuration.  At  the  "lx  tangent"  pitch,  the 
drag  was  reduced  by  0.65  counts.  Recall  from  the  nacelle  orientation  study  that 
the  wide  diverter  only  decreased  drag  by  0.14  counts  for  the  "lx  tangent"  pitch. 


1087 


Second  NASA/Industry  High  Speed  Research  Configuration 

Aerodynamic  Workshop 


Combined  Results 

Nacelle  Orientation  / Diverter  Shape 


It  was  decided  to  use  the  best  diverter  shape  found  in  the  diverter  shape  study 
and  combine  it  with  the  best  nacelle  orientation  observed  in  the  nacelle 
orientation  study.  It  was  predicted  that  the  new  geometry  would  combine  the 
drag  benefit  of  0.66  counts  from  the  nacelle  orientation  with  the  0.40  counts 
from  the  hybrid  diverter.  This  would  provide  a 1.06  count  reduction  in  drag. 
When  the  combined  geometry  (variation  16)  was  run,  the  result  was  a drag 
reduction  from  the  baseline  drag  polar  of  1.10  counts. 
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Nacelle  Orientation:  Conclusions 

• Minimum  Cross  Flow  Configuration  - Determined 

Inboard  Nacelle:  Primarily  Pitch  Up 
Outboard  Nacelle:  Primarily  Toe-in 

• Minimum  Drag  Configuration  - Determined 

Pitch  Change  - decreased  drag  0.14  counts 
Outboard  Nacelle  Toe  Change  - no  significant  change 
Inboard  Nacelle  Toe  Change  - dec.  drag  0.54  counts 
Diverter  Height  Increase  - increased  drag 

• Method  of  Isolated  Variable  Effects  - Validated 

ACD(pitch)  + ACD(toe)  = ACD(pitch  & toe) 


NOftTHROI*  GRUMMAN 


The  nacelle  shape  study  determined  the  minimum  drag  and  minimum  inlet 
distortion  configurations  for  the  TCA  geometry.  The  minimum  cross  flow 
configuration  was  obtained  by  primarily  two  adjustments:  an  inboard  nacelle 
pitch-up  and  an  outboard  nacelle  toe-in. 

The  minimum  drag  configuration  can  be  achieved  by  pitching  both  nacelles  up 
and  a toe  out  adjustment  of  the  inboard  nacelle.  This  will  reduce  the  drag  by 
approximately  0.7  counts. 

Separate  pitch  and  inboard  toe  variations  both  reduced  drag.  When  a 
configuration  was  modeled  that  contained  both  the  pitch  and  inboard  toe  change, 
the  resulting  drag  decrease  was  nearly  identical  to  the  sum  of  the  separate  drag 
reductions.  Thus,  the  isolated  effects  of  each  variable  can  be  considered  additive 
for  small  variable  increments.  Future  studies  could  use  this  approach  provided 
the  magnitude  of  the  changes  were  similar  to  those  used  in  this  study. 
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Diverter  Shape:  Conclusions 

• Effect  of  Diverter  Shape  on  Drag 

-Thin  With  "Pinched"  Flow  (Ox  and  lx  tangent): 
increased  from  baseline  equivalents 

-Thin  Without  "Pinched"  Flow  at  Inbd.  Nacelle  (2x  tangent): 

0.36  cnt.  reduction  from  wide  equivalent 
-Hybrid: 

0.40  to  0.56  cnt.  reduction  from  wide  diverter  equivalents 

• Effect  of  Nacelle  to  Wing  Trailing  Edge  Spacing  on  Drag 

-Best  Position  for  Wide  and  Hybrid  Diverters:  Tangent 
0.14  cnt.  reduction  for  wide  diverter  from  baseline  polar 
0.65  cnt.  reduction  for  hybrid  diverter  from  baseline  polar 
-Best  Position  for  Thin  Diverter:  2x  Tangent 
0.1 1 cnt.  reduction  from  baseline  polar 

Combined  Studies:  Conclusion 

Decreased  Drag  By  1.10  counts 


This  study  showed  that  the  thin  diverter  configuration  had  higher  drag  than  the 
wide  diverter  configuration  when  the  flow  in  the  diverter  channel  was  pinched 
off.  When  the  aft  end  of  the  nacelle  was  pitched  enough  to  eliminate  the 
pinch-off,  the  thin  diverter  reduced  the  drag  by  0.36  counts  from  a wide  diverter 
at  the  same  pitch  angle.  The  hybrid  diverter  proved  to  be  the  diverter  that 
resulted  in  the  least  drag,  providing  a reduction  of  0.40  to  0.56  counts  from  the 
wide  diverter  cases. 

The  best  pitch  angle  observed  for  the  wide  and  hybrid  diverters  was  the  "lx 
tangent"  value.  At  this  position,  the  wide  diverter  reduced  the  drag  by  0.14 
counts  and  the  hybrid  by  0.65  counts.  The  best  pitch  angle  observed  for  the  thin 
diverter  was  at  twice  this  pitch  angle.  The  resulting  drag  reduction  was  0.11 
counts  from  the  baseline  drag  polar. 

When  the  hybrid  diverter  was  modeled  with  the  best  observed  nacelle  pitch  and 
toe  angles  from  the  Nacelle  Orientation  and  Diverter  Shape  studies,  the  new 
configuration  reduced  the  drag  by  1.10  counts  from  the  baseline  TCA.. 
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Recommendations  for  Future  Work 


Investigate  the  Drag  Effects  of  Inlet  Spike  with  Optimum  Nacelle 
Orientation  and  Diverter  Shape 

Evaluate  the  Drag  Effects  of  Inlet  Spillage  and  Inlet  Bleed  Flows. 


NORTHROR  GRUMMAi*  \ 


This  study  made  use  of  a flow  through  nacelle  without  an  inlet  spike.  The 
effect  of  the  spike  on  drag  after  the  nacelle  orientation  and  diverter  shape  have 
been  optimized  should  be  investigated. 

Likewise,  the  drag  effects  of  inlet  spillage  and  bleed  extraction/exhaust  should 
be  evaluated  for  design  and  off— design  conditions. 
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Outline 
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ii 


Each  configuration  and  condition  assessed  using  surface  pressures, 
integrated  forces,  inlet  flow  field  chararacteristics. 
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No-slip  walls  are  all  turbulent. 
Central  difference. 
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bomnc  hsr.  1997  Aerodynamic  Performance  Workshop 

HSCT  High  Speed  Aerodynamics 

TCA  Wing/Body  OVERFLOW  Grid 


Grid  points  removed  for  clarity. 
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HSR  1997  Aerodynamic  Performance  Workshop 


Flight:  Mach  2.4,  Mid-cruise,  altitude^^SOCHt 

ReMAC  = 195  million, 

NOTES:  Three  equidistant  points  from  surface. 

Stretching  <1.2 

Wing  surface  has  106  points  from  LE  to  TE: 

Asle  = .01%  to  .02%  chord,  AsTE  ■ .4%  to  .5%  chord. 
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Flight:  Mach  2.4,  Mid-cruise,  altitude  = 56.500  ft 

^®mac  = 195  million, 

NOTES:  Three  equidistant  points  from  surface. 

Stretching  <1.2 

Wing  surface  has  106  points  from  LE  to  TE: 

Asle  = .01%  to  .02%  chord,  AsTe  = .4%  to  .5%  chord. 
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amics 

Wing  / Nozzle  Transition  Geometry 


Inboard  Nacelle  Installation 
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Inboard  Nacelle  Installation 
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Mach  0.90,  ReMAC  = 10.8  million 
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'namlcs 

TCA  Baseline  PAI  Assessment  Summary 
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macn  nignt  Keynoids  number  OVERFLOW  analysis  of  W/B/N/D. 

Continue  development  of  spillage  modelling  on  isolated  nacelle. 

Transonic  wind  tunnel  and  flight  analysis  of  W/B/N/D  with  and  without  spillage. 
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amics 

Inlet  Mach  Number  for  TCA  and  Optimized  TCA 
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1191 


Mike  Elzey 
Brian  Nishida 
Christine  Titzer 

Ross  Sheckler  (Dynacs  Engineering) 


HSR  Aerodynamic  Performance  Workshop 


1 


i 

s 

0) 

1 

C 

5 

£ 

1 

(0 

K* 

O 

S: 


LU 

CO 

CO 

LU 

CO 

CO 

< 

mJ 

o 

GC 


O 

O 

O 

z 

< 

>- 

t 

ffi 

2 

CO 

< 

o 

H 

CM 

■ 

■ 

■ ■ ■ 
00  w 

§ 

« -5 

m « 

5 5 


£■0 

g>  c 

_0  (0 

o 

^ .2 
o w 
a>  >> 
© 

CD  g 
r-  CO 


O CO 

0)  g 
.2-2 

a>  o. 

o E 
co  o 
o 


CD 

£ L. 

+*  +*  <D 

w © ‘jz 

a>  g <d 

(o  c a 

CO  O x 

< O © 


2c 
c w 
0 o 


T3 

C 

CO 


o 

CO 

c 

CO 


CO 


1! 

8- 

Eg 

0 o . 

>*:©  D) 

!q  © © 
© > © 

co  — 

■O  C g 

® O)  g 

® © 3 
CL  <i)  *- 

• -S  -O 


5 « 

£ CO 

<0  c 
o 

h-  3* 
O CO 
Z 

© 2 
« > 
© 

n c 

O) 

© © 
+2  © 


__  ■© 
© c 
c n 
c 

3 Q) 
*-  C 

■o  © 

£ m 

T © 

£ n 

TJ  © 

c £ 
© ** 


Q 

LL 

O 


ra  © © 


s 4 

M—  V 
© 

© ® 

£ © 


® © 

© 3 

© © 

c n 
© 52 

0 o 


o 

© 

© 

Q. 

3 

© 


D) 

Q. 

© 


■t-  CO 

© o *- 

■c  *3  3 


© 

© 


© 

O 

■MM 

4-* 

(0 


a> 

4-« 


a> 

E 

a 

o 

© 

> 

© 

•o 

c 

o 

© 


o 

© 

© 


c 

o 

o 

■o 

c 

© 


O) 


o 
o 

>» 

O) 

o 

0 

c 

€ 

£ 

© _ 

1 § 


X3 

© 

© 

§ 

o 

L. 

a 

E 


© 


© 

z 

© 

> 


£ 5> 
© 
© 
73 

Q. 

O 

© 
> 
© 
Q 


h- 

o> 

CD 


© 

© 

O 

© 

■ MM 

© 

« 

n 

■■M 

V. 

© 

© 

© 

© 

O 

© 

O 

© 

© 

© 

ly 

4-4 

£ 

£ 

c 

© 

O 

© 

© 

O 

E 

■^M 

■ ■ 

□ 

O 

00 

0) 

© 

Q. 

w 

© 

LL 

O 

O 

■MM 

X 

10 

O 

« 

co 

0 

00 

00 

■D 

(/) 

0) 

0 

0) 

> 

< 

3L. 

<D 

< 

O 

< 

O 

N 

O) 

H 

© 

1- 

o> 

CL 

O 

3 

Q. 

Cl 

O 

1 

© 

> 

E 

© 

> 

CO 

© 

0 

© 

o> 

O 

O 

Q 

o> 

T" 

• 

• 

• 

1192 


o 

<0 

0 

n 

0 

HM 

0 

Q 

0 

c 

0) 

E 

■■■■ 

5 

Q. 

X 

LU 

O 

66 

</) 

< 

o 

H 

Q. 

O 

<D 

> 

0) 

Q 


0 

~o 

■0 

c 

0 

E 

0 

0 

P 

0 

5 

a. 

4-* 

c 

0 

D 

0 

O 

< 

■ MM 

c 

0 

O 

0 

V. 

0 

0 

H— 

a 

c 

3 

0 

0 

'£5 

k. 

0 

B^M 

c 

■ MM 

H— 

*fr- 

o> 

c 

0 

73 

‘E 

0 

c 

■O 

0 

O 

a 

E 

■a 

■0 

c 

c 

0 

0 

0 

0 

c 

c 

0 

0 

E 

E 

£ 

0 

1- 

3 0 

3 

o’ 0 

cr 

0 Si 

0 

c 

£ 

■-  CD 

■ MM 

■0 

■O  . 

0 w 

0 ^ 

0 .2 

m 0 
0 (v| 

Q.  C 

Q-1- 

O 8 

0 0 

« 0 

•—  "O 

r c 

t 0 

0 0 
a.  J= 

0 ^ 
CL  ^ 

• 

• 

■D 

C 

0 

t 

o 

a 

a 

3 

0 

o> 

c 

B^H 

4-J 

0 

a> 

.2  (0 
C (0 

o 0 

0 ^ 
c g> 

(o  2 

* a. 

o.E 

= c 

3 £ 

g 3 

O « 

0 4S 

■ m\ 

&1  O 

4-»  .. 

0 ■{- 

3 CO 

^ C 
c 0 

p 

1 8 

a.  T3 


CO 

0 

£ « 
3 t- 
CO  m 

O ^ 

><£ 

■n  CD 

0 g 

1 j. 

*•  0 

CO  0) 

u 4-* 

0 co 

SI 

.£  3= 
c 0 

§■§ 

1 = 

CO  -T 

c 0 c 

c .2 
0 0 

E 3 
0 a. 
.i:  o 

0 *D 
*“  C 
C 0 

■ mm 

o>P 

C ■ 

■ MM  ■ 

CO  CO 
.2-^ 
.2  co 
tr  co 

co  > 
a.  c 


1193 


2/26/97 


0)  - 
CD  0 

"O  0 

»§■ 

*-  0 

I « 

1 .£ 

- — (!) 

8 S 

xjS 

S 3 

CO  0) 
0 0 

cS  — 
■o  <c 

"5  £ 

•H  V 

S| 

£ « 
■z  Q. 
0 X 
Q.  0 

g £ 


Q 0 
LL  Z 

O CO 

0 CO 
0 Jl 
3 O 


w 0 

T3  O 
O '4= 

£ ■- 

2 S 
r o 

co  2 
0 0 
C £ 

= o 
o>  o 

C +3 

g 8 

i : 

c ts 

fl)  0 

ii 

o 8 

£ o 

!>« 

O)  ■*- 

■—  0 
*r  ■o 

o O) 

•H  V- 
0 r- 
0 ^ 

® i 

j>  0) 
(0  c 

3 2 

£ E 

O 8 


1194 


2/26/97 


& 

■c 

I 

i 

Q) 

O 

s 

E 

€ 

£ 


6 

<8 

I 

i 

cc 

S: 


I 


I 

i 

I 

i 

1 

c 

o 

o 

<* 

t 

s 

a 

to 

I— 

o 

a: 


O 

© 

O 

Q. 

Q. 

< 


■a 

o 

JO 

© 

£ 

© 

c 

© 

a. 

CM 

0 

in 

< 

U)  . 

.£  aT 
a>  o 

3 C 

(o  £ 
© ® 

zi 

si 

1 ^ 

B o 


■g  © 

A ■* 

o > 

CO  "m 

c c 

.2  | 

o 8 

5 o) 

t >* 

O -U 

o o 

0) 

- CB 

*N 

a ■ 

— O) 
CB  „r 

O >£- 


> 

73 

£ 

© 

o aj 

o -S 

I« 

b>  © 

£■  3 

.2  LU 

if 

o co 
O 3 


co  II 

.2  2 

<0  J= 
T CB 

© g 

o ^ 

S D> 

© fl> 
° ^ 

* 1 

2o  8 

TJ  ^ 
'5)  CB 


0 <D 

ff 

£ Q. 
O X 
O © 


3 > 

_r  jq 


(0 

< 

-I 

o 
111 
< 

o> 
c 

'co 

3 

■e  ^ ^ co 

#n  ■ ■ 


0)  Ml 

E 2 

CM  LL 


O 

B^B 

4-« 

o 


<UJ1 


|5 

TJ  0) 

£ » 
co  CB  - 

5 f o> 


CB 


(D 


a -g 


CO 

3 


c 

g CO 

£ c 
a>  o 

o E 

£ 73 

— c 
■O  o 

’5)  ° 


.2  E co 
o „ © 


w w 

SZ  CO  *n 
3=  CO  - 
0)  CB 
O P 
o c 

d>  73 
O 
C 


. S>  co 

m C0  — C0 


O 
c 

a> 
o 

k.  J5 
CB  % 


O © _ 

T ? S 

« § £ 
So! 

JO  JQ  g 
© © _ 

© £ fe 

+■*  «£  TO 

© » a> 

— © c 

3 <|)  .= 
O _l 

75  = 

O . . 


© 


c c 


© 

o 


3=  ■o 
■“  © 


o 

3 


o _ 

E - - 

2 © "O 
>1  3 l 
■g  O >- 

P 3 © 

© ^ *2 
< (/>  o 


© 

N 

■ w^m 

k. 

© 

© 

c 

BBB 

T 

£ 

o 

c 


» 

CO 

§ 


© 


■o 

IB 


o 

© .£ 
■F  +■* 

© © 

■o  £ 

V.  ___ 

m CO 

g £ 

.£  o 

c °> 
o o) 
£ £ 
© 

■r  73 

O Jj 

Z §■ 

s ® 

o »- 

^ © 
O k. 

© © 
u £ 

o £ 

° w 
o c 
33  O 

2 += 
© o 

© £ 

CB  O 
© O 

.£  o 

mm  im 

4-* 

■U  w 
73  J5 
< © 


c 

o 

B^B 

© 

3 

£ 

© 

c 

■ tmm 

■D 

© 

© 

3 

© 

£ 

o 

</> 

X 

a 

LL 

To 

c 

o 

© 

c 

© 

£ 

■ Bi 

■o 

■ 

£ . 
.£  o 

S o 

© 

■°  £ 
© " 

m 10 
© Q) 

3 

■°  £ 
© 5 

h-  ^ 


1195 


HSR  Aerodynamic  Performance  Workshop 


1196 


Variables  in  database  (not  all  variables  for  each  term): 

Mach,  alpha,  beta,  control  surface  deflections  (~11),  dynamic  pressure, 
mass  case 
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Symbolic  Development  of  the  Adjoint  Method 

The  next  two  slides  show  the  basic  analysis  behind  the  adjoint  approach  as  outlined  under  the  Rationale  section  above.  The  following 
references  provide  a more  complete  treatment  and  cover  the  actual  application  to  the  governing  Euler  equations: 

1)  J.  Reuther  and  A.  Jameson.  Aerodynamic  shape  optimization  of  wing  and  wing-body  configurations  using  control  theory.  AIAA 
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Symbolic  Development  of  the  Adjoint  Method 
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Single-Block  Design  Approach  ( SYN87-SB ) , continued 
Numerical  Optimization  Packages 

Unconstrained  Ouasi-Newton  Algorithm  (ONMD1F7 ) 
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Single-Block  Design  Approach  (SYN87-SB) , continued 
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Single-Block  Design  Approach  (SYN87-SB) , continued 
Single-Block  Wing/Body  C-H  Mesh  Generation,  continued 

Grid  Perturbation,  continued 
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Single-Block  Design  Approach  (SYN87-SB) , continued 
Pseudo  Nacelle  Effects  (Supersonic  Applications  Only) 
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Single-Block  Design  Approach  (Continued) 
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• Linear  Constraints 

- Wing  Section  T/C  (or  Y)  at  Specified  X/C 

” Cabin  Section  Radius  at  Specified  Polar  Angle  (Body  Camber  Fixed) 


Pseudo  Nacelle  Effects,  ccontinued 
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points  on  the  upper  aft  nozzles  project  into  the  upper  half  of  the  C-H  mesh,  whereas  the  scheme  must  pick  the  nearest  grid  cell  trom 
the  lower  half  only.  Any  such  search  must  check  every  cell  of  the  lower  half  grid  before  concluding  that  the  target  point  is  not 
contained  in  any  of  them,  and  this  is  expensive.  Updating  the  buoyancy  corrections  for  each  line-search  solution  takes  about  1.8  CPU 
seconds,  but  basically  after  setting  it  up  the  pseudo  nacelle  simulation  has  negligible  further  cost. 
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Single-Block  Design  Approach  (SYN87-SB) , continued 


2 

e 


CQ 

U 

-w 

CO 

G 

O 

u 

u 

CQ 

cy 

fl 


sc 


<u 

G 

O 

T3 


si 


CO 

" G *G 


D T)  T3 
> G 

'3  « 

G on 
1)  co 

a,  ja 

x S 

(D  ^ 
«.2 
E ■= 


« js* 

-*—>  CO 
cO  on 

o U 
$ H 
w a 

<z>  .S 


. § $> 
>:§  § 
dfc3*5 
.2  3 8 g « 

u M u £ 

- .2  •§  IT  £* 

•o-g  73 

g *5.9*3 

CO  g X)  3 03 

.2  -a  3 E .22 

.tS  es  ,e  « 
c.S-  •■- 


CO 

3 

cr 


<U 

«*o 

<D 


CO  •* 
a>  t3 
c 0) 

.G  N 


.2  o 

H c 
o ^ 


co  1 

G <D  £ 

. ^ cl 

o 2® 

«5  <0.2.9 

sS-S  &•* 

5 S o *-*  <l> 
*C  ^ H fi 

CO  iu  > 

> .fa  l co  2 

W)3  0-2- 

C C X >> 

■3  2 o .2  P 

g o'"0  5 
fa  C >>  > 

s| 

35 


a) 

g 

o 

~ <u 

W)  two 


CO 

<U 

•a 

*> 

o 


<u 

a 

c 

o 

co 

cO 


CO 

g 

o J& 

G ** 

CO  y .S 

G *G 
O ~ 2 

• — * p-G 

co  a ~ 


_ G O 

0.4=  a.  “ g 
m 3 ° x?  § 

Vi  ^-2  b « 

*>  « a -g 


r- 

oo 

Z 


- w c 
.EEs 
?a  o 


f->  o __  r* 

>*  to  « U £ 

&0  «2  O c ^ 


J 

O 


C 

o . 

CJ  CO 

w ih« 

00  G X) 
Dm  g CO 

Zs''" 

- § 

G 


cO 


ft, 

CO 


co 

C o 
O O 

- °.E 


G 

C 


cO 


CO 

C 

O 


2 O 


tH  CO  <4-h 


^ c 

G CO 
cO  O __ 

9 “ u . 

- . . .2  G 

r o ~o  v 
d g „rS 

O U CO 
G G S"1  H 


D,  <D 

U-i 

O X 

.g  0/2 

Il'i'SJS 

5 E 

c O -r 
.9  X)  & 


«.2x 

o " 


2 £5 


cO 


G D cO 

S 3 £ 

hco 


G3  co 

- .2  ^ 
.S3  c-1 
««2 


1223 


Single-Block  Design  Approach  (SYN87-SB) , continued 
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- Fuselage  Area  Ruling  (Fixed  Section  Shape) 

- Fuselage  Section  Shape 

- Cabin  Floor  Kink  Point  Locations 


Single-Block  Design  Approach  (SYN87-SB) , continued 
SYN87  Design  Variables 
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The  body  section  area  variables  preserve  section  shape,  using  an  iterative  procedure  to  achieve  the  specified  streamwise  area  ruling 
The  fuselage  volume  may  be  conserved  implicitly  by  a choice  of  2 streamwise  perturbations  which  add  zero  net  cross-sectional  area. 
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Improvements  to  SYN87-SB  since  7/96 
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Improvements  to  SYN87-SB  Since  7/96 
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mprovements  to  SYN87-SB  Since  7/96  (Continued) 


on 

<D  w 

£ w) 

1 1 
V-i 

fc;  S 

43  H 


o 

P<10 

j+h 

T3  PJ 

<D 

*5 

• t-h 

T3  > 
O O 

00  S u 

| | s 

^ ti 

^ S3  <D 

S o > 

£ g a 

o o w 

J *s  ^ 

43  ^ 43 

o o 73 

a W .5 

S 9 ^ 

'T  S « 
- 5 T3 
ss  5 u 

e S z 

<u  — 

> a >» 

•r-i  r*  »-( 

£ -a  | 

0 g ! 

CO  W ^ 

mi  E D 

T3  O Vh 

W 43  B 

h m 

8 a.  > 
£*  rj  *Sh 
o<  ^ n 

^ u fed 
t-  c3  « 
§ So 

^ § & 
< £ z 

1 I 


1231 


Improvements  to  SYN87-SB  since  7/96,  continued 
Fortified  Body  Design  Capability 
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Improvements  to  SYN87-SB  Since  7/96  (Continued) 
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• Generalized  Lofting  of  Wing  Section  Perturbations 

- Hicks-Henne  Shape  Function  Lofting  + Original  Polynomial  Lofting 


Improvements  to  SYN87-SB  Since  7/96  (Continued) 
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- Option  to  Read  Grid  + Flow  Field  Speeds  Drag  Polar  Calculations 

- Nacelle/Diverter  Lift/Drag  Breakdown  + Linear  Constraint  Details 

- PostScript  Plots  May  Include  Spanwise  and  Chordwise  Surface  Cuts  + 
Overlays  + Cabin  Constraint  Polygons 


Improvements  to  SYN87-SB  since  7/96,  continued 
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• New  routine  ELLIP3D  smooths  the  sub-block  interiors  with  full  spacing  and  orthogonality  control  as  for  ELLIP2D. 

For  subsequent  perturbed  grids,  WARPQ3D  and  WARP3D  are  applied  to  the  various  sub-block  faces  and  interiors  (very  efficient). 


Improvements  to  SYN87-SB  Since  7/96  (Continued) 
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* Cabin  Floor/Wing  Spar  Clearance 
Generalized  Lofting  of  Wing  Section  Perturbations 
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Improvements  to  SYN87-SB  since  7/96,  continued 
Improved  Input/Output 
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Improvements  to  SYN87-SB  since  7/96,  continued 
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The  perturbation  scheme  required  a refinement  for  the  same  reason  as  indicated  above  for  the  single-block  wing  surface  grid — to 
overcome  unintended  spanwise  shifts  in  the  nacelle/diverter  locations  caused  by  changes  at  the  wing  root,  particularly  if  the  fuselage 
diameter  has  been  modified.  Hence  an  artificial  boundary  was  introduced  (at  span  station  125”  for  the  TCA  designs),  and  the 
perturbation  procedure  now  absorbs  all  wing-root  changes  in  the  resulting  near-body  panel. 


Improvements  to  SYN87-SB  Since  7/96  (Continued) 
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- Stage  3:  Ensure  Exact  (x,y,z)  Matching  at  All  Boundaries 

- TCA  Optimization  Made  Use  of  5 Sets  of  Nacelle  Data  (4  Updates) 


SYN87-SB  Gradient  Comparisons  (Adjoint  vs.  Finite  Differences) 
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SYN87-SB  Gradient  Comparisons  (Adjoint  vs.  Finite  Differences,  continued) 
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TCA-6  HSCT,  Wing/Body,  Mach  = 2.4  1 1/21/96 


SYN87-SB  Gradient  Comparisons  (Adjoint  vs.  Finite  Differences,  continued) 
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TCA-6  HSCT,  W/B/N/D,  Mach  = 2.4  1 0/24/96 


SYN87-SB  Gradient  Comparison 

L/D  Objective  for  fixed  Alpha  = 3.6  deg.;  2x10  sine  bumps  modifying  camber 

A Pseudo-nacelles,  forward  differences  (converged  7 orders,  h=0.0001) 
Pseudo-nacelles,  adjoint  method  (converged  6 & 3 orders,  h=0.0001) 
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TCA-6  HSCT,  W/B/N/D,  Mach  = 2.4 


SYN87-SB  Gradient  Comparisons  (Adjoint  vs.  Finite  Differences,  continued) 
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SYN87-SB  Gradient  Comparison 

D/L  Minimization  at  Fixed  Alpha;  17  SIN2  bumps  at  K = 9 modifying  camber;  exponents  7 

I Forward  differences  (converged  6 orders,  h=0.00015) 

X Adjoint  method  (converged  6 & 3 orders,  h=0.0001 5) 


1248 


Ames  12-31,  W/B/N/D,  Mach  = 2.4 


SYN87-SB  Gradient  Comparisons  ( Adjoint  vs.  Finite  Differences,  continued) 

Figure  5 (lower  surface,  but  otherwise  the  same  as  for  Figure  4)  shows  better  agreement  in  the  nacelle  region,  possibly  because  these 
gradient  elements  are  not  as  small. 
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Ames  12-31,  W/B/N/D,  Mach  = 2.4  01/31/97 


SYN87-SB  Gradient  Comparisons  ( Adjoint  vs.  Finite  Differences,  continued) 

Figure  6 is  for  the  fixed  lift  case,  showing  three  configurations  far  from,  part-way-towards,  and  near  to  the  final  Ames  design.  As  for 
the  fixed  Alpha  situation,  the  moderate  discrepancies  are  in  the  nacelle  region  and  they  increase  somewhat  as  those  gradient  elements 
decrease  in  magnitude  with  the  geometry. 
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W/B/N/D  Cases,  6 & 3 Orders  Convergence,  0.00015  Differencing  Intervals,  Mach  = 2.4  02/22/97 


Further  SYN87-SB  Developments 
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revealed  by  the  above  comparisons,  elrorts  should  be  made  to  retest  the  two  alternative  se 
routines  in  the  adjoint  solver  which  are  currently  present  but  not  active  in  the  production  code. 


Further  SYN87-SB  Developments,  continued 
Multitasking 
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Further  SYN87-SB  Developments,  continued 
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Future  SYN87-SB  Developments! 
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Ames  Optimized  TCA  Configuration 
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Team  Effort 

Wing  Design:  Susan  Cliff 
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OVERFLOW  Computations:  Goetz  Klopfer 

David  Baker 
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Design  Approach 

Optimization:  SYN87-SB  (Euler,  single  block)  + NPSOL/QNMDIF 
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Periodic  explicit  smoothing  in  circumferencial 
direction 
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Wing  Defining  Stations  in  SYN87-SB 


C n a 2 

w>  § .2  ^ 

2 o jg  S 

o 73  « S 

•2  g W)^3 

J5  0.S3 

3 ‘3  C *G 

p , 03  W> 

E 2 

° w ? ^ <u 
^ ^ o S ^ 

-C  O _5  o c 

2 g i £ 

■g  J *13  § uj 

O e3  JB  5 '■§  C 

rv  w t3  <3  2 c3 


_ m-h  ■>. 

•§  ° 


£ §^-0  ^ 
S S op  *> 


u n iu  h 

*p£S  i-si 

.2  — - > d 5 ^ 

z'a'^sZ  § 

0.^3  D T3  3 ° 
03  33  <L>  o 

o to  8 H Z * 
££^2o 
Hh  c u .5,  c 


1268 


Wing  Defining  Stations  in  SYN87-SB 

AIRPLANE  Cp's  for  AMES  1-03 
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Optimization  Procedure 
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705  Comparison  of  Lower  Surface  Cp 

The  lower  half  of  this  figure  shows  the  AIRPLANE  lower  surface  grid  and  Cp  contour  li 
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Design  Statistics 

220  Design  changes,  40%  were  determined  unsuccessful 
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used  in  later  designs 
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Euler-Based  Optimization  of  TCA-6 

Delta  L/D  Improvements  Between  Design  Runs 
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SYN87  Calculations  (193  x 49  x 65  Mesh),  Mach  2.4 

Pseudo-Nacelle  Pressures  from  AIRPLANE  01/24/97 


Euler-Based  Optimization  of  TCA-6  (Cumulative  L/D  Improvements  Between  Runs) 

This  figure  displays  a running  sum  of  the  performance  improvements  shown  in  the  previous  figure.  Several  points  of  interest  are  better 
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SYN87  Calculations  (193  x 49  x 65  Mesh) 
Pseudo-Nacelle  Pressures  from  AIRPLANE 
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SYN87  Optimization  of  TCA-6:  Adjoint-Based  Gradients  vs.  Finite  Differencing 

Cumulative  L/D  Improvements  from  Design  Runs 

Successful  Design  Changes  Retained 

• Adjoint-Based  Gradients 


SYN87  Calculations  (193  x 49  x 65  Mesh) 
Pseudo-Nacelle  Pressures  from  AIRPLANE 
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Intermediate  Design  Cps:  AIRPLANE 
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Intermediate  Design  Cps:  AIRPLANE 

Mach  2.4,  alpha  3.756 


C\l 

o 

CO 

o 

CO 

a 

CL 

O) 

O 

‘co 

CD 

<D 

Q 

CO 

T3 

lo 

CL 

c 

u_ 

O 

CO 

CM 

o 

1 

1 

CO 

Q- 

o 

O) 

CO 

CL 

o 

0) 

Cps 

15 

"O 

CL 

3 

15 

T3 

CL 

ID 

CO 

lO 

O 

o 

1 

CM 

r- 

-T— 

O 

CO 

CO 

CO 

CL 

CL 

O 

O 

0> 

O) 

15 

15 

“O 

"O 

CL 

CL 

to 

3 

3 

T- 

co 

OO 

o 

1 

<t 

O 

o 

1— 

1— 

1285 


NASA  Ames  1-03  Wing/Body  (front  view) 
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NASA  Ames  1 -03  Wing/Body 
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Baseline  Wing/Body 
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NASA  Ames  1 -03  Wing/Body 
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Baseline  Wing/Body 
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Wing  Trailing  Edge  Spar  Depth  Constraint 
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Thickness/Chord  Ratio  Constraint 
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Landing  Gear  Box  Constraint 

The  spanwise  cross-sectional  cuts  shown  in  this  figure  were  taken  at  fuselage  stations 
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B2  COCKPIT  AREA  CONSTRAINT 


Fuselage  Curvature  Constraints 
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Fuselage  Curvature  Violations 
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Ames  01-03:  Pressure  Vessel 


Euler  Computations  - Wing/Body/Nacelle/Diverter  - Lift  Curves 
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Euler  Computations,  Wing/Body/Nacelle/Diverter 

M=2.4,  no  internal  or  base  nacelle  forces,  entire  fuselage 


Angle  of  Attack 


Euler  Computations  - Wing/Body/Nacelle/Diverter  - Drag  Polars 


c n 

3 

o 

*> 

8 3 

o-E 

jg'R 

•5  o 


o 


X 
P 

N <7  . 

c9  = 
.§«"  o 


O ^ i3  on 


C 3 

£ CO 

o *- 

X W 

1/5  z 

P < 

3 Ph 


CU  £ 

Dh  c 


C 

3 

ON 

l> 


m 

o 


C O H 

X c/5 


P 

C/3 

3 
g X) 
3 p 
X 


D -C 

|S 

o 3 
0-3 
O cj 


I 

< 

- 2^ 
P X 
c/3  — ' 

Ph  £ 
P 

c/3 

c 

o 


bJQ  § ^ 
3 P X 
,P  cS 


5tS  < S P ^3  g 


c 

§ ^ 
3 *p 
3 3 


3 


W)  P<  j£ 


*g  E 

§ 8 

o 


C X)  +-  X 
^ ^ u T3 

WC?  ^ 
1^1.. 
111 


o 

8 
o 
o 
o 

c/3 

g X 

lc 

3 o 
2 gp 
o-  S 

3jj 
3 


"O 

P 

N 


P 
O 

.s  « 

i o ’S 


3 22  *3  X 


^5 


X 
X 

. <r>  r-  O 


8"9 

•g  JS 


3 3 


X 

c 

3 


p.a 

W).ts  CO  Z 


C/3 

p 

T3 

o < 
o 

o 

£ 


c3 

Dh 

3 


P 

C 

*p 

V5 

3 

X 

P 

X 


£ 8-2 

g X O 


u r-  ^ 

oo  CO 


C/3  CQ 


0) 

€* 

u 3 

X £ 3 

P 'H— ( 

E ° 

3 

05 


o S*> 

§>CQ 


3 CO 

8 ^ X 
'53  J§  00  ‘ ‘ 


o 


W> 

c 

3 


;5oo 


3 X 


O 

C/5 

jj  o 
•5  x 

O C/3 


c CO 

? -0 
C 
3 


>/o  >>i2 

g >,  x 

0-0  0 

c 


C/3 

3 

O 


JJ 

”3 

W 

o 

o 


O OQ 

S.| 
■£  00 


O ^ - 

c ^ .H 

c^VOT3 

jS  ^ 

Si 

I*  < 

&jOi— 3 
«J  fiU 


C/3 


8; 

Ou 


3 


P, 

S* 

C/3 

C 

o 


p 

p ^ 

x ~ 


3 


. 3 

o£ 

C 

"8  « 
P X 
> X 


_ J>  8 _ 

^ nX  M 

e o^-ss 


o 
o o 

T3  ^ 

o OC 
O 

W>  2 

<D  3 


0) 

S.'Q  M 
2 ^ 

|< 

— >>&  ■£  « 

O £> 

00 

X 73 

a)  o 4 

■S  'S 

H P X-*  X C/3 

c &S?  g 2 

§ gl  & 

1 Ph  U—i 

Sf?  p 
P <n 


<u 

O X 


JS 

3 


_ o 

3 X 
JJ 

&jo  3) 


<D 
: T3 


po  w>  2 ~3 

H X T3  > 


3 

Uh 

3 

X C 
C 3 

o o 

o o 


T3 

P 

C/3 

3 

3 

3 

T3 


P 

C/3 

3 

O 


1308 


Euler  Computations,  Wing/Body/Nacelle/Diverter 
M=2.4,  no  internal  or  base  nacelle  forces,  entire  fuselage 


Effect  of  Fuselage  Fairing  - Wing/Body/Nacelle/Diverter 

?Ue  t0  the  aerodynamic  fairing  on  the  side  of  the  fuselage  near  the  trailin' 
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Effect  of  Fuselage  Fairing,  Wing/Body/Nacelle/Diverte 
M=2.4,  no  internal  or  base  nacelle  forces,  entire  fuselage 

* AIRPLANE  1 -03  - Pressure  Vessel  Fuselage 
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Euler  vs.  Navier-Stokes  - Wing/Body/Nacelle/Diverter  - Lift  Curves 
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Lift  Curve  Comparisons,  Wing/Body/Nacelle/Diverter 
M=2.4,  no  internal  or  base  nacelle  forces,  entire  fuselage 
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OVERFLOW  Computations,  Wing/Body/Nacelle/Diverter 
M=2.4,  no  internal  or  base  nacelle  forces,  entire  fuselage 
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Pressure  Drag  Comparisons,  Wing/Body/Nacelle/Diverter 
M=2.4,  no  internal  or  base  nacelle  forces,  entire  fuselage 
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Baseline  TCA:  AIRPLANE  Solution 
Mach  2.4  a -3.6 
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Ames  1-03:  AIRPLANE  Solution 
Mach  = 2.4  a = 3.756 
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Baseline  TCA:  AIRPLANE  Solution 
Mach  = 2.4  a = 3.6 


1327 


.v.»«:npv  j |0|imi  ) ..j  i <1  j|  .1  ijikii  i';  uoi|i>ui  ■ n 11  | 

of  0 sro  000 


Surface  Pressures  - Ames  1-03  - AIRPLANE 


00  > 

§ 3 

.2  o <D 

> c/5  X X 

a.  g,-a  » 


X 


s 

o 

c if 

•—i  X 

<l>  H3 

a & 


a3  cd  <3 

^ U<  <u 

P x 


3 

X) 


„ S>  a 

x S 2? 
g<x  x 

8 § § 
c s g< 

Oh  <d 

iD  c/3 

«->  X 

3 


a -8 
O o 

ui  a 
<D 


C/3 


OJ  _ 

’S  o 

2 2 -S  a .g 

c3  *3  ^ "g 

j-i  .g  ^ G C 
uh  H-h  X'  3 3 

3 ^ x ->.22 

8 2 — 73 
£ o D 
O »j  X 

O *c. 

o ^ c 

<U  g *- 1 £h 

g .5?  ° a.  B- 

§ S g 
1/5  -a  «3 
g <u  <u 

c 53  § 


C/3 

3 

O 


8 8 


>%  D 

"S  3 


o 

3 


<D  ex!  X X 
D C ^ <U  3 

g E'g’g  c 

I > 1 3 

■S'p  Mu  S 


» O S > M 

8 £.2  J.S 

5 ^ 3 

< c r 

^ 5 32  D ^ 

^ <u 


<u 

x 

» 

U< 

,o 

C/3 

<u 


3 


3 

O 


3 <u 

C/3  Crt 


3 3 

CX,  <U 

B v 

o 

CJ 


o D 5 


3 g 

8 | 

D 

c/3  D s— ' 

B £ g 

3 O £ 


C/3 


<U 


<2  §m<2,m 

» — < r/“i  ’ '"O  - ^ Cj__i 


<U 

Ut 

3 


5 -2i 


C/3 

3 C/3 


^ o3  U o * 


x> 

& 

o 

W 

z 

< 

iJ 

cu 

Pi 


cxx  « 

8 

o . *- 

O p d 
('  ■ fa  X 

S)H 


3 C/3 

<U  3 
O 


i- 

£ 

<D 
c /D 

o 

cx 

G 

3 

3, 


•2  2 >- 

< g g§>’§ 

<U 

X 


3.X 


2o« 


a .5 


L .h  X v/ 
H ^ X o 


5-H 

CX 


1328 


Ames  1 —03:  AIRPLANE  Solution 
Mach  = 2.4  a = 3.756 
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M = 2.4,  CL  = 0.1 
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M = 2.4,  CL  = 0.1 
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Approach  for  Complex  Geometries 


1354 


The  current  optimization  efforts  employ  the  Euler  equations  to  represent  the  flow-field.  The  following  is  a standard  presentation  of  the 
Euler  equations.  The  AIAA  paper  97-0103  (Constrained  Multipoint  Aerodynamic  Shape  Optimization  Using  an  Adjoint  Formulation 
and  Parallel  Computers)  and  associated  references  give  a complete  treatment  of  the  solution  algorithm. 
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Multiblock  Euler  Flow  Solver,  Continued 
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1358 


Supersonic  Analysis  of  Three  Optimized  Designs 

Jh/oughout  the  recently  concluded  optimization  efforts  by  NASA  Ames  on  the  TCA-6  geometry,  the  multiblock  Euler  flow  solver 
(rL087-MB)  was  successfully  applied  to  the  baseline  and  intermediate  designs.  The  full  configuration  consisting  of  the  wine 
f^age  nacelles,  and  diverters  was  modeled  using  a mesh  with  180  blocks  and  1.5  million  cells.  The  drag  increments  predicted  with 
FL087-MB  were  consistent  with  the  results  from  the  AIRPLANE  code. 
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MDC  Cycle2:  SYN87-MB 
Mach  = 2.4  C.  = 0.1 
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Supersonic  Analysis  of  Optimized  Designs,  Continued 

MDC  Lower  Surface  Iso-Cps 

The  lower  surface  iso-Cp  contour  plots  for  the  MDC  design  continue  to  display  the  smoothness  of  their  pressure  contours.  The 
severity  of  the  shock  interactions  in  the  nacelle  regions  have  been  ameliorated  slightly. 
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MDC  Cycle2:  SYN87-MB 
Mach  = 2.4  C = 0.1 
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Boeing  WB8WT:  SYN87-MB 
Mach  = 2.4  C = 0.1 
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Supersonic  Analysis  of  Optimized  Designs,  Continued 
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Boeing  WB8WT:  SYN87-MB 
Mach  = 2.4  C = 0.1 


1370 


Supersonic  Analysis  of  Optimized  Designs,  Continued 
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Supersonic  Analysis  of  Optimized  Designs,  Continued 
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Supersonic  Analysis  of  Optimized  Designs,  Continued 
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Lift-Drag  Polar  Comparison 
Wing-Body-Nacelle-Diverter  (No  Interior  Forces)  SYN87-MB 
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MultiblockAdjoint  Method  for  the  Euler  Equations  (ADJ87^ 
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Multiplying  by  a Lagrange  Multiplier  ip  and  integrating  by  parts  gives 
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The  cost  function  thus  becomes 
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Multipoint  Adjoint  Method,  Continued 
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The  Adjoint  Solver  Uses  the  Same  Multiblock  Strategy  as  the  Flow  Solver 
• Discretization  and  Time  SteDnin^ 
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Coarser  meshes  in  the  multigrid  cycle  have  the  same  double  halo 
structure  (Permits  proper  multigrid  convergence  acceleration) 


Geometry  Modification  and  Constrained  Optimization 
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Geometry  Modification  and  Constrained  Optimization,  Continued 

The  imposition  of  linear  and  nonlinear  constraints  in  the  optimization  procedure  has  been  handled  by  coupling  the  multiblock 
approach  to  the  NPSOL  algorithm  of  Gill,  Murray,  Saunders,  and  Wright.  NPSOL  is  a sequential  quadratic  programming  (SQP) 
method  in  which  the  search  direction  is  calculated  by  solving  the  quadratic  subproblem  where  the  Hessian  is  defined  by  a quasi- 
Newton  approximation  of  an  augmented  Lagrangian  merit  function.  The  Langrange  multipliers  in  this  merit  function  serve  to  scale  the 
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Multiblock  Mesh  Perturbation  Algorithm  (WARP3D-MB) 
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Multiblock  Mesh  Perturbation  Algorithm  (WARP3D-MB) 
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For  a more  detailed  explanation  of  the  3-stage  WARP3D  grid  perturbation  scheme,  see  AIAA  paper  96-0094  (Aerodynamic  Shape 
Optimization  of  Complex  Aircraft  Configurations  via  an  Adjoint  Formulation),  1996. 
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• WARP3D 
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Parallel  Implementation 

The  main  strategies  used  to  accomplish  the  parallelization  of  the  design  code  are:  a domain  decomposition  model,  a Single  Program 
Multiple  Data  (SPMD)  strategy,  and  the  Message  Passing  Interface  (MPI)  standard  for  message  passing.  The  choice  of  MPI  was 
determined  by  the  requirement  that  the  resulting  code  be  portable  to  different  computing  platforms  as  well  as  to  homogeneous  and 
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Parallel  Implementation 


U 


For  the  Application  of  Design  Variables  and  Mesh  Perturbation  No 
Message  Passing  is  Needed 
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Multipoint  Constrained  Transonic  Design 
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- Design  Point  3:  M = 0.83,  CL  = 0.25 


Multipoint  Constrained  Transonic  Design,  Continued 
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Transonic  Business  Jet  Configuration 
Iso-Cp  Contours 
M ss  0.82,  Cj  = 0.30 


Optimized 
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Business  Jet  Configuration  Configuration.  Drag  Minimization  at  Fixed  Lift. 

M = 0.82,  CL  = 0.3 

90  Hicks-Henne  variables.  Spar  Constraints  Active. 

— , Initial  Pressures 
— , Pressures  After  5 Design  Cycles. 


1400 


XJC  w 

c:  span  station  z — 0.665  d:  span  station  z — 0.856 


Business  Jet  Configuration.  Multipoint  Drag  Minimization  at  Fixed  Lift. 
Design  Point  1,  M = 0.81,  Cl  = 0.35 
90  Hicks-Henne  variables.  Spar  Constraints  Active. 

- - Initial  Pressures 
— , Pressures  After  5 Design  Cycles. 
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c:  span  station  z = 0.665 


d:  span  station  z = 0.856 


Business  Jet  Configuration.  Multipoint  Drag  Minimization  at  Fixed  Lift. 
Design  Point  3,  M = 0.83,  Cl  = 0.25 
90  Hicks-Henne  variables.  Spar  Constraints  Active. 

— , Initial  Pressures 
— , Pressures  After  5 Design  Cycles. 
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Business  Jet  Configuration.  Multipoint  Drag  Minimization  at  Fixed  Lift. 
Design  Point  1 , M = 0.81,  CL  = 0.35 
90  Hicks-Henne  variables.  Spar  Constraints  Active. 

- - Pressures  for  Design  Test  Case  1 . 

— , Pressures  After  5 Design  Cycles. 
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c:  span  station  z = 0.665  d:  span  station  z = 0.856 


Business  Jet  Configuration.  Multipoint  Drag  Minimization  at  Fixed  Lift. 
Design  Point  2,  M = 0.82,  Ch  = 0.30 
90  Hicks-Henne  variables.  Spar  Constraints  Active. 

— , Pressures  for  Design  Test  Case  1. 

— , Pressures  After  5 Design  Cycles. 
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d:  span  station  z = 0.856 


Business  Jet  Configuration.  Multipoint  Drag  Minimization  at  Fixed  Lift. 
Design  Point  3,  M = 0.83,  CL  = 0.25 
90  Hieks-Henne  variables.  Spar  Constraints  Active. 

- - Pressures  for  Design  Test  Case  1. 

— , Pressures  After  5 Design  Cycles. 
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Constrained  Supersonic  Design,  Continued 
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Design  Variables 

“ 144  Hicks-Henne  Functions  at  8 of  1 
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Constrained  Supersonic  Design,  Continued 
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Supersonic  Transport  Configuration 

Iso-Cp  Contours 
M = 2.2,  CL=  0.105 

Upper  Surface  Lower  Surface 


Conclusions 

This  presentation  has  shown  how  the  complicated  design  of  both  transonic  and  supersonic  aircraft  configurations  including 
airframe/nacelle  integration  effects  can  be  accomplished  in  a routine  fashion  for  multiple  design  points  with  the  inclusion  of 
constraints  utilizing  a multiblock  adjoint-based  design  method. 

This  mrthnHnWv  is  advantageous  in  two  respects.  First,  the  utilization  of  a multiblock  scheme  allows  the  treatment  of  complex 
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Conclusions 
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Unstructured  Meshes 

Constrained  Multipoint  Optimization  is  Demonstrated  with  a View 
Towards  Eventual  MDO  Problems 
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Currently  Available  Aerodynamic  Shape  Design  Methods 
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CDISC  Langley  Euler/NS  82-Pres.  Full  Actual  Inverse  Type  Method 
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Developing  Navier-Stokes  Flow  Solver  Technologies 
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• J Coarsening  Multigridding 

• Block  Jacobi  Preconditioning 
Turbulence  Models  ??? 


*-«  i-*  <D  b£  co 
n ^ is  it/ii 


^ o 

X SISK'S  CO  £ oi_e  H 


S c 

c o 


isi 

is-2 


(1)  (L)  c«  <D  5“H  Ol, 

« SP*  y & ac 

u. 


as 


-2 

3 


<D  ^ C_  o 
cfl-  . rt 
J)X1  JJ'O 

^ i wa 


co  co  OJD^c 


•2  3 n o 


cd  *0 


c J—  3 ^ 
« 3 ^ ^ 

2>  t4_H  X X 

rrt  I —* 

O 

u as 


rt 

E 


3 S.2  « 3 td  3 g 


« s 

4>  X) 

J*J  - 
a3 


’O 

C v_ 

as  < 

W)o" 


o co> 


CO 

E 2 D >< 


cd  o 
CJ 


1/5  ,n 

S « «5 

C cd  .£ 

4)  _ 

bX>3 


os  S /2  ^ -2 

^ 3 *3  ,S  ”§ 

2 oa£2 

on  bS  T3  £ £5  *> 


£ D 

SoSHu 


rrt  03 
«.v  ^ 

^ x 2 o 

i!  Ot2  £ 

^ *-*  0^  <L>  ~ 

cfl’O  O^ 

.2  8 S ^ 

5P  3 .2 

-P  co  £ 


£ <D 


o " 

.§■§ 

»c 


11 
11 
CO  ^ 

c is 

5 -2 

p 3 


o ^ ^ 

W)  c * 
3^o 
a 2 
2 S 

X u 


3 

3 

a- 

co 

a> 

x 


co  ”3 

«* 

« Su 

^ £ ~±i 

Cd  p 3 *- 


o*x 

3 o 


CO  U 

E 


& 

S?X  CO  4- 

3 *w 

g<o  gx 
*2^x  2 

><'i  S « 

•8  S2P<  ” 

o x 3 

CO  CO  an 

a aJ= 


cd 
O 
X!  D 

3 +5 


E *o 

S a 

o 4> 

c 

P 
o 


£ S|S 

°il| 

co  o b ^ 

O 

Scd  T3 

g CJ 

co  ^ __,  •> 

.a  c .g  c £ 

o 2** 


* £ 


*—  — 

a3  X 


O 
O 

,E3 

u «? 
o 3 
o 


g w2 

^ -i  J—  o 

I 8-S 

§ § ^^^5“ 
« §•- 
■g  I ®p.2f>  H £9 

. m r-1  C w 

»S1 

c *C  c 

O bow 

*4— > * J— 

■s-i.2  2 

% E > x 

b '*-'  "2 

S-o  S 53 

« C-r 
O a3  — 
co 


Is 


M.S2  8 22 


'S 


P 

So 

8 o 

“ Eo 


-S  o 


CO 

co 
<D 
CJ 

O 

U(  V3 

D-  co 

O 

o 
o 


cSO 
,0-0 
O-  2 
% oc  C 

^ ' O 

oo  E 
2 -H 
S- 

c o e 
So- 

> OO 


iu  — > - c 

Oh  W^T* 


CO 


3-g^CO  ctJ 

g 3 3 t xs 
•-  i5  <N  7- 
3 Cu  4^ 


>,  5 
c > 

as  & 
o J> 

co  -C 
3 ‘ 


3 w-t  -a 


s “ M S 1 

S 8.s  ^ s> 


co  C 

a> 

W).t5 

r C 4— * co 
^ w X u 


W>^ 
c 2 
o -g 
E 

b co 
as  v 

co  a 

o 


§s 

'"m 


(D 

O 

3 

3 

3 

-e 

3 


E 2 


O •- 
C X 

,ss^  8i? 
C!S  £ 

^ — -co 


a>3  n 
S £ "3  u-° 


o 


co 


cS 


o 


<D 


i x3 


as 


o o -2  ■£ 

S^| 

wn£E“'3_ 
|oS«^2og- 
^04^  .xspj-iJx: 


(D 

- £ z 

aS  ”0  . 4—> 

m (U  1)  a) 


ct£ 

4>  Mh 


og,^e^cS£a-s£--o-oO 
■Siooo.2,a"1§,s,»2->f  ? 


t«  u is 


S 8-S 

O | ^ 3 

D ^3 

><  P 


'O  ^i— » w 

1 O 

S 2 


> 

— 3 

CO 

c o 
o w 
o c 


hJ  -5 

Ph  tt  - 

co  _ 
u-  a)  3 

>'3>^ 
o5  SJ 

CO  O *Z3 

Dh  as 

O ~ 


„ *o 

^5  -2  § 

S)  g*T> 

C 1/3  C -C 

3 4)  CN 


Oh  C 

E P 


£ 

o 

C 


o 0^  £ T3  X5 

■ « Qj  

CX  2r^  . "o  •“ 


0.52 

'S  E 

>».«  M 
in  g "O  ■ n 3 jz 

ae  8 
f « 2 

g t«  >>  CO 
co  > (j  X 

0^0“ 


2f)  o ^ 

^ OX)  o ^ 
co  cd  t3S 
S to  ^3 
D 


X 

O 
co 

OX) 
_ C 

oig 

3 3 

t. 

to  OX) 
C * — 

o 


fl  ' 

O, 


o 

o 


- ^3 

-e  s- 

cd 


2 4> 


w ^ o i s c 

3 3 c^’oi 


OX)X 

c ^ 


as 

‘5bij  3 

■-  <u  3 

o o 

o -o 

cd 


CO 


S3 

^ £ <o  >. 

8 sf  z 

2 ^ w T3 
jS  o O C 


ci  o 
o 


_ -o  _ 

© o .S3  4)  -g  5 (U 

o — • cd  <U  J3 

bX)X  rr  bX)t3 

f § 5 “ 

0^-0  as 

Oh 

aT  P P 


o> 

aS  <4 
co 


H 

.11 

e 3 

£.22 

u. 

3 U« 

o o 

o X 


s 

I- 

Oh— « 

S § 

o .£ 

> o 

0) 


a g 


C 
3 

E 
E -o 
o 
o 


X 
as 

E 

. o 

22  i-J 

ii 

co  ^ 
P^X3 
W aS 
to  PQ 

§ « 

!h 


bx>< 

C 

X X 
3 0) 

Sn 

Ijj 


o <u 

2 x 

04H 

P • 

04^ 

X 


TJ  p 

uh  ro 
0 34 
O4  3 

a g 

3 P 


co  c 


CO 

as  *S 


C 

cj  2 

2 B 

aS  p 

^-2 

3 ° 
p'S^ 

22  -8 


U 


5 x x 

sX  Cd  CO 


c g 
X § 


o O s 2 3 s 

co  X as  co  X cd  X1 


a) 

> _ 
O 03 


s-g 


« O 

t-H  X 


= <u 

Q> 

ll-SS-2 
»§jl| 

°S;'S-c  3 

OO  > X 


Oh  £ E 
O 4> 


co  0 

Is 

co  o 

CO  ^ E 


1427 


FL01Q7-MB  Multiblock  Navier-Stokes  Solver 
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- Residual  Smoothing 

- Full  Multigridding 

- Preconditioned  J-Coarsening  Multigridding  (under  development) 


Analysis  Carried  Out  on  Parallel  Computers 

- Code  ported  to  IBM-SP2,  SGI-R8000/R10000,  CRAY-J90,  and  HP 
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Computational  Costs 
- Up  to  10  Analysis 


Low-Wing  Transport  Configuration:  Trailing  Edge 

Navier  Stokes  Solution 

Mach  = 0.86  CL=  0.6  Re  = 50xl06 

1.5  Million  Mesh  Points  48  Blocks 
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Computational  Costs 

- Each  Analysis  was  Run  for  500  Multigrid  Cycles 

- The  Average  Residual  was  Converged  4 Orders  of  Magnitude 
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As  a final  example  of  the  analysis  capability  of  FLO107-MB,  a complete  configuration  transonic  business  jet  is  analyzed.  The  mesh 
took  1.5  weeks  of  interactive  time  using  GRIDGEN.  The  objective  of  this  study  was  to  determine  the  limits  of  current  N-S  capability 
to  extend  into  the  transonic  buffet  regime  and  to  predict  drag  rise  accurately.  The  configuration  included  wing,  body,  nacelle,  pylon, 
vertical  tail  and  horizontal  tail.  The  mesh  featured  240  blocks  and  5.8  million  computational  cells.  Each  analysis  was  run  for  300 
multigrid  cycles  on  32  processors.  The  average  residual  was  converged  4.3  orders  of  magnitude  and  took  3.25  hours  of  wall  clock 
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Transonic  Business  Jet  Configuration 
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- 1.5  Weeks  Initial  Mesh  Generation  Time 


Computational  Costs 

- Each  Analysis  was  Run  for  300  Multigrid  Cycles 

- The  Average  Residual  was  Converged  4.3  Orders  of  Magnitude 
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Transonic  Business  Jet 
FLO107-MB  Solution 
Baldwin-Lomax  Turbulence  Model 
Mach  = .82  240  Blocks  5.8  Million  Mesh 
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all  of  both  the  inviscid  and  viscous  analysis  methods  applied  to  the  TCA-6  configuration  have  agreed  well  on  the  pressure  drag. 
However,  when  the  optimized  designs  are  analyzed,  the  inviscid  and  viscous  methods  show  discrepancies.  In  the  case  of  viscous  drag, 
the  discrepancies  are  considerably  greater,  and  they  are  only  clouded  further  by  the  trip  drag  and  Alpha  uncertainty  issues  apparent  in 
all  wind  tunnel  tests.  It  is  well  within  the  HSR  charter  to  seek  resolution  of  these  issues. 


Design  Method  Requirements] 
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Finally,  as  seen  from  a practical  point  of  view,  a constructive  analytical  study  of  relevant  physical  phenomena  may  lead  to  the 
formulation  of  objective  functions  that  are  better  suited  to  supersonic  design.  It  is  fortunate  that  for  transonic  flows  an  intuitive 
understanding  of  the  shape  of  the  pressure  distribution  can  lead  to  the  development  of  very  efficient  designs.  While  the  Natural  Flow 
Wing  technology  takes  a stab  at  this  type  of  approach  for  supersonic  configurations,  it  would  be  beneficial  if  a far  more 
comprehensive  understanding  of  the  flow  physics  could  be  established.  Such  an  understanding  may  point  to  new  objective  functions 
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* Nacelle  Shape  Design 

* Structures  and  Materials? 
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In  a view  towards  making  the  design  tools  both  current  and  future  to  be  more  versatile  and  easy  to  use,  improvements  are  needed  in 
the  areas  of  geometry  treatment,  mesh  generation,  and  mesh  perturbation. 
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Geometry  and  Mesh  Requirements 
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— Offers  Greater  Geometric  Complexity 

- Reduces  Required  Number  of  Mesh  Points 

“ Automatic  Mesh  Generation/Perturbation  Seems  Feasible 


One  apparent  realization  that  is  clear  about  all  of  the  design  optimization  methods  currently  available  for  HSR  applications  is  that  they 
are  not  easy  to  use  by  anyone  except  the  developer(s)  or  those  that  are  closely  associated  with  the  developer(s).  To  change  this  trend,  it 
will  be  necessary  to  build  a tool  package  that  interfaces  better  with  these  sophisticated  modem  design  methods.  It  is  ironic  that  we  are 
living  in  the  era  of  the  information  revolution  yet  our  best  aerodynamic  design  tools  fade  further  and  further  into  user-unfriendliness 
owing  to  their  extreme  complexity.  Four  areas  needing  attention  are:  design  space  visualization,  constraint  space  visualization,  a 
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Visualization  Techniques  and  Tools  for  Engineers 
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AN  ANALYSIS  OF  CFD  AND  FLAT  PLATF  PREDICTIONS  OP 
FRICTION  DRAG  FOR  THE  TCA  W/B  @ SUPPRsAnI^ 

Scott  L.  Lawrence 
NASA  Ames  Research  Center 
Moffett  Field,  California 
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This  paper  presents  results  of  a study  which  attempted  to  provide  some  understanding  of  the 
relatonship  between  skin  friction  drag  estimates  produced  by  flat  plate  methods  and  those  produced 
by  Navier-Stokes  computations.  A brief  introduction  is  followed  by  analysis,  including  a fiat  plate  qrid 
study,  analysis  of  the  wing  flow,  an  analysis  of  the  fuselage  flow.  Other  results  of  interest  are  then 
presented,  including  turbulence  model  sensitivities,  and  brief  analysis  of  other  configurations 
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SUMMARY  OF  TCA  COMPONENT  SKIN  FRICTION  DRAGS 


Wind  Tunnel  Conditions:  Re  = 4 x 106/ft,  1.675%  scale  model 

0.008 
0.007 
0.006 

CDF 

0.005 
0.004 
0.003 
0.002 
0.001 
0.000 


Integrated  Skin  Friction  Drags 


OBJECTIVES: 

-UNDERSTAND  WHY  N-S  METHODS  TEND  TO  UNDERPREDICT  FLAT 
PLATE  CDvON  THE  WING 

-DETERMINE  IF  THERE  IS  A “BEST’  FLAT  PLATE  METHOD(OLOGY) 


Lawrence- HSR/CA-  - 3 of  5$” 
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This  study  was  primarily  motivated  by  two  existing  situations:  1)  the  substantial  variation  in  flat  plate 
predictions  produced  by  the  various  organizations  involved  in  the  optimization  effort,  and  2)  the 
substantial  variation  between  typical  Navier-Stokes  predictions  for  friction  drag  and  flat  plate 
predictions. 

This  figure  illustrates  both  of  these  issues.  Friction  drag  is  plotted  for  six  different  prediction  methods, 
broken  down  by  component.  The  flat  plate  methods  were  applied  consistently  within  the  NASA  Ames 
friction  drag  code  developed  by  Hicks,  and  may  not  represent  the  identical  method  of  application  used 
by  the  associated  organizations.  Friction  drag  is  computed  for  the  fuselage  with  a single  evaluation  of 
average  skin  friction  for  the  Reynolds  number  based  on  the  length  of  the  vehicle.  The  wing  is  broken 
spanwise  into  strips  and  average  skin  friction  is  evaluated  on  each  strip  using  the  Reynolds  number 
based  on  the  local  chord  length.  The  flat  plate  wing  and  fuselage  drag  estimates  have  been  adjusted 
so  that  the  associated  areas  agree  with  the  areas  as  demarcated  in  the  CFD  grid. 

A four  drag  count  discrepancy  is  observed  between  the  highest  (van  Driest  II)  and  lowest  (Prandtl)  flat 
plate  predictions  for  the  wing/body  total  friction  drag.  On  a percentage  basis,  the  discrepancy  is 
greater  on  the  fuselage  than  on  the  wing.  This  is  caused  by  the  failure  of  the  Prandtl  method  at  higher 
Reynolds  numbers  as  will  be  seen  in  a subsequent  figure. 

Results  of  th9  parabolized  Navier-Stokes  code  UPS  and  the  OVERFLOW  Navier-Stokes  code  are 
observed  to  fall  within  the  spread  of  flat  plate  methods  for  the  fuselage,  but  are  2-5  counts  below  the 
flat  plate  results  on  the  wing.  The  main  objective  of  this  study  was  to  understand  this  discrepancy.  A 
secondary  objective  was  to  shed  light  on  the  spread  in  flat  plate  predictions  and  obtain  information 
which  might  lead  toward  flat  plate  method  standardization. 
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APPROACH 


• IDENTIFY'  FLAT  PLATE  METHOD  THAT  AGREES  WITH  GRID- 
INDEPENDENT  CFD  SOLUTIONS  ON  A FLAT  PL  ATE 

-TURBULENT  INCOMPRESSIBLE  FORMULAS 
-COMPRESSIBILITY  CORRECTIONS 
-FLAT  PLATE  GRID  STUDY 

•APPLY  THE  CFD-COMPATIBLE  FLAT  PLATE  FORMULA  TO  THE  TCA 
WING/BODY 

-LINE  INTEGRAL 

-SURFACE  INTEGRAL 

•ATTEMPT  TO  ACCOUNT  FOR  DIFFERENCES  BY  QUANTIFYING 
VARIOUS  NON-FLAT-PLATE  PHYSICAL  EFFECTS 

-WING 

-FUSELAGE 


LaWrence-HSH/CA  - b of  55  ” “ — — Appiiad  Computational,  Aerodynamics  Branch 


Results  ofthe  UPS  parabolized  Navier-Stokes  code  were  the  focus  of  the  present  analysis  because  of 
the  authors  familiarity  with  the  code.  Similar  analysis  can  obviously  be  applied  to  any  N-S  solution. 

The  first  step  in  the  analysis  was  to  try  to  remove  what  might  be  called  bias,  i.e.,  differences  in 
computed  and  flat  plate  skin  friction  on  a case  for  which  the  flat  plate  assumptions  are  valid.  A wealth 
of  flat  plate  methods  are  available  from  the  various  incompressible  formula/compressibility  correction 
permutations.  A method  was  sought  that  would  provide  the  best  agreement  with  UPS-generated  skin 
friction  coefficients  for  Reynolds  numbers  ranging  from  wind  tunnel  conditions  to  flight  conditions.  A 
grid  study  was  performed  so  that  the  comparisons  would  be  performed  on  grid-independent  solutions. 

The  UPS-compatible  method  could  then  be  applied  to  the  TCA6  wing/body  using  both  the  line- 
integration  method  described  earlier  as  well  as  by  surface  integration  of  the  local  skin  friction 
coefficients.  The  latter  method  allowed  for  point-to-point  comparison  with  UPS  results. 

Finally,  use  of  the  surface  integration  approach  allowed  for  some  modification  of  the  flat  plate  method 
to  include  some  non-flat-plate  effects,  and  this  allowed  these  effects  to  be  quantified  to  some  extent. 
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TURBULENT  FLAT  PLATE  FORMULAS 

INCOMPRESSIBLE  FORMULAS: 


Method 

Local  Skin  Friction 

Average  Skin  Friction 

Prandtl 

df  - 0.058 /Rex02 

dF  = 0.074 /Rex°2 

Prandtl- 

Schlichting 

C?f  = 0.39/(log/?^)15S 

&F  = 0.455 /( log /f^)2*58 

Karman- 

Schoenherr 

-L  = 4.15]og(C‘,*eJ)+  1.7 

•4=  = 4.131og(4*0 
& 

White 

df  = 0.455/(ln(0.06/?fJ))2 

dF  = 0.523/(ln(0.06/?tfx))2  ? 

Sivells-Payne 

df  = 0.088Iog(/?<ri/233.7)/(log(/?c/31.62))3 

dF  = 0.088/(log(i?ex/31.62))2 

NOTES: 

-KARMAN-SCHOENHERR  IS  THE  ONLY  ITERATIVE  APPROACH 
-THE  WHITE  CF  FORMULA  DOES  NOT  APPEAR  TO  BE  CONSISTENT 
WITH  THE  C,  FORMULA 


Lawrence-HSWCA  - 1 ot  5T" 
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Compressible  flat  plate  formulas  tend  to  consist  of  an  incompressible  formula  that  is  stretched  and 
scaled  with  compressibility  effects.  The  incompressible  formula,  to  a large  extent,  provides  the  basic 
shape  if  the  skin  friction  distribution.  Some  of  more  well  known  of  the  incompressible  formulas  are 
shown  in  this  figure,  along  with  a method  which  may  not  be  as  well  known  (Sivells-Payne)  but  which 
will  be  relevant  to  this  study.  Each  method  for  local  skin  friction,  Cf,  can  be  integrated  from  the  leading 
edge  of  the  plate  to  x to  provide  a function  for  the  average  skin  friction,  Cp 

The  Sommer-Short  and  van  Driest  II  compressibility  corrections  are  typically,  though  not  necessarily, 
associated  with  the  Karman-Schoenherr  incompressible  formula,  the  only  formula  requiring  iteration. 
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T£®  ln^mpreSSible  forrTMJ,as  described  on  the  previous  page  are  plotted  on  the  left  in  terms  of  local 

S ZlTr?  COefhaT-  °"  ^ rightK  formulas  are  vvfth  the  KarmlnSKShJrJ fcSa 

h ,nnLr  r f P®rcenta9®  vanabon.  Not  surprisingly,  the  formulas  compare  fairly  well  in  the  “wind 

yn°  ?S  n.U^be/  range  The  White  formula  is  higher  than  the  others  for  all  Reynolds  numbers 
The  other  most  notable  feature  of  these  plots  is  the  deterioration  of  the  Prandtl  formula  at  high  Re 
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COMPRESSIBLE  FLAT  PLATE  FORMULAS 

STRETCH  AND  SCALE  INCOMPRESSIBLE  FORMULA  BY: 

CcARex)  = ±rC){FxRex)  AND  (fF(Rex)  = jrC‘F(FxRex) 

t c r c 


SOME  RELEVANT  APPROACHES  (FOR  ADIABATIC  WALLS): 


Uwrence-HSTVCA  -TTofbS’ 


Method 

Fx 

Fc 

Nicolai 

i 

(1+0.144  Af/)0'65 

Sommer- Short 

KFJt)-re 

1 + 0.1 157  A/,2 

Spalding-Chi 

(TS  Tf™ 

. (Tw/Te-  lyf 

7 ;/rf  J 

Van  Driest  II 

V(Tt) 

tw/t4~  1 

MTw)Fe 

. /'TW/T4~lSl 

“l  ) 

White-Christoff 

HTt) 

TJT'-  1 

MTJ- J(TW/T,)FC 

“l  tw/t,  } 

Aoptiad  Computational;  Aarodyrtamic*  Branch 


As  stated  previously,  compressible  skin  friction  formulas  are  generally  obtained  through  modification  of 
an  incompressible  formula.  The  incompressible  formula  can  be  said  to  be  stretched  (evaluated  at  a 
scaled-down  Reynolds  number)  and  scaled  according  to  functions  of  the  edge  Mach  number  and  wall 
temperature  ratio.  These  functions  are  shown  in  this  figure  for  some  of  the  more  well-known 
compressibility  corrections.  The  formulas  shown  here  have  been  simplified  by  assuming  an  adiabatic 
wall. 

The  method  given  in  Nicolai  (Fundamentals  of  Aircraft  Design)  is  applied  to  the  Prandtl  incompressible 
formula  to  obtain  the  results  ascribed  to  NASA  Ames. 
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J™™LTnr*  **«*"  were  Perfo^ed  to  insure  that  a reasonably  grid  independent  solutions 
flat  P|ate  method.  The  first  was  performed  for  moderate  Reynolds  number  and 
tne  second  for  high  Reynolds  number.  The  wind  tunnel  study  is  shown  here. 

hSfnS  wer«°^ed  on  1®  different  grids  using  both  the  Baldwin-Lomax  and  Spalart-Allmaras 

f F°T  r f °1 9ridS  distinguished  by  *»  wall-normal  stretching  function  used 
were  studied,  with  four  levels  of  refinement  within  each.  Then  within  each  family,  starting  with  a grid  of 

iTni°min?in?rd  3 "S*  °!  y+a* the  wal1  of  025 ' **  wa"  sF»®ng  was  doubled  at  each  level  of 
f 5 T L°,  ^intS  W8S  Cut  in  half‘  OVERFLOW  solutions  were  also  computed  on 
three  gnds  using  both  central-differencing  and  upwinding  for  the  inviscid  fluxes. 

Skin  friction  coefficients  computed  at  a Reynolds  number  of  10  million  are  plotted  in  terms  of  percent 
variation  from  the  flat  plate  method  of  Sivells-Payne/Spalding-Chi  (SP/SC),  against  four  meshP 
characteristics:  1)  wall-spacing  in  terms  of  y+,  2)  spacing  at  the  boundary-layer  edge  in  terms  of  at  v+ 

^ 4)  *"  ’™*9‘  »*  «““*>«  «*"  •» 

Some  interesting  trends  are  evident.  Rrst,  the  results  tend  to  asymptote  from  above  the  qrid- 

I8000.!’ the  S:A  m^del  generates  consistently  higher  skin  friction  (3-4%)  than 
Slw  HPQ  r ’ l ('S  Rey^in,umber-  Third,  OVERFLOW  results  tend  to  fall  approximately  1% 

! m f eaf1,Lurbulence  model-  Rnal|y- the  UPS  skin  friction  predictions  appear  to 

mcrease  quadrahcally  with  the  average  grid  spacing  in  the  boundary  layer.  This  may  be  afunction  of 

^ereafovP^pTX^U?J^  °VERFLOW  “P™nd  results  (higher)  show  a similar  behavior, 
whereas  OVERFLOW  central-differencing  results  (lower)  appear  somewhat  less  sensitive. 
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The  equivalent  study  at  a flight-like  Reynolds  number  is  shown  here.  A similar  trend  with  average  grid 
spacing  is  observed  here  to  a greater  extent.  Note  that  the  skin  friction  coefficient  is  virtually 
converged  with  200  points  in  the  boundary  layer. 
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Experimental  results  as  compiled  by  Hopkins  and  Inouye  are  plotted  on  the  left  in  generalized  form, 
i.e.,  unstretched  and  unsealed  using  three  different  compressibility  approaches  into  incompressible 
curves.  The  Karman-Schoenherr  formula  is  plotted  for  comparison.  The  data  consist  of  adiabatic  flat 
plate  results  for  freestream  Mach  numbers  ranging  from  1.5  to  5.8  with  the  bulk  taken  between  Mach 
1 .5  and  2.5.  UPS  results  for  400  point  grids  with  the  Baldwin-Lomax  turbulence  model  are  plotted  in  a 
similar  generalized  form,  but  are  plotted  against  Re*  rather  than  Ree. 

Despite  the  different  abscissa  units,  similar  trends  are  observed  in  the  experimental  and  computed 
data.  For  example,  both  sets  show  good  agreement  with  Sommer-Short  at  moderate  Reynolds 
numbers  and  good  agreement  with  van  Driest  II  at  high  and  low  Reynolds  numbers.  That  is,  both  the 
computed  and  measured  data  indicate  somewhat  higher  curvature  than  the  Karman-Schoenherr 
formula  produces.  The  experimental  data  appears  to  exhibit  slightly  more  curvature  than  the  CFD 
results,  though  it  should  be  noted  that  the  data  at  the  lowest  Reynolds  number  was  taken  at  Mach  5.8. 
Recognizing  that  the  Spalart-Allmaras  model  produces  3-4%  higher  skin  friction,  the  scatter  in  the 
data  (10%  spread)  makes  it  difficult  to  choose  a turbulence  model  based  on  accuracy  on  a flat  plate. 
Other  problems  experienced  with  Spalart-Allmaras,  to  be  discussed  later,  led  to  the  decision  to  focus 
on  the  Baldwin-Lomax  results. 
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CHOOSING  A FLAT  PLATE  FORMULA 


A 


Generalized  Skin  Friction  (Mach  2.4) 


Prandtl  - 1/7th  Power 

- - Prandtt-Schlichdng 

Sivells-Payne 

White 


o — o UPS:  Sommer-Short  (Hi-Re) 
d — oUPS:  Spalding-Chi  (Hi-Re) 
x — " UPS:  Van  Driest  II  (Hi-Re)' 

• UPS:  Sommer-Short  (Lo-Re) 
■ — • UPS:  Spalding-Chi  (Lo-Re) 

* — « UPS:  Van  Driest  II  (Lo-Re) 


•COMPATIBLE  PAIRS: 


Lawrence- HSFVCA  - 19ol55 


(UPSBL:  SP/SC),  (UPSSA:  White/SC), 

(OFBL:  SP/SS),  and  (OFSA:  SPA/D) 

- Applied  Computational;  Aarodynamtcs  Branch 
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Given  a grid-converged  solution  at  both  wind  tunnel  and  flight  Reynolds  numbers,  the  next  step  is  to 
obtain  a flat  plate  method  that  gives  satisfactory  agreement  for  a wide  range  of  Reynolds  numbers. 
These  figures  show  computed  skin  friction  coefficients  in  generalized  form,  as  on  the  previous  page, 
but  plotted  in  terms  of  a percent  variation  from  Karman-Schoenherr. 

The  curvature  observed  on  the  previous  page  is  accentuated  here,  and  for  corrected  Reynolds  number 
greater  than  106,  the  computed  results  using  Baldwin-Lomax  all  (including  OVERFLOW  results,  shown 
in  red)  exhibit  a common  slope.  Similar  behavior  is  exhibited  by  the  Spalart-Allmaras  model.  The  slope 
in  the  B-L  results  at  moderate  to  high  Reynolds  number  agrees  fairly  well  with  that  of  the  Sivells-Payne 
formula.  The  different  compressibility  corrections  appear  mainly  to  shift  the  curves  by  different 
amounts. 

Finally,  it  can  be  seen  that  the  UPS-BL  results,  corrected  using  the  Spalding-Chi  compressibility 
formula,  agree  with  the  Sivells-Payne  incompressible  formula  to  within  2%  for  Re  > 105,  and  to  within 
1 % for  5x1 06  < Re  < 1 08.  Thus,  the  SP/SC  flat  plate  approach  was  used  for  the  remainder  of  the  study. 
The  best  flat  plate  method  for  the  S-A  model  appears  to  be  Spalding-Chi  applied  to  the  White 
incompressible  formula.  For  the  OVERFLOW  results:  using  B-L,  Sivells-Payne  with  Sommer-Short 
compressibility,  and  using  S-A,  Sivells-Payne  with  van  Driest  II  compressibility. 
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INTEGRATED  COMPONENT  FRICTION  DRAGS 


Including  UPSBL-Compatible  Method  (Sivells-Payne/Spalding  Chi) 


This  figure  shows  the  bar  chart  of  Figure  2,  with  results  of  the  Sivells-Payne/Spalding-Chi  average  skin 
friction  formula  included.  It  is  observed  to  fall  approximately  halfway  between  the  van  Driest  II  and  the 
Sommer-Short  results  for  both  the  wing  and  fuselage. 

The  UPS  results  were  computed  on  a four  zone,  axially  patched  mesh  with  a constant  wall-normal 
spacing  of  0.01  inches,  airplane  scale.  This  spacing  yields  y+  values  at  the  wall  of  less  than  0.5  over 
the  fuselage  and  wing  acreage  and  less  than  3 over  virtually  all  the  surface.  The  average  number  of 
wall-normal  points  in  the  boundary  layer  is  between  45  and  70  for  the  four  zones.  The  total  number  of 
wall-normal  points  used  was  120. 

The  Baldwin-Lomax  results  are  less  than  a half  count  lower  than  SP/SC  on  the  fuselage,  but  are  still 
3.5  counts  below  flat  plate  on  the  wing. 
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The  method  in  which  the  flat  plate  formula  is  applied  can  significantly  affect  the  resulting  flat  plate 
friction  drag  estimate.  Two  alternative  methods  are  shown  schematically  in  this  figure.  The  line- 
integration  approach  shown  on  the  right  has  been  used  to  generate  the  flat  plate  estimates  shown  in 
the  previous  bar  charts.  As  described  earlier,  this  method  uses  a line  integration  of  the  average  skin 
friction  coefficients  along  the  trailing  edge  for  the  wing  estimate,  and  a single  evaluation  of  the  average 
skin  friction  at  the  aft  end  of  the  body  for  the  fuselage  estimate.  The  areas  of  the  spanwise  wing  strips 
are  computed  using  the  strip's  average  arc  length.  A similar  approach  is  taken  on  the  fuselage,  before 
subtracting  the  airfoil  area  of  the  innermost  wing  section. 

With  the  surface  integration  method,  the  flat  plate  formula  for  local  skin  friction  is  used,  along  with  a 
function  that  defines  the  wing  leading  edge  (X|e  as  a function  of  span  location),  to  perform  an 
integration  of  local  skin  friction  coefficient  over  the  surface  of  the  configuration.  This  is  essentially  the 
same  method  as  that  used  to  evaluate  friction  drag  from  a CFD  solution. 

There  is  a difficulty  with  the  flat  plate  formula  integration  that  the  CFD  integration  doesn’t  have, 
however.  Given  a value  of  local  skin  friction  coefficient  on  a surface  element,  it  is  not  clear  in  which 
direction  the  skin  friction  is  acting,  i.e.,  what  fraction  of  the  surface  element  area  should  be  included. 
The  two  simplest  approaches  are:  total  area  (dStd  = dS),  and  two-dimensional  ( dStd  = dSplan , 

wing  and  dStd  = dSr  , fuselage).  These  methods  produce  significantly  different  results  near  the 

wing  leading  edge,  and  the  effect  on  integrated  drag  is  significant  as  well.  The  first  method  is  the  most 
like  the  line-integration  approach  on  the  wing.  The  second  method  essentially  assumes  a chordwise 
flow  direction.  Although  the  assumption  of  chordwise  flow  is  not  accurate  at  the  leading  edge,  it  does 
seem  to  be  the  logical  choice  when  using  x - xle  as  the  characteristic  length.  There  would  appear  to 

be  some  cancellation  of  errors  with  this  approach,  as  well.  That  is,  the  underestimation  of  the 
component  of  the  surface  velocity  in  the  drag  direction  (especially  along  the  leading  edge)  is 
compensated  by  an  overestimation  of  the  skin  friction  coefficient  caused  by  assuming  new  boundary 
layer  is  being  formed  at  the  leading  edge. 
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FLAT  PLATE  APPROXIMATIONS 

FLAT  PLATE  WALL  SHEAR:  zFJall  = P„vlcf(M„,  Re„x) 

APPROXIMATIONS: 

-(W)  EDGE  CONDITIONS  = FREESTREAM  CONDITIONS 
-(W)  FLOW  IS  2D,  i.e.,  s(x,z)  = x - x,e 

-(W)  NO  STREAMWISE  PRESSURE  GRADIENT 
-(W/B)  NO  INTERFERENCE  EFFECTS 
-(B)  NO  BODY  GROWTH  IN  THE  AXIAL  DIRECTION 
-(B)  NO  SPANWISE  CURVATURE 
-(W/B)  ETC. 

SOME  POSSIBLE  MODIFICATIONS  TO  FLAT 
-MFPI:  t - p ■y,C,(M„Re„lx), 

- “FP":  - py,CtW,.Re,Ax). 

- MFPIII:  Cf'"  = p 

of  55  ' — App*e<?  Computation*!;  Aftxtymmc*  Br*nch 


PLATE  WALL  SHEAR: 
ax  = x - xte 

AX  = X - Xte 

s = surface  streamline  length 


In  addition  to  the  ability  to  account  for  some  geometric  effects,  the  surface  integration  method  allows 
tor  the  inclusion  of  local  flow  variations  in  the  estimate  of  local  wall  shear. 

The  line-integration  method  forces  the  use  of  freestream  quantities  to  evaluate  compressibility  effects 

Skiri  fnc?!?n  coefficient  int0  wal1  shear  (the  author  attempted  to  include  some  angle- 
of-attack  effects  by  using  oblique  shock  and  Prandtl-Meyer  theory,  but  the  effects  were  found  to  be  less 
than  25  /«,  of  the  angle-of-attack  sensitivity  observed  in  CFD  solutions).  This  simple  formula  comes  at 
the  expense  of  neglecting:  1 ) the  effects  of  local  inviscid  flow  variations  (edge  conditions)  2)  three- 
dimensionality,  3)  effects  of  streamwise  pressure  gradient,  4)  interference  effects,  5)  effects  of 
increasing  body  circumference  in  the  streamwise  direction,  6)  spanwise  curvature  effects,  and 
probably  many  others.  The  various  effects  are  labeled  above  according  to  whether  they  are  believed  to 
T°[f  ®J9nrficant  on  the  wing,  on  the  body,  or  both.  Although,  spanwise  curvature  is  associated  with 
the  body  here,  there  is  a significant  spanwise  curvature  effect  on  the  TCA  wing  caused  by  the  landina 
gear  box,  as  will  be  seen.  y y 


The  first  of  these  effects  is  believed  to  be  the  most  significant  Fortunately,  it  is  straightforward  to 
include  edge  conditions  into  the  flat  plate  formula  for  local  skin  friction.  They  can  be  included  at  two 
levels:  first,  the  local  edge  dynamic  pressure  can  be  used  to  convert  the  skin  friction  coefficient  into 

W?T.  ®ar>  fGcond>  local  edge  Mach  number,  Reynolds  number,  and  temperature  can  be  used 
within  the  calculation  of  the  skin  friction  coefficient  itself. 


A third  effect  was  also  investigated.  The  effect  of  three-dimensional  surface  streamlines  was  studied 
by  computing  streamlines  backward  in  time  from  the  given  location  on  the  wing  to  where  they 
encounter  the  leading  edge.  The  lengths  of  these  streamlines  were  computed  and  compared  with  the 
local  x — X[e,  and  the  skin  friction  coefficient  adjusted  accordingly.  It  is  not  practical  to  perform  this 

operation  for  each  surface  element,  so  spot  checks  at  three  axial  locations  were  performed  to  obtain  a 
general  feeling  for  the  level  of  importance  of  this  effect. 
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The  effects  of  some  of  the  modifications  discussed  on  the  previous  two  pages  are  quantified  in  this  bar 
chart  of  the  wing  friction  drag.  The  notation  is  as  follows: 

-FP/LI  is  Sivells-Payne/Spalding-Chi  applied  using  line-integration 

-FP/SI  is  Sivells-Payne/Spalding-Chi  applied  using  2D  surface-integration 

-MFPI  is  FP/SI  scaled  locally  by  qedge/q„ 

-MFPII  is  SP/SC(M edge,  Reedge,  Tedge),  scaled  locally  by  qedgt/q„ 

-UPSBL  is  UPS  with  Baldwin-Lomax  turbulence 

Significant  geometric  effects  are  observed  as  exhibited  by  the  one  count  difference  between  FP/LI  and 
FP/SI.  This  drop  is  almost  entirely  associated  with  the  wing  leading  edge  curvature. 

Even  more  significant  is  the  effect  of  boundary-layer-edge  dynamic  pressure,  accounting  for  another 
1.9  count  reduction.  The  average  dynamic  pressure  on  the  surface  tends  to  be  5-10%  below  the 
freestream  value  and  has  a surface  distribution  very  similar  to  the  pressure,  though  scaled  by  a 
function  of  the  ratio  of  edge  to  freestream  Mach  number.  Thus,  the  upper  surface  tends  to  show  lower 
wall  shear  than  the  lower  surface  of  the  wing. 

The  effect  of  including  the  edge  conditions  in  the  evaluation  of  the  local  skin  friction  coefficient  is  small 
compared  with  the  dynamic  pressure  effect.  Also,  it  is  not  clear  to  the  author  whether  the  use  of  local 
edge  conditions  (virtually  extending  a flat  plate  at  the  local  flow  conditions  to  the  leading  edge) 
provides  a more  physically  realistic  result  than  simply  using  freestream  quantities. 

On  the  right,  the  bar  chart  is  scaled  up,  and  figures  are  included  indicating  the  percentage  difference 
between  the  associated  flat  plate  estimate  and  the  CFD  result.  Note,  nearly  80%  of  the  difference 
between  the  original  flat  plate  application  and  the  CFD  can  be  explained  by  geometric  and  inviscid  flow 
effects.  Also  note,  these  numbers  do  not  correspond  exactly  to  the  results  present  earlier  because  a 
slightly  different  line  of  demarcation  was  used  to  separate  the  wing  and  fuselage  on  the  lower  surface. 
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To  provide  a qualitative  idea  of  the  potential  effects  of  three-dimensional  flow  on  the  wing  surface, 
surface  streamlines  are  plotted  in  this  figure.  The  surface  is  colored  by  the  normalized  difference 
between  the  UPS  results  and  the  flat  plate  result  generated  by  MFPII  (full  use  of  edge  conditions)  The 
color  scale  extends  from  -20%  (blue)  to  +20%  (white).  A white  contour  line  is  shown  for  the  value  of 
zero  (UPS  and  MFPII  give  the  same  wall  shear).  Values  are  normalized  by  FP/SI.  Three  locations  are 
indicated  where  the  skin  friction  coefficients  will  be  examined  in  more  detail,  including  an  attempt  to 
quantify  three-dimensional  streamline  effects.  Finally,  a line  plot  is  included  that  shows  the  axially 
accumulated  difference  in  skin  friction  drag  on  the  upper  and  lower  wing  surfaces. 

Streamlines  on  the  lower  surface  indicate  a fairly  two-dimensional  flow,  with  a slight  outboard  flow  in 
the  forward  part  of  the  wing.  Note  that,  due  to  the  leading  edge  sweep,  an  outboard  flow  tends  to 
produce  a streamline  length  somewhat  larger  than  * - xle  while  an  inboard  flow  tends  to  produce  a 

streamline  length  less  than  x — x . Inboard  flow  is  observed  on  the  forward  part  of  the  upper  surface. 

Based  on  the  streamline  shown  here,  however,  the  effect  of  three-dimensionality  would  not  appear  to 
be  as  significant  as  the  inviscid  effects  previously  discussed. 

The  colored  surface  indicates  regions  on  the  upper  surface  associated  with  the  recompression  shock 
where  th©  CFD  results  fall  below  flat  plate  and  closer  to  the  leading  edge  on  the  lower  surface. 
Regions  where  the  viscous  computations  give  higher  shear  than  flat  plate  are  observed  on  the  upper 
surface  toward  the  inboard  trailing  edge,  on  the  biconvex  section,  and  near  the  leading  edge  just 
inboard  of  the  break.  On  the  lower  surface,  there  are  regions  of  relatively  high  shear  toward  the  trailing 
edge  and  at  the  edge  of  the  landing  gear  box.  Some  oscillatory  behavior  is  seen  in  the  UPS  results 
toward  the  trailing  edge.  This  behavior  is  the  result  of  variation  in  the  wall-normal  spacing  generated 
be  the  HYPGEN  hyperbolic  grid  generator  in  response  to  streamwise  clustering  of  the  surface  grid 
around  the  wing-tip  corner. 

The  line  plot  indicates  that  0.4  counts  of  the  friction  drag  difference  lies  on  the  upper  surface,  while  0.3 
counts  lies  on  the  lower  surface.  The  slight  rise  in  both  curves  near  the  front  of  the  wing  is  associated 
with  the  fact  that  the  flat  plate  result  goes  to  zero  at  the  leading  edge  due  to  the  2D  assumption, 
whereas  the  CFD  results  indicate  a significant  spanwise  flow  at  the  leading  edge  (especially  inboard). 
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The  first  location  to  be  studied  is  at  x=1900”,  which  is  in  the  mid  portion  of  the  subsonic  wing 
section. UPS  results  for  local  skin  friction  coefficient  are  plotted  in  comparison  with  four  editions  of  the 
flat  plate  local  skin  friction  coefficients.  The  upper  left  figure  shows  the  flat  plate  method  with  no  local 
flow  corrections,  the  upper  right  figure  shows  the  MFPI  result,  the  lower  right  shows  the  MFPII  result, 
and  the  lower  left  shows  the  MFPII  corrected  to  account  for  local  streamline  length.  Each  figure  has  an 
associated  figure  in  red  indicating  the  integrated  difference  at  this  x location  between  the  CFD  result 
and  the  flat  plate  result. 

The  pure  flat  plate  method  is  observed  to  provide  good  agreement  with  the  CFD  on  the  lower  surface, 
at  this  location,  but  significantly  overpredicts  the  CFD  on  the  upper  surface,  where  the  dynamic 
pressure  is  significantly  lower  than  freestream.The  bump  in  the  UPS  results  for  the  inboard  lower 
surface  is  associated  with  the  landing  gear  box. 

Scaling  the  flat  plate  skin  friction  coefficient  by  the  normalized  edge  dynamic  pressure,  qedge/ 

provides  surprisingly  good  agreement  with  the  CFD  results  on  both  the  upper  and  lower  surfaces.  The 
flat  plate  results  fall  consistently  lower  than  the  CFD  results  and  there  is  a slight  dip  observed  in  the 
CFD  results  that  is  apparently  not  explained  by  dynamic  pressure  effects. 

As  was  observed  in  the  integrated  friction  drag  values,  very  little  effect  is  observed  when  local 
conditions  are  used  to  evaluate  the  local  skin  frictions.  A slight  tendency  is  observed  toward  lower 
values  than  MFPI  on  the  upper  surface  and  higher  values  on  the  lower  surface. 

Finally,  adjusting  for  streamline  length  appears  to  have  little  effect  on  the  average  skin  friction  at  this 
location.  Some  improvement  in  the  agreement  with  CFD  is  observed  on  the  lower  surface,  but 
agreement  on  the  upper  surface  is  worse. 
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Similar  results  are  shown  downstream  near  the  leading  edge  break. 

aom« ' ho ttor^n  Sfi  Plat®  te^t?v!!ard  the  lower  surface  CFD  skin  friction.  The  MFPI  results 

agree  better  on  the  lower  surface  at  this  location,  except  toward  the  leading  edge.  On  the  uoDer 

surface,  an  inflection  is  observed  in  the  CFD  skin  friction  which  is  present  to  a much  lesser  dearee  in 
the  flat  plate  result.  Inboard  the  agreement  between  the  theory  and  the  computation  appears  to  be 
!™Pr°vingh  Even  le®f  influence  of  edge  Mach,  etc.,  is  observed  at  this  location  than  at  x=1900”  At  this 
location,  however,  the  integrated  effect  of  using  streamline  length  is  significant.  The  influence  is 
primarily  locatednear  the  lower  surface  leading  edge,  where  the  outboard  flow  was  observed  in  the 

SS'ISJ?0'-  1^ou?n"r’  a?aln’  ’an,ia  10  incraasa »» 

the  Reynolds  number  at  which  the  skin  friction  function  is  evaluated.  " " 
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Finally,  results  are  presented  along  a line  which  passes  near  the  front  of  the  wing  tip.  This  was, 
perhaps,  an  unfortunate  location  to  study,  since  it  includes  results  from  within  the  oscillatory  region  of 
the  CFD  solution.  The  computed  skin  friction  data  at  this  location  appears  to  be  slightly  lower  than  that 
slightly  upstream  and  downstream,  generally,  as  can  be  seen  in  the  following  figure. 

Somewhat  similar  trends  are  observed  as  in  the  previous  two  figures.  Again,  the  relatively  low  CFD 
data  indicated  outboard  of  the  break  is  believed  to  be  anomalous  to  this  axial  station.  Inboard  results 
show  good  agreement  between  the  UPS  results  and  all  three  modified  flat  plate  theories. 
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In  this  figure,  the  surface  is  again  colored  as  in  the  streamline  figure,  but  contours  of  surface  Cp  are 
overplotted. 

It  appears  possible  to  qualitatively  correlate  regions  of  relatively  low  CFD-predicted  skin  friction  with 
adverse  pressure  gradient  regions,  and  regions  of  relatively  high  CFD-predicted  skin  friction  with 
favorable  pressure  gradient  regions.  The  primary  exception  to  this  rule  is  the  “hot  spot"  along  the 
landing  gear  box  which  is  generated  by  the  lateral  convex  curvature.  The  pressure  gradient  also 
influences  the  surface  streamlines  in  such  a way  as  to  compound  its  effect,  i.e.,  favorable  pressure 
gradient  near  the  leading  edge  pulls  the  streamlines  inboard,  and  adverse  pressure  gradient  near  the 
leading  edge  pushes  the  streamlines  outboard. 

The  effects  of  pressure  gradient  on  the  local  skin  friction  coefficient  have  not  yet  been  quantified. 
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Similar  analysis  was  performed  for  the  TCA  fuselage  flow  field  and  the  “bottom-line”  results  are  shown 
in  this  figure. 

The  trend  for  the  flat  plate  modifications  is  similar  to  that  observed  for  the  wing.  One  may  wonder, 
however,  about  the  relatively  large  difference  between  the  line-integrated  and  surface-integrated 
results,  since  there  is  no  leading  edge  effect  in  the  surface-integration  of  the  fuselage.  The  reason  for 
the  difference  is  the  effect  of  the  shrinking  nose  circumference.  The  surface  integration  effectively 
produces  a weighted  average  skin  friction  coefficient,  with  the  high-local-Cf  nose  region  weighted  less 
because  of  the  lower  circumference.  Similarly  on  the  wing,  the  surface  integration  weights  the  high- 
local-Cf  leading  edge  less  because  the  flow  direction  is  not  in  the  drag  direction. 

As  on  the  wing,  the  average  surface  dynamic  pressure  is  lower  than  freestream  and  the  edge-modified 
theories  show  a lower  integrated  friction  drag  than  the  freestream-based  methods. 

In  the  case  of  the  fuselage,  however,  the  CFD  result  is  significantly  higher  than  those  of  the  modified 
flat  plate  methods,  falling  close  to  the  line-integrated  result.  The  following  figures  will  attempt  to,  at 
least  qualitatively,  explain  this  trend. 
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TCA  FUSELAGE  SURFACE  PRESSURE  w/  ACf 

Mach  2.4,  Re  = 4 Million/ft,  « = 3.6 


ACCUMULATED 


Lawrence-HSK/CA  -41  or  55 


This  plot  shows  the  fuselage  surface  colored  using  the  same  normalized  skin  friction  difference  as 
shown  in  previous  figures  focusing  on  the  wing  flow  field.  Here,  the  axial  distance  is  compressed  to 
highlight  the  fuselage  flow  and  the  wing  is  cut  out.  The  color  scale  here  extends  from  -50%  to  +50%, 
rather  than  the  narrower  scale  used  in  the  previous  figures.  Surface  pressure  contours  are  overplotted 
Also  included  is  a plot  of  axially  accumulated  friction  drag  difference. 

The  line  plot  indicates  that  more  than  90%  of  the  difference  between  the  modified  flat  plate  and  the 
UPS  friction  drag  estimates  is  contained  on  the  forebody,  upstream  of  the  wing.  At  the  point  of 
maximum  fuselage  radius,  the  curve  sharply  flattens. 

Other  “hot  spots”  on  the  fuselage  are  associated  with  the  flow  of  wing  boundary  layer  onto  the 
fuselage  behind  the  wing  shock.  The  wing  shock  provides  an  adverse  pressure  gradient  lowering  skin 
friction  and  diverting  flow  away  from  the  wing,  and  the  aft  wing  expansion  sweeps  across  the  fuselage 
creating  a relatively  high  C(  region.  The  low  speed  flow  from  the  wing  boundary  layer  is  observed, 
superimposed  on  the  fuselage.  Finally,  some  colorful,  and  perhaps  non-physical,  effects  are  observed 
on  the  aft  fuselage.  Although  more  pronounced  in  this  figure  than  the  forebody  coloring,  it  should  be 
remembered  that  the  normalization  used  here  tends  to  accentuate  differences  at  higher  Reynolds 
number.  Thus,  the  relatively  small  effect  on  accumulated  friction  drag  observed  in  the  line  plot. 
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Flying  Fuselage  Test  Case 

PRIMARY  REASON  FOR  4%  DIFFERENCE  BETWEEN  UPS-BLAND  SI**,: 
-TURBULENT  CONE  FLOW  EFFECT  FOR  -40%  OF  BODY  LENGTH 


x,  inches 

^wrenceTtSR/CX  'rc  5IT5 


White:  For  turbulent  flow, 

0.087  < c-one-  - 1 <0.176 

^ f,  plate 


CFD  shows  elevated  Cf  where  ^>o 

? /ir 
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In  an  effort  to  understand  the  forebody  flow  field,  an  axisymmetric  flow  case  was  computed  using  the 
axial  radius  distribution  corresponding  to  the  TCA  fuselage.  This  is  referred  to  as  the  flying  fuselage 
case.  The  radius  distribution  is  plotted  at  the  bottom.  Skin  friction  computed  by  UPS  is  compared  with 
flat  plate  (MFPII)  in  the  middle  plot,  and  the  percentage  difference  in  axially  accumulated  friction  drag 
is  plotted  at  the  top. 

Similar  behavior  is  observed  in  this  axisymmetric  case  as  was  seen  in  the  3D  case  - significantly 
higher  Cf  on  the  portion  of  the  body  where  the  radius  in  increasing  in  the  axial  direction.  This  level  of 
excessive  skin  friction  is  consistent  with  that  experienced  in  cone  flow,  where  the  growth  of  the  body 
produces  a boundary  layer  that  is  thin  relative  to  a flat  plate  or  cylinder.  White  (Viscous  Fluid  Flow,  2nd 
Ed)  provides  estimated  bounds  on  the  level  of  the  effect  of  between  8.7%  and  17.6%  for  turbulent  flow. 
At  maximum  radius,  the  integrated  CFD  friction  drag  is  observed  to  be  8.65%  higher  than  the  flat  plate 
theory.  This  produces  an  overall  increase  for  the  body  of  3.94%,  which  is  in  fairly  good  agreement  with 
the  level  of  4. 3-4.7%  observed  in  the  wing/body  case,  considering  the  reduced  fuselage  area  in  the 
wing/body  case. 
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In  the  category  of  odd  and  ends,  this  figure  shows  the  difference  in  skin  friction  produced  by  the  UPS 
code  using  two  different  turbulence  models:  Baldwin-Lomax  and  Spalart-Allmaras.  A black  contour  line 
at  zero  is  included  to  separate  regions  of  higher  B-L  skin  friction  from  those  of  higher  S-A.  Also  shown 
is  a plot  of  the  axially  accumulated  differences  on  the  wing  and  fuselage. 

As  indicated  by  the  flat  plate  study,  the  S-A  model  produces  higher  skin  friction  over  most  of  the  wing 
and  fuselage.  However,  there  are  a few  regions  where  the  B-L  model  produces  higher  Cf,  most  notably 
near  the  wing  leading  edge.  The  S-A  model  was  observed  to  produce  a significant  undershoot  in  Cf 
near  the  leading  edge,  especially  on  the  outboard  section.  Some  hint  of  this  behavior  was  experienced 
in  the  flat  plate  case  as  well,  where  attempts  to  trip  the  boundary  layer  at  arbitrarily  low  Reynolds 
numbers  failed.  This  effect  is  also  evident  in  OVERFLOW  solutions  of  the  TCA  wing/body  and  wing/ 
body/nacelle/diverter  (low  Cf  rings  are  seen  near  the  nacelle  leading  edges)  to  a somewhat  greater 
degree  than  in  the  UPS  solutions,  perhaps  because  of  a coarser  axial  grid  spacing  than  that  used  in 
the  UPS  solutions.  The  integrated  effect  is  observed  in  the  wing/fuselage  breakdowns  shown  in  the 
first  two  bar  charts.  Interestingly,  at  this  Reynolds  number,  the  leading  edge  effect  almost  exactly 
cancels  the  tendency  toward  higher  flat  plate  skin  friction.  At  flight  conditions,  the  leading  edge 

undershoot  was  not  observed  and  S-A  produced  higher  skin  friction  drag  on  both  the  wina  and 
fuselage. 


1474 


Finally,  the  analysis  described  here  for  the  TCA  was  also  applied  to  other  cases.  The  integrated  results 
for  the  TCA  wing  at  flight  conditions,  the  1-03  wing,  the  Ref.  H wing,  and  the  7-04  wing  are  shown  in 
this  bar  chart. 

The  trends  for  the  TCA  wing  at  flight  conditions  and  the  1 -03  wing  are  very  similar  to  the  trends  for  the 
TCA  wing  at  tunnel  conditions.  At  flight  conditions,  the  MFPII  friction  drag  is  2.3%  higher  than  the  UPS 
result  compared  with  1 .9%  at  tunnel  conditions.  Thus,  the  method  appears  relatively  insensitive  to 
Reynolds  number. 

The  results  for  the  Ref.  H and  Ref.  H-based  bodies  are  somewhat  different  in  that  the  UPS  predicts 
higher  friction  drag  than  the  modified  flat  plate  theories.  A surface  plot  of  skin  friction  differences 
indicates  the  presence  of  chordwise  flow  features  produced  by  the  UPS  simulations  near  the  wing 
fuselage  junction.  These  features  are  reduced  in  intensity  on  the  7-04  configuration,  but  are  not 
removed.  At  present,  it  is  not  known  whether  these  are  physical  flow  features  associated  with  the 
higher  leading  edge  sweep  on  the  inboard  Ref.  H wing  section,  or  whether  they  are  numerical  oddities 
produced  by  the  Baldwin-Lomax  model. 
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SUMMARY 


FLAT  PLATE  STUDY 

- S-A  MODEL  GIVES  3-4%  HIGHER  SKIN  FRICTION  THAN  B-L 

- UPS  GIVES  HIGHER  SKIN  FRICTION  THAN  OVERFLOW  (1%) 

- QUADRATIC  IMPROVEMENT  WITH  AVERAGE  BOUNDARY-LAYER  Ah 

UPSbl-COMPATIBLE  FLAT  PLATE  METHOD 

-SIVELLS-PAYNE  INCOMPRESSIBLE  FORMULA 
-SPALDING-CHI  COMPRESSIBILITY  CORRECTIONS 

DIFFERENCES  BETWEEN  UPS  AND  FLAT  PLATE  QUANTIFIED  (WING): 

-TOTAL  DIFFERENCE  = 3.5  ct 

• LEADING  EDGE  GEOMETRIC  EFFECT  = 1 ct 

• BOUNDARY-LAYER  EDGE  DYNAMIC  PRESSURE  = 1 .8  ct 

QUALITATIVE  EXPLANATIONS  FOR  REMAINING  VARIATIONS  (WING  AND  BODY) 
-PRESSURE  GRADIENT  EFFECTS  AND  CONVEX  CURVATURE  (WING) 
-CONE  FLOW  AND  INTERFERENCE  EFFECTS  (BODY) 

TURBULENCE  MODEL  SENSITIVITIES  OBSERVED 

L*wr«nc*-HJ>«A^A  - of  S5  * Computation*;  Afoctymmcs  Brunch  -S 


A summary  of  the  present  study  is  provided  in  this  chart. 


A flat  plate  study  was  performed  to  understand  the  grid  sensitivities  of  the  Baldwin-Lomax  and  Spalart- 
Allmaras  turbulence  models,  and  to  obtain  a true,  grid-independent  UPS  solution  to  the  flat  plate 
problem  at  cruise  Mach  number.  F 


The  Sivells-Payne  incompressible  skin  friction  formula  in  conjunction  with  the  Spalding-Chi 
compressibility  correction  was  found  to  give  good  agreement  with  the  UPS  B-L  results  over  a wide 
range  of  Reynolds  numbers. 


observed  to  account  for  approximately  80%  of  the  difference 
between  the  flat  plate  and  CFD  friction  drag  estimates  on  the  TCA  wing. 


The  physical  location  of  pressure  gradients  on  the  wing  were  observed  to  be  correlated  with  the 
location  of  the  most  noticeable  local  differences  between  the  modified  flat  plate  theory  and  the  CFD 
results  On  the  body,  effects  of  the  increasing  body  radius  over  the  upstream  30%  of  the  body  could 
account  for  the  elevated  CFD  prediction  for  fuselage  friction  drag. 


The  Spalart-Allmaras  turbulence  model  tends  to  produce  higher  skin  friction  when 
state,  but  tends  to  revert  to  laminar  flow  at  low  local  Reynolds  number. 


in  a fully  turbulent 
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CONCLUSIONS 


FLAT  PLATE  APPROXIMATIONS 
-NO  “BEST”  METHOD 

• SOMMER  AND  SHORT  WORKS  WELL  FOR  TUNNEL  CONDITIONS 

• VAN  DRIEST  II  WORKS  WELL  AT  FLIGHT  CONDITIONS 

• SIVELLS-PAYNE/SPALDING-CHI  COMPARES  WELL  WITH  CFD 
FOR  A WIDE  RANGE  OF  RE 

FLAT  PLATE  APPLICATION 

-80%  IMPROVEMENT  IN  FLAT  PLATE  PREDICTIONS  CAN  BE 
OBTAINED  USING  SURFACE  INTEGRATION  WITH  AN  INVISCID  SOLN 

N-S  SOLUTIONS  NEEDED  TO  DETERMINE  INFLUENCE  OF: 

-INTERFERENCE 
-PRESSURE  GRADIENT 
-LATERAL  CURVATURE 
-FUSELAGE  CONE-FLOW  EFFECT 


Oh»t»nc»HbMra-T;r 


otst 
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Some  general  conclusions  are  given  here. 

The  study  failed  to  determine  a universally  best  flat  plate  method.  Various  methods  work  better  for 
different  conditions.  For  comparison  with  CFD,  it  seems  appropriate  to  choose  a method  which  agrees 
with  a grid  independent  solution  of  the  chosen  N-S  solver  on  a flat  plate. 

Flat  plate  method  application  appears  to  be  as  important  as  the  particular  method  chosen.  In 
particular,  surface  integration  using  in  viscid  dynamic  pressures  produces  much  better  agreement  with 
CFD  than  line  integration.  This  would  allow  first-order  friction  drag  sensitivities  to  be  included  within 
the  design  process. 

Second-order  effects  which  still  require  Navier-Stokes  analysis  are  listed. 
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The  work  to  be  described  was  performed  at  the  NASA  Langley  UPWT  (4-ft  supersonic), 
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Method  (1) 
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• According  to  theory  of  creeping  flow,  the  slope  will  be  approximately 
constant  near  the  leading  edge  of  a line  of  oil  ( dy/dx  = (i  / (x  t)  ). 


A line  or  dot  of  transparent  silicone  oil  (DC-200  Fluid,  Dow-Corninq,  Mi 
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Method  (2) 
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First  of  all,  it  needs  to  be  established  that  the  proposed  method  works  at  all  under  the 
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Note  inclusion  of  the  machinist's  scale  to  provide  a length  and  axis  reference. 


Measured  and  Computed  Skin  Friction  Coefficients 
^ef  H (Upper  Surface,  Outboard)  at  Mach  = 2.40,  Re  = 3 million/ft 

B Turbulent:  Cf_x  from  UPS  Code  (Fuselage  Station  2282") 
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Measured  and  Computed  Skin  Friction  Coefficients 
Ref  H (Upper  Surface,  Outboard)  at  Mach  = 2.40,  Re  = 3 million/ft 
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Span  Station  (inches) 
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Measured  and  Computed  Skin  Friction  Coefficients 
Ref  H (Upper  Surface,  Midwing)  at  Mach  = 2.40,  Re  = 3 million/ft 


Span  Station  (inches) 
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Rough  estimates  of  the  impact  on  Cf  of  the  various  error  sources  are  collected  together  to 
guide  the  planning  of  a follow-on  "real"  experiment.  It  should  be  noted  that,  in  contrast  to 
forces  and  pressures,  a measurement  of  local  skin  friction  values  to  even  ±20-25%  is  not 
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allow  routine  Cf  measurements  on  wind  tunnel  models, 
new  capability. 
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The  computational  fluid  dynamics  (CFD)  comparisons  being  presented  are  compared  to  each  other  and  to  wind  tunnel 
(WT)  data  on  the  baseline  TCA.  Some  of  the  CFD  computations  were  done  prior  to  the  tests  and  others  later.  Only  force 
data  (CL  vs  CD)  from  CFD  will  be  presented  as  part  of  this  report. 
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Baseline  TCA  Model  2a 
Nacelles  Installed  Without  Fairings 
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Baseline  TCA  Model  2b 
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Baseline  TCA  Model  2b 
Nacelles  Installed  with  Fairings 
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This  figure  is  a close  up  of  the  nacelles  installed  on  model  2b.  Notice  the  build  up  of  the  fairings.  The  fairings  were  built 
up  using  templates  to  mold  bondo  onto  the  surface.  This  process  is  not  very  exact  and  takes  about  an  hour  to  perform. 
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WB  Only 

Viscous  Drag:  Baldwin-Lomax 


Each  organization  was  given  an  iges  file  of  the  baseline  TCA  OML  that  the  Technology  Integration  ITD  provided  to  the 
HSR  organization.  The  CFD  results  presented  were  computed  at  the  WT  Reynolds  number  and  forces  were  integrated 
up  to  the  location  the  wind  tunnel  model  ended.  A configuration  feature  each  organization  modeled  differently  was  the 
nacelle  installation.  The  only  two  organizations  that  modeled  the  nacelle  fairing  was  McDonnell  Douglas  in  CFL3D  and 
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Wind  Tunnel  vs  Full  Potential  and  Euler 
M=2.4,  RN=4  M/ft 
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Wind  Tunnel  vs  Full  Potential  and  Euler 
M=2.4,  RN=4  M/ft 


1518 


This  figure  shows  the  drag  polar  zoomed  in  around  cruise  of  CL=0.1  for  comparison  of  Full  Potential  and  Euler  codes. 
There  is  a total  difference  in  pressure  drag  of  2 counts  for  WB  and  WBND.  The  CFD  drag  polars  are  to  the  right  of  WT 
data.  The  CFD  data  shows  a total  drag  difference  of  3 counts  for  WB  among  codes  and  3 counts  for  WBND. 
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When  comparing  OVERFLOW  to  WT,  the  differences  are  smaller  than  the  Full  Potential  and  Euler  results.  For  WB,  the 
drag  from  OVERFLOW/Ames  is  less  than  1 count  from  WT.  For  WBND,  the  WT  data  is  between  both  OVERFLOW 
solutions. 
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Wind  Tunnel  vs  CFL3D 
M=2.4,  RN=4  M/ft 


NASA  ARC/Cappuccio 


This  figure  compares  CFL3D  with  WT  data.  CFL3D  gives  higher  drag  than  WT  data  for  both  Euler  and  N-S  by  4 counts. 
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This  figure  shows  a comparison  of  CFL3D,  OVERFLOW/Boeing,  and  WT  data. 
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Wind  Tunnel  vs  N-S 
M=2.4,  RN=4  M/ft 
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Wind  Tunnel  vs  N-S 
M=2.4,  RN=4  M/ft 
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Wind  Tunnel  vs  CFD 
WB  at  M=2.4,  RN=4  M/ft 


NASA  ARC/Cappuccio 


The  total  drag  for  WB  ranges  over  5.5  drag  counts  between  all  the  CFD  and  WT 
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Full  Potential,  Euler,  and  N-S  CFD  Comparisons 
WB  Pressure  Drag 
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NASA  ARC/Cappuccio 


The  pressure  drag  among  all  the  codes  are  more  consistent  than  the  total  drag.  The  pressure  data  for  WB  ranges  over 
0.5  drag  count  for  N-S  and  2 drag  counts  for  Full  Potential  and  Euler. 
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WB  Viscous  Drag  Comparisons 
M=2.4,  RN=4  M/ft 


NASA  ARC/Cappuccio 


The  figure  shows  the  variation  of  viscous  drag  computation  for  the  WB.  The  spread  of  the  data  is  over  5.25  drag  counts. 
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Wind  Tunnel  vs  CFD 
WBND  at  M=2.4,  RN=4  M/ft 


NASA  ARC/Cappuccio 


The  total  drag  for  WBND  ranges  over  6.5  cts  between  all  of  the  CFD  and  WT. 
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Full  Potential,  Euler,  and  N-S  CFD  Comparisons 
WBND  Pressure  Drag 
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NASA  ARC/Cappuccio 


The  total  pressure  drag  for  WBND  ranges  over  all  the  CFD  codes  is  2 drag  counts. 
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WBND  Viscous  Drag  Comparisons 
M=2.4,  RN=4  M/ft 
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NASA  ARC/Cappuccio 


The  major  differences  lie  in  the  viscous  drag,  which  is  6 counts. 
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T-1679*  T-l  671  * T-1671*  AIRPLANE  SYN87MB  TRANAIR  CFL3D  CFL3D  OVERFLOW  OVERFLOW  OVERFLOW 

no  fairing  Euler  N-S  ARC  BCA  BCA  (IDP) 

Method 


o 


“O 

c 


o> 

c 


0 .== 
£ Si « 

0 C 3 
o>  0 ’q_  eg 

fps 

|lp 

l-s  i 5 

> ffi  -o 

^ CD  CD  ^ 

C C -i 
C © ® 5 


CD  © *Z 

o t;  -x  c 
c 3 i2  2 


© _ 
® © 

~ © 
TO  )- 
_ CD 

<5  > 


m C 

£ O 

£ § 

°-o 
sz  o 

O CO 


© 

c •*- 
.2  £ 
is 

S <D 

E E 

•^d  2 ® 

£*ls 


C 

o 


-O  o ® Ji 

s«.| 

■o  O « § 

£ S i rl 


— © 

« C 3 

^0(3: 
© ° © © 
3 © £ £ 
o -c  “ ~ 

o ^ ^ CO 
0 c c c 
> ® CO  CO 

Z o \ Z © 

|£iE 


CO 


CD  - 
£ O <D 
x:  x: 

i-  c/>  +* 

o * c 


I 

CO 

*°  CO 
p « CO 

E 

o 


o 

T2  3 


O) 

© 


S.SZs'? 

3 ® © © « 
— O C j3  2 

§ 5 ffl  2 8 

;*©£*-•> 
O S o a) 

Td  © £ 

© o r o 
I - Eg 

CD  .O 


* a ■=  1 1 

1 “ f i = 

® "O  0>  (ft 

sfs££ 

■=§©  = § 


CO 


1 £ 

c 

D 


^ 05 


C © 
® © 
o © 

E S’ 

o o © 

0 0 5= 


c CO  ® B 

2 u!  © 

o • a > 
O T3  „ o 
a>  ® ^ 
c 


> -C  CD 

*is£q 

cd  _c 


o 

CD 


c 

CO 


3 0 0 

~ « “P 


■O 

c 


c 

CD 

0 


C 

0 

E 

© 

u. 

O 

c 


i sS 


0 

X) 

0 

0 


0 
5 
o 

-C 
0 0 

CO  <S 
c CO 

0 "O 
£ 0 

fl 


0 “O  = -rj 
« .2  ® © 

o ® O 


c 

© 

E 

£ © 

.£  3 

£ 8 

8<? 
© t_ 
c •*- 

co  2 
© CO 

■°  d 

Q 

LL  C 

O g 

© «i= 
£ © 


|»l 

2 § o 
o o a 

£ ° r- 

CO  00  1- 


o 

0 

0 

0 

E 

0 

O) 

c 

0 
k— . 

-4—* 

c 

0 

E 

© 


•o  o o ® 

© -p  *;  O 
r-  C C c 

S 0 0 •== 

bee® 

c 2 £ © 

0 .E  c 


0 0 


_ 0 
— c 

C — 


QL  0 0 0 

1 s 

O 


ii 

“i »? 

£ « co  ■$ 

ll±I 


a. 

E 

o 

o 


€Qo 


0 
c 

lea  8 

111  £ CM 


C 

O O 


© 


> 5 


P O o 


1547 


Recommendations 
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Recommendations  for  both  CFD  and  WT  data  are  as  follows. 

One  item  that  is  still  not  being  modeled  in  CFD  like  being  tested  in  the  WT  is  where  the  boundary  layer  is  being  tripped 
to  go  turbulent.  Most  CFD  codes  are  modeling  fully  turbulent  flow  which  is  definitely  not  happening  in  the  WT. 

CFD  should  attempt  to  model  the  wind  tunnel  nacelle  internal  duct. 
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.551 


HSR  Aerodynamic  Performance  Workshop 
NASA  Langley  Research  Center 
February  25-28, 1 997 


Overall  Test  Program  Objectives 
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Program  Definition  Approach 


1553 


The  group  has  been  in  the  program  definition  stage  for  the  past  10  months.  The  approach  used  was  to  get  a large 
group  of  technical  people  together  to  define  the  objectives  of  this  milestone  and  define  the  PIE  technology  items.  These 
people  came  form  the  different  disciplines  needed  to  work  on  this  task.  They  came  from  propulsion,  inlet,  nozzle, 
configuration  aerodynamics,  and  technology  integration.  The  group  identified  these  PIE  technology  items,  assessed 
them  as  a function  of  key  parameters,  ranked  and  prioritized  them.  The  group  chose  the  highest  priority  items  for  the 
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Propulsion  Induced  Effects 
Technology  Objectives 
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The  technology  items  identified  by  the  larger  group  were  group  in  terms  of  inlet,  nozzle,  and  nacelle  flow.  Under  inlet 
flow  came  items  associated  with  the  installation  of  the  inlet  on  the  aircraft.  They  include  inlet  spillage  drag,  spillage 
interference  drag,  impact  of  inlet/diverter  design,  ram  drag,  inlet  flow  field,  and  inlet  unstart.  Under  nozzle  flows  came 
items  associated  with  the  installation  of  the  nozzle  on  the  aircraft.  They  included  isolated  nozzle  boattail  drag,  installed 
nozzle  boattail  drag,  impact  of  nozzle  design,  plume  effects  on  the  afterbody  and  tail,  climb/ejector  interactions,  and 
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Technology  Objective  Selection  Criteria 
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Priority  (Perf+S&C)  * Risk 


The  criteria  used  to  select  which  items  had  higher  priority  than  others,  was  based  on  many  things.  The  criteria  were 
airplane:  performance  and  stability  & control  impacts;  flight  regime;  impact  on  certification  requirements,  inlet/nozzle 
selection,  and  take  off  gross  weight.  The  assessment  for  each  of  these  categories  were  set  by  1=low  and  5=high.  Zero 
meant  not  applicable.  The  impact  was  either  a yes  or  a no.  Once  this  was  done,  the  group  had  to  decide  what  was  the 
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Propulsion  Induced  Effects  Technology  Items 
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Bleed  and  bypass  interference  was  ranked  priority  one.  Inlet  spillage  was  ranked  second.,  followed  by  inlet  unstart, 
climb/ejector  interactions,  nozzle  boattail,  and  afterbody  and  tail  interactions. 
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il 


Summary  of  Test  Objectives 
Presented  to  CA 
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Interactions 


Of  these  items  the  group  decided  on  inlet  bleed  & bypass  interference,  inlet  spillage,  and  installed  nozzle  boattail  as  PIE 
items  to  work  on  under  this  test  program.  Transonically  the  issues  are  inlet  spillage,  installed  nozzle  boattail,  and 
plume/tail  interactions.  Supersonically  the  issues  are  inlet  bleed  & bypass  interference,  and  plume/tail  interactions. 

The  priority  for  each  item  is  CFD  validation  and  experimental  data  is  required  to  meet  this  objective.  In  the  case  of  bleed 


1562 


Outcome  of  Review  to  CA 
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A review  of  what  the  PIE  team  had  decided  on  for  the  test  program  and  the  preliminary  design  by  MicroCraft  was 
presented  to  the  CA  team.  The  outcome  of  the  review  was  to  concentrate  on  the  transonic  testing  of  inlet  spillage  and 
jet  effects  on  the  nozzle  boattail  and  tail  interactions.  CA  said  to  incorporate  the  supersonic  testing  in  the  model  design 
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Current  Summary  of  Test  Objectives 
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So  now  these  are  the  current  test  objectives.  MicroCraft  has  been  asked  to  remove  the  bleed  design  from  the  model. 
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Test  Requirements 
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1569 


The  balance  of  choice  at  this  time  is  the  NASA  Ames  6 inch  diameter  flow  through  balance.  Based  on  the  TCA 
performance  data  acquired  to  date,  the  range  of  data  to  be  measured  is  within  the  balance  rhombus.  If  data  is  only 
required  at  the  cruise  angle  of  attack,  it  would  be  preferable  to  redesign  and  fabricate  a new  balance  where  data  could 
be  measured  over  the  full  range  of  the  balance  to  get  the  best  data  possible  for  PIE  increments. 
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Preliminary  Model  Design 
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Due  to  the  instrumentation  and  airflow  needed  to  simulate  inlet  and  nozzle  internal  flow,  the  model  was  scaled  at  5.7% 
to  be  able  to  fit  in  the  Ames  11x11  transonic  wind  tunnel.  The  model  scale  is  such  that  the  data  that  comes  out  of  the 
test  would  be  incremental  and  not  absolute,  due  to  shock  reflections.  This  scale  forces  the  model  to  be  tested  in  either 
the  Lewis  10x10  or  AEDC  16S  supersonically.  An  alternative  to  the  Ames  1 1x1 1 is  AEDC  16T.  The  model  is  a full  span 
model,  all  metric,  with  4 flowing  inlets  and  powered  nozzles.  Vertical  and  horizontal  tails  would  be  modeled  on  an 
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Model  Components 
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Model  components  would  include  a wing  body  made  entirely  of  aluminum  except  for  the  balance  block.  The  fuselage 
can  hold  up  to  100  ESP  psi  modules. 
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Flow  Through  Nacelle 
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1575 


The  aerodynamic  flow  through  nacelles  would  model  the  actual  nozzle  flaps  and  side  walls.  They  would  have  quick 
disconnects  at  the  diverter  interface  to  the  wing  so  pressure  instrumentation  would  be  easy  to  change  with  nacelle 
configuration  changes.  The  internal  duct  would  have  a constant  area,  from  inlet  to  exhaust,  as  the  aerodynamic 
performance  TCA  models  that  have  been  tested  to  date. 
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Inlet  Nacelle 
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The  inlet  nacelle  would  have  the  same  quick  disconnects.  The  internal  duct  would  be  modified  so  that  a centerbody 
spike  could  be  installed  to  model  inlet  spillage.  Basically  3 spike  positions  would  be  built:  Mach  < 1.65,  Mach=2  and 
2.4.  A spacer  would  adjust  the  position  of  the  spike.  In  addition  to  the  inlet  centerbody,  bypass  doors  would  be 
incorporated  into  the  inlet  model.  These  doors  can  be  set  to  the  different  positions  necessary  depending  on  the  Mach 
number  tested  at. 
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Propulsion  Nacelle 


The  propulsion  nacelle  is  the  one  that  would  be  used  for  the  jet  effects  part  of  the  test.  This  nacelle  would  have 
interchangeable  nozzle  parts,  quick  disconnects,  and  an  inlet  fairing.  The  nozzle  team  is  very  interested  in  this  test 
program  to  get  installed  nozzle  boattail  drag,  and  to  study  side  wall  effects.  They  plan  to  participate  and  fund  the  nozzle 
testing  on  this  model. 
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Model  Air  Delivery  System 
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The  model  air  delivery  system  will  take  on  board  high  pressure  air  through  the  balance  center,  go  through  flow  meters 
and  be  adjusted  by  metering  valves.  This  schematic  shows  3 lines  being  fed  to  each  nacelle.  This  would  be  correct  if 
we  were  going  to  be  modeling  bleed  in  addition  to  the  bypass  doors  and  nozzle,  since  bleed  has  been  eliminated,  only 
one  line  needs  to  be  fed  to  each  nacelle.  One  line  per  nacelles  would  be  used  to  supply  air  through  either  the  bypass 
doors  or  nozzle. 
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Nacelle  Configurations 
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The  nacelle  configurations  that  need  to  be  tested  are  the  aerodynamic  reference  flow  through  nacelle.  As  described 
earlier,  this  nacelle  would  be  similar  to  what  gets  tested  on  the  aerodynamic  performance  TCA  models. 
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Inlet  Spillage  and  Bypass  Nacelles 
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Nacelle 


The  inlet  nacelles  would  Include  an  inlet  spillage  configuration  and  an  inlet  spillage  and  bypass  door  configuration.  In 
both  cases,  a choke  plate  would  be  installed  down  stream  to  control  the  mass  flow  through  the  nacelle.  Total  pressure 
rakes  and  static  pressures  would  be  used  for  measurements  to  compute  the  mass  flow  through  the  nacelle.  Each 
nacelle  would  have  to  be  calibrated  in  a facility  outside  of  the  wind  tunnel  to  measure  the  internal  drag  for  making 
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Nozzle  Boattail  Nacelles 
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The  propulsion  or  nozzle  boattail  nacelles  are  made  up  of  a reference  nacelle  to  be  able  to  increment  off  the 
aerodynamic  reference  flow  through  nacelle,  a transonic  and  supersonic  reference  nacelle.,  and  operating  nacelles 
transonically  and  supersonically.  In  all  cases,  choke  plates  are  need  to  get  the  correct  nozzle  pressure  ratios  desired. 
Total  and  static  pressure  would  also  be  measured.  These  nozzles  would  need  to  be  calibrated  in  a facility  prior  to  the 
wind  tunnel  test. 
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Number  of  Nacelle  Types 
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Total  Number  Nacelles  to  be  Fab’d,  Calibrated,  & Tested  = (10  + 8 -1)  x 4 = 68 


There  would  be  a total  of  68  nacelles  that  would  need  to  be  fabricated,  calibrated,  and  tested.  The  table  describes  ail  of 
the  nacelles  that  would  have  to  be  tested.  The  inlet  spillage  nacelle  would  have  to  be  tested  at  4 mass  flow  ratios 
transonically.  4 operating  nozzles  shapes  transonically  would  also  have  to  be  tested.  In  total  10  types  of  nacelles 
would  be  tested  transonically. 
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Issues 
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The  issues  for  this  kind  of  a test  program  are  magnitude  of  PIE  to  be  measured.  This  is  one  of  the  reasons  bleed  was 
removed  from  the  list  of  objectives.  Bleed  is  such  a big  unknown.  Currently  Northrop-Grumman  is  contracted  under 
Boeing  to  look  at  this  problem  using  CFD.  The  size  of  these  holes  are  so  small,  that  it  might  not  be  possible  to  model 
them  at  the  wind  tunnel  scale  talked  about  here. 
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Schedule 


1593 


The  schedule  put  together  for  this  project  shows  some  critical  milestones.  One  is  the  program  feasibility  assessment. 
This  is  where  the  group  decides  if  the  nacelles  are  scaled  and  designed  properly  to  be  able  to  model  the  PIE  items  of 
interest  via  a pilot  nacelle.  If  successful,  the  pilot  nacelle  would  be  one  of  the  four  nacelles.  If  not,  changes  would  be 
made  for  the  4 nacelles  to  be  built.  After  this  milestone  the  model  will  be  fabricated,  nacelles  will  be  calibrated, 
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ii 


Cost  Estimate  & Phasing 


1595 


In  addition  to  the  1 .75  million  dollars  the  model  will  cost,  there  would  be  additional  costs  for  industry  planning  support, 
pilot  nacelle,  balance  fabrication  and  calibration,  nacelle  calibrations,  industry  engineering  test  support,  and  model  and 
instrumentation  test  support.  The  total  cost  would  be  at  4 million. 
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Facilities 
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Not  included  in  this  cost,  is  tunnel  entry  cost.  If  the  NASA  Ames  11x11  foot  Transonic  Wind  Tunnel  does  not  come  up  in 
time  for  the  test,  testing  could  be  done  at  AEDC  16T  at  an  estimated  cost  of  3.3  million.  Supersonic  testing  could  be 
done  at  16S  at  an  estimated  cost  of  3 million,  if  the  NASA  Lewis  10x10  foot  Supersonic  Wind  Tunnel  is  not  available. 
Test  location  is  all  an  unknown  right  now.  The  plan  is  to  do  the  transonic  testing  in  the  Ames  11x11,  since  this  is  what 
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What’s  Next? 
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So,  what  is  next?  A meeting  was  held  in  Seattle  between  the  group  on  February  19.  The  group  was  informed  of  the 
outcome  of  the  CA  review.  The  group  went  to  check  out  the  Boeing  Calibration  Facility  for  the  nacelle  calibrations. 
Individual  members  were  going  to  work  out  details  with  MicroCraft  on  the  next  step  of  the  model  design. 
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Program  Evolution 


1601 


The  program  will  evolve  from  the  program  planning  stage,  working  it’s  way  to  program  execution  and  completion  over 
the  next  2 years. 
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Program  Planning 
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Program  Feasibility 
Assessment 


The  program  planning  includes  the  preliminary  model  design  which  is  done,  putting  together  a test  matrix,  selecting 
equipment,  test  facilities,  identifying  test  methodologies,  supporting  CFD  and  simulation  requirement.  All  these  items 
will  feed  into  the  program  feasibility  assessment  along  with  the  output  of  the  pilot  nacelle  calibration. 


Q. 

| 3 

g>8 
p ® 

&8 

CO  D 

S w 

I-  CO 


1604 


it 


Inlet  Spillage  Drag  Predictions  using  the  AIRPLANE  Code 
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Mach/vector  and  pressure  contours. 

Oscillatory  nature  of  the  solutions. 

Plots  of  lift  and  drag  as  functions  of  the  mass  flow  ratio. 
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AIRPLANE 
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Transonic  Spillage:  AIRPLANE  Mesh 

Wing/Body/Nacelles/Centerbodies/Plugs 
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outboard  -1.196  2.401 

so  the  inboard  nacelle  is  nose-down  slightly,  and  outboard  is  nose-up  with  respect  to  the  freestream  velocity  vector. 
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Transonic  Spillage:  AIRPLANE  Mesh 
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Transonic  Spillage 
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Transonic  Spillage:  AIRPLANE  Solution 

Wing  and  Inboard  Centerbody,  Nacelle,  and  Plug 
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Transonic  Spillage:  AIRPLANE  Solution 
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Transonic  Spillage:  AIRPLANE  Solution 
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Transonic  Spillage:  AIRPLANE  Solution 
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Transonic  Spillage:  AIRPLANE  Solution 
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Transonic  Spillage:  AIRPLANE  Solution 

Plug  Position  X=15 
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Outboard  Nacelle,  Iteration: 


Transonic  Spillage:  AIRPLANE  Solution 

Four  Plug  Positions 
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Transonic  Spillage:  AIRPLANE  Solution 

Four  Plug  Positions 
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Transonic  Spillage:  AIRPLANE  Solution 
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Iteration 


Transonic  Spillage:  AIRPLANE  Solution 
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where  (rho*V)  is  the  local  momentum  vector  and  (dA)  is  the  vector  associated  with  each  element  of  area  on  Si. 
rhoinf*Vinf  is  the  magnitude  of  the  freestream  momentum  vector  and  Aref  is  a reference  area.  In  the  present  case  this  i 
the  area  of  the  inlet,  assumed  to  be  a circle  of  radius  28.3854  inches,  so  Aref  = 2531 .2784  square  inches. 
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Transonic  Spillage:  AIRPLANE  Mesh 

At  4 Degrees  Pitch:  Magenta  Faces  Forward,  Blue  Faces  Aft 
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Transonic  Spillage:  AIRPLANE  Mesh 
At  4 Degrees  Pitch:  Magenta  Faces  Forward,  Blue  Faces  Aft 
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Transonic  Spillage:  AIRPLANE  Solution 
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Transonic  Spillage:  AIRPLANE  Solution 
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Transonic  Spillage:  AIRPLANE  Solution 

Cp  on  a Chordwise  Nacelle  Cut  at  ZCUT=203 
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Transonic  Spillage:  AIRPLANE  Solution 

Cp  on  a Chordwise  Nacelle  Cut  at  ZCUT=359 
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Transonic  Spillage:  AIRPLANE  Solution 
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Transonic  Spillage:  AIRPLANE  Solution 

Drag  as  a Function  of  MFR 
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Transonic  Spillage:  AIRPLANE  Solution  Mach  0.9,  Alpha  4.0 
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Transonic  Spillage:  AIRPLANE  Solution 
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Transonic  Spillage:  AIRPLANE  Solution  Mach  0.9,  Alpha  4.0 
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Use  of  CFD  Results  in  the 
Excrescence  Drag  Estimation 

Sasan  Yaghmaee 
HSCT  Aerodynamics 
Boeing  Commercial  Airplane  Group 
Seattle,  WA 


1997  HSR  Aerodynamic  Performance  Workshop 
February  25-28,  1997 

NASA  Langley  Research  Center,  Hampton,  VA 


A Navier-Stokes  based  method  has  been  developed  to  assess  the  excrescence 
drag  contribution  to  the  flight  polar  build-up  process.  A major  step  in  the 
method  is  the  estimation  of  the  local  boundary  layer  thickness  and  edge 
conditions  from  the  calculated  CFD  solution.  Initial  application  of  the  method 
has  been  successful  in  identifying  drag  critical  areas  of  the  surface,  where 
deviations  from  the  aerodynamically  designed  smooth  sealed  surface  should  be 
avoided.  The  drag  of  a generic  excrescence  item  is  also  presented.  The 
developed  method  can  be  used  for  accurate  and  timely  assessment  of 
weight/drag  trade-offs  for  manufacturing  concepts.  The  detailed  boundary 
layer  data  extracted  from  the  viscous  Navier-Stokes  solution  can  be  a valuable 
tool  in  understanding  fluid  dynamics. 
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H 

SR  Technology  Integration  Airplane  Desi 
Process  Study 

ign 

TI  Task  20 

Airplane  Design  Process 
High  Speed  Drag 

NASA-Industry  "Flight  Polar" 

Linear  Design/Analysis  Code  Comparisons  Build-up  Process 

(MDA)  (BCAG) 

r 

-j  Benchmark  WT  Polars 

-T  WDES/AWAVE  ! 

j '1 

I" ■ 

j Benchmark  Data  Corrections 

-j  BRISTOW/HABS  j 

-|  Excrescence  Methods  & AssessmenJ 

1 

■j  Assessment  of  PD  Design  Assumptions  j 

This  presentation  reports  on  the  progress  made  in  developing  a methodology  to 
assess  excrescence  drag  contributions  to  flight  polar  build-up  process  in  the 
Technology  Integration  (TI)  task  20.  The  flight  polar  build-up  sub-task  was 
assigned  to  Boeing  Commercial  Aircraft  Group  (BCAG). 
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Definition  of  Excrescence  Drag 


The  excrescence  drag  is 

• defined  as  the  sum  of  all  deviations  from  an 
aerodynamically  smooth  sealed  external  surface 

• minimized  by  good  manufacturing  techniques  and 
design  practices 

• generally  due  to  surface  imperfections  smaller  than 
local  boundary  layer  thickness 

• traditionally  normalized  by  the  total  skin  friction 
drag  of  the  airplane 


Technically,  excrescence  drag  is  the  additional  drag  on  the  airplane  due  to  the 
sum  of  all  deviations  from  a smooth,  sealed,  aerodynamically  designed  external 
surfaces.  As  defined,  accurate  manufacturing  techniques  can  significantly 
reduce  the  excrescence  drag.  Design  can  also  reduce  excrescence  drag  by 
proper  placement  of  non-aerodynamic  surfaces. 

The  excrescence  drag  is  generally  due  to  surface  imperfections  that  are  much 
smaller  than  the  boundary  layer  thickness.  The  drag  of  an  imperfection  is  often 
represented  as  a function  of  the  local  skin  friction.  Thus,  to  correlate  the  data  on 
various  airplanes,  the  excrescence  drag  is  traditionally  normalized  by  the  total 
skin  friction  of  the  aircraft. 
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Sources  of  Excrescence  Drag 

• The  excrescence  drag  arises  from 

• Discrete  Items;  antennas,  masts,  and  lights 

• Surface  mismatch;  uneven  manufacturing  joints 

• Leakage  and  internal  flow;  gaps,  holes,  and  seals 

• Surface  roughness  and  waviness;  fasteners  and 

pressurization  “pillowing” 

• The  first  results  from  operational/certification  requirements 

• The  last  three  result  from  structural/manufacturing  trade-offs 


The  sources  of  excrescence  drag  can  be  grouped  into  four  categories,  drag  of 
discrete  items,  drag  due  to  surface  mismatch  at  manufacturing  joints,  drag  of 
surface  roughness,  and  drag  due  to  leakage.  Of  the  four  categories,  the  first 
usually  arises  from  operational/certification  requirements,  while  the  last  three 
result  from  manufacturing  cost/weight/drag  trade-offs. 

Discrete  items  refer  to  antennas,  masts,  etc.  that  are  exposed  to  the  external 
stream.  The  drag  of  the  discrete  items  is  determined  by  their  size,  shape, 
location,  and  orientation  on  the  airplane.  A detailed  knowledge  of  the  external 
flow  during  the  design  phase,  can  be  used  to  optimize  the  placement  of  the  item 
and  in  some  cases  the  shape  of  the  protuberance. 

Surface  mismatch  refers  to  the  steps  that  occur  at  skin  joints,  around  windows, 
control  surfaces,  access  panels,  etc.  A drag  force  is  associated  with  a pressure 
on  a forward  facing  area  and/or  suction  on  a backward  facing  area. 

The  internal  flow  caused  by  leaks  reduces  the  available  useful  energy  and  causes 
additional  drag  force. 

Surface  roughness  refers  to  the  effects  of  distributed  disruptions  related  to 
fastener  flushness  and  surface  finish  smoothness  and  waviness.  The  effect  is 
generally  to  increase  the  skin  friction  force  caused  by  the  air  flowing  over  the 
airplane  surface. 
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Cruise  Excrescence  Drag  Trends 


Airplane  Wetted  Surface 
Area  1000  ft2 


The  Excrescence  contribution  to  total  drag  decreases  with  an  increase  in  wetted 
surface  area  for  subsonic  transports.  The  data  suggests  an  excrescence  level  of 
7%  of  skin  friction  for  large  transports.  A detailed  estimate  of  the  excrescence 
drag  of  the  1970  US-SST  2707-300  configuration  led  to  a 6%  level,  which  is 
slightly  better  than  today’s  subsonic  transport  trend. 

The  2707  estimates  already  included  some  of  the  design  details  for  low 
excrescence  that  we  currently  assume  in  our  “technology  projection”, 
confirming  that  7%  is  a reasonable  excrescence  level  for  a current  technology 
“status”  airplane. 
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Excrescence  Drag  Estimation 


The  excrescence  drag  of  a protuberance  is  traditionally 
defined  as 

CD  = CD*  (A*/Sref)  q«/qi  m 

CD*  Measured  drag  of  isolated  excrescence  (Empirical) 
A*  Area  used  in  CD* 

Sref  Wing  surface  reference  area 

qe  The  equivalent  q experienced  by  the  protuberance 
q[  Local  boundary  layer  edge  velocity 

m Drag  magnification  factor  correcting  for  the  variation 

of  the  local  q from  that  of  the  free  stream 


The  excrescence  drag  of  a discrete  protuberance  is  traditionally  defined  as,  see 
reference  1, 


where  CD^  is  the  empirical  data  on  the  drag  of  individual  isolated 
protuberances,  and  A^  is  the  area  used  in  CD  Tabulation.  S is  wing  surface 
reference  area.  The  factor  q^  q]  corrects  for  the  average  dynamic  head 
experienced  by  excrescence  on  the  surface  when  CD  is  measured  for  the 
isolated  excrescence  exposed  to  a uniform  stream.  For  many  excrescence 
geometries  of  interest  the  drag,  CD  , is  measured  in  presence  of  flat  plate 
turbulent  boundary  layer.  In  this  case  the  dynamic  head  correction  is  set  to 
unity.  The  factor  m is  known  as  the  magnification  factor,  which  corrects  for  the 
variation  of  the  local  qj  from  that  of  the  free  stream. 
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Equivalent  q over  an  Excrescence 


The  equivalent  q is  defined  as  the  average  dynamic  head  experienced  by  the 
protuberance.  Assuming  a power  law  velocity  profile  and  corresponding 
temperature  distribution  for  the  compressible  turbulent  boundary  layer,  the 
definition  becomes 
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where  y is  the  ratio  of  specific  heats,  8 is  the  boundary  layer  thickness,  h is  the 
height  of  the  protuberance,  Me  is  the  local  boundary  layer  edge  Mach  number,  n 
is  the  power  law  exponent,  and  r is  the  recovery  factor. 

The  equivalent  q is  a function  of  the  power  law  exponent  n and  rMe2.  The 
exponent  n is  between  7 and  10  depending  on  the  Reynolds  number.  For 
typical  HSCT  flight  Reynolds  numbers  of  about  200  million  the  appropriate 
value  is  10. 
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The  Equivalent  q 


The  figure  above  shows  the  variation  of  effective  dynamic  head,  q^q,,  with  the 
protuberance  height,  h/8.  The  value  asymptotes  to  one  as  h/8  increases.  The 
assumptions  that  have  gone  into  the  above  universal  curve  are  based  on  the  flat 
plate  turbulent  boundary  layer  theory.  This  should  approximate  the  local 
boundary  layer  flow  about  the  slender  TCA  geometry  at  low  a lift  level  with 
reasonable  accuracy. 

An  alternative  approach  is  to  numerically  integrate  the  product  of  calculated 
CFD  velocity  squared  and  density  profile  at  a surface  location  for  a given 
protuberance  height.  This  method  may  be  more  accurate,  however,  it  is 
considerably  more  time  consuming  and  is  not  suited  for  preliminary  design 
activities. 
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The  Magnification  Factor 

• Accounts  for  the  pressure  gradients  from  excrescence 
location  to  the  downstream  far  field 

• Scales  the  effects  of  the  excrescence  on  the  subsequent 
viscous  flow  on  the  surface 

• Equals  unity  for  flat  plate  flow 

• For  compressible  flow  is  defined  by 
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The  magnification  factor  approximates  the  effects  of  the  excrescence  drag  on  the 
subsequent  development  of  the  flow  on  the  surface.  More  specifically,  the  drag 
of  the  excrescence  significantly  effects  the  subsequent  development  of  the 
momentum  thickness  associated  with  the  boundary  layer  on  a surface  that 
experiences  a non-uniform  pressure  field  . Nash  and  Bradshaw  (Ref.  2 ) 
developed  an  approximation  for  this  effect  for  two  dimensional  incompressible 
flow  by  representing  the  effect  of  excrescence  by  a sudden  change  in  momentum 
thickness.  For  compressible  flow  (Ref.  3)  the  approximation  for  the 
magnification  factor  is 
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Reference  4 compared  the  predictions  based  on  the  compressible  form  with 
subsonic  experimental  data  and  concludes  the  for  sub-critical  conditions  the 
difference  is  within  10%.  However,  for  supercritical  conditions  the  comparison 
significantly  deteriorates  due  to  the  presence  of  shock  waves.  For  transonic 
flows,  Reference  5,  limits  the  validity  of  the  approximation  to  shock  free 
conditions.  The  range  of  validity  of  the  approximate  form  for  supersonic  flows 
has  not  been  established  yet.  However,  it  is  expected  that  outside  of  the  non- 
linear transonic  flow  range,  where  small  perturbations  can  have  significant 
effects  on  the  inviscid  pressure  distributions,  the  approximation  should  be  valid. 
This  assumption,  however,  needs  to  be  verified. 
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CFD-Based  Estimation  of  B.L.  Parameters 

• Boeing  Aerodynamics  developed  post-processor 

• Uses  TLNSMB  calculated  flow  field 

• Estimates  boundary  layer  edge  from  vorticity  profile 

• Calculates  detailed  boundary  layer  data 

- The  local  edge  conditions;  density,  velocity 

- The  displacement  and  momentum  thickness’ 

• Provides  an  alternative  drag  calculation  method 

• Improves  understanding  fluid  dynamics 

• The  detailed  data  facilitates  excrescence  drag  calculation 


Boeing  Aerodynamics  Research  has  developed  a program  to  post-process  the 
TLNSMB  calculated  flow  field.  The  post-processor  judges  the  boundary  layer 
edge  based  on  the  local  vorticity  profile  and  calculates  detailed  boundary  layer 
data.  The  data  contains  the  local  boundary  layer  edge  conditions,  such  as 
density  and  velocity  vector,  as  well  displacement  and  momentum  thickness’. 

The  detailed  data  is  used  to  provide  an  alternative  drag  calculation  method  to  the 
traditional  force  integration.  A further  benefit  is  improved  understanding  of  the 
flow  field  and  increased  confidence  in  the  calculated  drag. 

The  detailed  data  made  available  by  the  post-processor  can  be  used  in  the 
excrescence  drag  calculation. 
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TLNSMB  Calculated  Surface  Pressures 

TCA  Configuration  at  Supersonic  Cruise 


This  plot  shows  the  TLNSMB  calculated  contours  of  surface  pressure  for  TCA 
configuration  at  supersonic  cruise  Mach  number  of  2.4  and  3.5°  angle  of  attack. 
The  salient  features  are  the  glancing  shock  on  the  inboard  upper  surface  of  the 
wing,  and  the  leading  edge  shock  on  both  surfaces  of  the  outboard  wing. 
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Estimated  Boundary  Layer  Thickness  from 
TLNSMB  Solution 

Full  Scale  TCA  Configuration  at  Supersonic  Cruise 


This  plot  shows  the  contours  of  estimated  boundary  layer  thickness  for  the  TCA 
configuration  at  supersonic  cruise  Mach  of  2.4  and  3.5°  angle  of  attack.  The 
calculations  were  made  at  wind  tunnel  Reynolds  number  of  6 million  per  mac 
to  flight  Reynolds  of  number  200  million  per  mac.  The  dimensions  are  in  inches 
full  scale  . For  comparison,  a flat  plate  of  length  of  the  mean  aerodynamic 
chord  of  1 140  in.  would  have  an  estimated  boundary  layer  thickness  of  15  in. 

The  boundary  layer  on  the  body  crown  is  seen  to  be  very  thick.  This  is  a 
classical  feature  of  conical  fore-body  flows  at  an  angle  of  attack,  where  the 
inviscid  pressure  differential  pushes  the  retarded  boundary  viscous  layer  from 
the  keel  to  crown. 
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This  plot  shows  the  contours  of  estimated  magnification  factor  for  the  TCA 
configuration  at  supersonic  cruise  Mach  of  2.4  and  3.5°  angle  of  attack.  The 
magnification  factor  is  seen  to  be  close  to  unity  for  the  most  part.  This  results 
from  the  slender  TCA  geometry  being  at  low  a lift  level,  that  leads  to  small 
deviations  from  a flat  plate  flow. 
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Initial  CFD-Based  Excrescence  Estimations 

• The  method  has  been  used  to 

• Identify  excrescence  critical  regions; 

Larger  values  of  (m  q^) 

• Establish  the  allowable  “sand  grain” 
roughness  height 

• Estimate  the  drag  of  a generic  excrescence; 
A block  of  0.1x6x12  in. 


The  method  has  been  used  to  identify  excrescence  critical  regions,  establish 
allowable  “sand  grain”  roughness  height,  and  estimate  the  drag  of  a generic 
excrescence. 

The  excrescence  critical  region  is  judged  as  the  areas  where  the  term  “m  q^q/’  is 
large,  since  the  term  scales  the  isolated  drag.  The  q^q,  component  of  the  term 
depends  on  the  height  of  the  excrescence,  and  a height  of  0.1  in.  was  used  in  the 
initial  calculations. 

This  allowable  height  is  the  limit  below  which  the  surface  would  be  classified  as 
hydraulically  smooth.  This  height  is  of  the  order  of  laminar  sub-layer  and  is 
defined  by 

y+  = h v*/v  < 5 

where  v*  = (t/p)0-5  , x is  shear  stress,  and  p is  the  density.  All  fluid  properties 
are  evaluated  at  wall. 

To  validate  the  calculation,  the  drag  of  a 0.1  in.  high  block  was  estimated.  The 
block  was  6 in.  long  and  12  in.  wide.  For  the  sake  of  excrescence  drag 
calculation,  the  block  was  represented  by  the  sum  of  a forward-facing  and  an 
aft-facing  steps  of  0.1”  height.  The  experimental  data  on  the  isolated  drag  of 
these  geometries  in  supersonic  flow  are  found  in  Ref.6. 
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Assessment  of  Excrescence  Critical  Regions 

TCA  Configuration  at  Supersonic  Cruise 


High  drag 


Low  drag 


This  plot  shows  the  contours  of  the  quantity  (m  the  TCA  configuration  at 

supersonic  cruise  Mach  of  2.4  and  3.5°  angle  of  attack  for  a generic  excrescence 
of  0. 1 in.  height.  The  higher  levels  of  drag  define  the  critical  regions.  The 
outboard  wing  is  identified  as  excrescence  critical.  On  the  upper  surface  of  the 
wing,  the  area  upstream  of  the  shock  wave  is  also  excrescence  critical. 

The  low  drag  area  near  in  the  wing  body  juncture,  especially  near  the  leading 
edge  on  the  upper  surface,  should  be  viewed  with  suspicion,  since  the 
assumption  of  excrescence  estimation  may  be  grossly  violated.  The  juncture 
flow  is  highly  three  dimensional.  At  the  leading  edge,  the  introduction  of  an 
excrescence  may  effect  the  shock  propagation  and  associated  wave  drag. 
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The  Allowable  “Sand  Grain”  Roughness  Height 

Full  Scale  TCA  Configuration  at  Supersonic  Cruise 


This  plot  shows  the  contours  of  estimated  allowable  sand  grain  roughness  height 
for  the  TCA  configuration  at  supersonic  cruise  Mach  of  2.4  and  3.5°  angle  of 
attack.  The  admissible  roughness  height  is  seen  to  about  0.003  in.  outside  the 
leading  edge.  This  is  consistent  with  estimations  based  on  flat  plate  turbulent 
boundary  layer  flows  ( at  a distance  of  600  in.  admissible  roughness  is  .004  in.). 
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Drag  Sensitivity  to  0.1”  High  Block 

TCA  Configuration  at  Supersonic  Cruise 


0.005  counts 


0.004 

0.003 

0.002 

0.001 


This  plot  shows  the  estimated  the  drag  of  a 0.1x6x12  in.  block  centered  at  a 
location  on  the  TCA  configuration  at  supersonic  cruise  of  Mach  2.4  and  3.5° 
angle  of  attack.  This  generic  excrescence  will  add  about  a thousandth  of  a count 
per  foot  in  the  non-critical  region. 
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Summary  and  Future  Work 


• Progress  made 

• Validated  CD^  value  of  7%  of  CFas  reasonable 

• Developed  a Navier-Stokes  based  excrescence  drag 
estimation  process 

- Identified  the  excrescence  critical  regions 

- Estimated  allowable  “sand  grain”  roughness  height 

- Estimated  the  drag  of  a generic  excrescence 

• Future  work 

• Estimate  drag  due  to  surface  skin  waviness  and 
pressurization  pillowing 

• Utilize  in  weight  and  drag  trade-offs 

• Update  of  technology  projections  for  excrescence  drag 


A Navier-Stokes  based  method  has  been  developed  to  assess  the  excrescence 
drag  contribution  to  the  flight  polar  build-up  process.  Initial  application  of  the 
method  has  been  successful  in  identifying  drag  critical  areas  of  the  surface, 
where  deviations  from  the  aerodynamically  designed  smooth  sealed  surface 
should  be  avoided.  The  drag  of  a generic  excrescence  item  is  also  presented. 
There  has  been  much  progress  on  this  sub-task,  however,  there  are  a few  items 
left  to  complete  the  study.  The  drag  due  to  surface  skin  waviness  and 
pressurization  pillowing  will  be  estimated  in  the  near  future.  The  method  will  be 
utilized  in  weight  and  drag  trade-offs  in  evaluation  of  structural  and 
manufacturing  concepts.  As  the  drag  of  typical  excrescence  items  are  estimated, 
an  update  on  the  technology  projections  for  excrescence  drag  will  be  made. 
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ABSTRACT 

This  paper  presents  an  alternative  method  for  the  calculation  of  wind  tunnel 
model  trip  drag  increments  for  application  to  the  build-up  of  full  scale  airplane 
flight  polars.  Under  Technology  Integration  Task  20,  a benchmark  data  set  of 
Reference  H wind  tunnel  data  was  created  for  test-theory  comparisons  and 
increments.  Several  corrections  were  applied  to  increase  the  fidelity  of  the 
benchmark  data,  one  of  which  was  a correction  for  model  trip  disk  drag. 

A technique  for  the  estimation  of  model  trip  drag  is  derived  based  on  traditional 
excrescence  drag  methods.  The  method  relies  on  both  the  experimentally 
determined  drag  coefficient  of  a cylinder  and  CFD  estimations  of  the  boundary 
layer  thickness  and  magnification  factor.  For  the  proper  application  to  full  scale 
airplane  predictions,  a laminar  flow  correction  is  also  derived  and  incorporated. 
Validation  of  the  proposed  method  is  done  by  comparison  to  the  Braslow  variable 
roughness  size  method  using  wind  tunnel  test  data  from  both  the  NASA  2.7% 
Reference  H model  and  preliminary  data  from  a recent  1.7%  TCA  test. 

Agreement  between  experimental  and  excrescence  methods  is  best  at  the  baseline 
trip  height  of  both  models,  but  diverge  as  trip  height  is  increased.  Possible 
sources  for  the  differences  include:  test  data  uncertainty,  curve  extrapolation,  data 
applicability,  magnification  factor  estimate,  and/or  boundary  layer  characteristic 
estimates.  Recommendations  are  made  for  additional  work  to  validate  the 
excrescence  and  other  trip  drag  estimation  methods  at  other  test  conditions  and 
suggestions  are  made  for  additional  checks  and  assessments. 
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Presentation  Outline 
— 

• Objective  - Technology  Integration  Task  20 

• Review  of  trip  drag  correction  methodology 

• Excrescence  method  description 

• Sample  results 

• Conclusions 

• Recommendations  for  future  work 


The  content  of  this  presentation  is  organized  in  the  following  manor.  First  a 
description  of  the  1996  Technology  Integration  task  under  which  the  current 
work  was  performed  is  reviewed.  Next,  a brief  discussion  of  the  current 
experimental  methods  available  for  determining  trip  drag  is  covered,  reviewing 
the  pro’s  and  con’s  of  each.  The  emphasis  of  the  presentation  is  then  devoted  to 
the  review  of  excrescence  methodology  and  its  application  to  the  problem  of  trip 
drag.  Each  term  in  the  excrescence  equation  is  reviewed  in  detail  and  sample 
results  are  presented.  Based  on  the  results  of  this  exercise  important  conclusions 
and  insights  are  drawn  from  the  data  and  possible  sources  of  error  are  identified. 
Finally,  conclusions  about  the  significant  findings  of  the  study  are  made  and 
recommendations  for  future  work  are  identified. 
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The  FY  1996  Technology  Integration  Airplane  Design  Process  Study  (Task  20) 
was  composed  of  two  separate  sub-tasks.  The  first  was  a comparison  of  linear 
design  and  analysis  codes  conducted  by  McDonnell  Douglas  Aerospace.  For  this 
task,  MDA  compared  the  capabilities  of  the  three  major  linear  theory  analysis 
codes  used  by  NASA,  Boeing,  and  MDA.  The  process  involved  design  of  a 
linearly  optimized  TCA  and  then  a reanalysis  and  evaluation  of  the  resulting 
configuration  in  each  of  the  three  codes  under  consideration. 

The  second  sub-task,  conducted  by  Boeing  Commercial  Airplane  Group,  was  a 
detailed  assessment  and  improvement  of  the  prediction  process  of  full  scale  flight 
performance  of  the  TCA.  To  calibrate  these  prediction  methods,  a single  set  of 
corrected  wind  tunnel  data  was  created  to  serve  as  a common  “benchmark”  for 
calibration  of  prediction  methods.  Several  corrections  were  made  to  the 
benchmark  data  to  account  for  wind  tunnel  to  flight  differences,  one  of  the  most 
significant  being  a correction  for  model  trip  drag. 

TCA  wind  tunnel  data  was  not  available  until  late  in  fiscal  1996,  therefore  the 
benchmark  data  for  this  task  is  based  on  Reference  H data  taken  between  1993 
and  1995.  ADP  work  at  Boeing  in  FY97  will  include  updating  TTs  performance 
polars  for  refined  technology  projections  and  TCA  wind  tunnel  test  data,  as  well 
as  improving  the  capability  of  predicting  the  drag  of  multi-surface  control 
concepts. 


Corrections  to  Benchmark  Wind  Tunnel  Data 


Wind  Tunnel  Data  Corrections: 

• Model  Forces 

• Model  Cavity  Pressures 

• Nacelle  Base  Pressures 

• Nacelle  Internal  Forces 

Missing  Corrections: 

• Aeroelastics/Reynolds  No. 

• Model  Fidelity 

• Model  Trip  Drag 

• Leading  Edge  Laminar  Flow 


Before  wind  tunnel  results  can  be  applied  to  full  scale  flight  performance 
predictions,  the  benchmark  data  tables  must  be  corrected  for  all  significant  wind 
tunnel  induced  effects.  This  step  is  required  because  the  limitations  in  either 
analysis  capability  or  testing  capability  (model  fidelity,  test  flexibility,  facility 
Reynolds  number  capability,  etc.)  often  introduce  differences  between  the  tested 
model  and  the  analyzed  configuration.  Force  data  from  a test  is  corrected  for 
angle  of  attack,  cavity  and  nacelle  base  pressures,  and  nacelle  internal  forces. 
Corrections  not  applied  to  test  data,  but  required  for  full  scale  airplane 
performance  predictions  include: 

• Model  aeroelastics/Reynolds  number 

• Model  fidelity  and  facility  selection 

• Effect  of  boundary  layer  tripping  devices 


To  more  accurately  represent  the  flow  experienced  by  the  full  scale  vehicle  at 
flight  Reynolds  number,  it  is  common  practice  in  wind  tunnel  testing  to 
artificially  transition  the  flow  to  fully  turbulent  using  a boundary  layer  trip.  In 
recent  years,  the  trend  in  testing  has  been  to  use  transition  disks  made  of  epoxy 
resin  instead  of  the  traditional  “grit”  composed  of  finely  screened  sand  particles. 
Trip  disks  are  easier  to  apply,  measure,  and  maintain  and  are  more  durable  than 
their  grit  predecessors.  Their  disadvantage  is  that  they  introduce  a significant 
drag  increment  into  the  data  which  to  date  has  never  been  satisfactorily 
quantified.  The  goal  of  this  study  is  to  quantify  the  trip  drag  of  transition  disks 
using  excrescence  methods.  A proper  accounting  for  leading  edge  laminar  flow  is 
also  necessary  to  correctly  evaluate  trip  drag. 
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Trip  Drag  Correction  Methods  and  Issues 

— — 


• Variable  roughness  size  Cd 

Baseline  Trip>^  x 

- Curve  fit  A CD 

_ _ - - -*■  * 

- Applicable  data  range  ~ 

~ 1 Rn=C 

- Extrapolation  range  ' 

Trip  Ht. 

• Variable  Reynolds  no.  qd 

Basel  ir«  Trip 

- Facility  selection  AC^’ 

ht.=c 

- Ae roe lastic  corrections 

- Extrapolation  range  - 

Rn 

• Excrescence  calculations 

- Theoretical 

- Cheap  C. 


To  date,  calculation  of  drag  increments  for  transition  disks  on  HSR  models  have 
followed  a methodology  outlined  by  Braslow  (Reference  1)  known  as  the 
variable  roughness  size  method.  This  method  is  based  on  grit  type  trips  and  their 
second  order  variation  of  drag  based  on  the  characteristic  "area"  of  a grit  (k2), 
where  k represents  the  height  (or  width)  of  the  grit.  Lines  (Reference  2) 
describes  a method  similar  to  Braslow's  variable  roughness  size  method  except  he 
proposes  a characteristic  area  (k),  based  on  disk  height  This  linear  variation  in 
drag  with  trip  height  was  proposed  to  better  represent  disk  type  trips  in  transonic 
flow  which  extend  beyond  the  edge  of  the  laminar  boundary  layer. 

A second  method  proposed  by  Braslow  known  as  the  variable  Reynolds  number 
method,  offers  an  alternative  method  for  the  calculation  of  trip  drag.  This  method 
relies  on  extrapolation  of  data  taken  for  a constant  height  trip  over  a large 
Reynolds  number  range  where  the  value  for  fully  turbulent  flow  is  known. 

These  experimental  determinations  of  the  trip  disk  drag  increment  lend  insight 
into  the  magnitude  of  the  trip  drag  correction,  but  remain  inconclusive  because  of 
their  inherent  uncertainties.  Uncertainties  in  the  variable  roughness  method 
include:  which  characteristic  "area"  is  appropriate  for  a cylinder  in  supersonic 
flow  (hence,  which  order  curve  fit  was  most  appropriate  for  the  Cj,  vs.  roughness 
height  curves,)  large  extrapolation  regions,  and  what  range  of  heights  are 
appropriate  to  curve  fit  for  accurate  results.  The  variable  Reynolds  number 
method  is  limited  by  facility,  extrapolation  range,  and  model  loading. 

Calculation  of  trip  drag  using  excrescence  methods  is  proposed  as  an  alternative 
to  the  experimental  methods  above  without  the  aforementioned  uncertainties. 
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Excrescence  Method  - Trip  Drag  Application 

*5^ 

A Cn  =Cn  x^-xmxK 

utrip  udisk  C 

°ref 

= Trip  drag  increment  correction 
= 3D  cylinder  on  flat  plate  in  supersonic  flow  (empirical) 

= “Magnification  factor”  of  turbulence  due  to 
local  Cp  and  boundary  layer  profile 
= Total  cylinder  reference  area/wing  reference  area 

= corrects  for  differences  between  original  experiment 
and  application  (1  for  experimentally  derived  CDm  ) 


The  generalized  excrescence  equation,  discussed  in  Reference  3,  can  be  rewritten 
to  represent  the  incremental  drag  caused  by  wind  tunnel  model  trip  disks  in 
transonic  or  supersonic  flow  and  is  presented  above  as  AC Dtrip.  CDdisk  is  the  drag 
coefficient  of  a cylinder  on  a flat  plate  based  on  local  q and  total  disk  reference 
area,  . Sref  is  the  model  reference  area  and  m is  the  drag  magnification  factor 
due  to  the  non-uniform  pressure  distribution  effects  on  the  boundary  layer 
downstream  of  the  roughness  element. 

The  factor  K,  applied  to  , allows  for  the  differences  between  the  geometry 
and  orientation  of  the  excrescence  as  originally  tested  in  the  wind  tunnel  and  as 
located  on  the  wind  tunnel  model.  If  there  are  no  differences,  then  is  truly 
representative  and  K = 1.  The  original  wind  tunnel  test  experiment  used  to  derive 
the  drag  coefficient  of  a three  dimensional  cylinder  appears  to  properly  capture 
dynamic  pressure  effects  and  should  accounts  for  the  shock  off  of  the  portion  of 
the  cylinder  which  extend  beyond  the  edge  of  the  boundary  layer  into  the 
ffeestream.  Therefore,  for  HSCT  class  wind  tunnel  models  and  the  excrescence 
method  proposed  for  calculating  CDdisk , the  cylinder  data  is  applicable  and  K is 
assumed  to  equal  1. 
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Presented  above  is  CFD  solution  of  the  drag  magnification  factor  for  the 
Reference  H configuration  at  wind  tunnel  conditions.  The  viscous  drag 
magnification  factor,  m,  for  an  excrescence  in  compressible  flow  is  defined  in 
Reference  4 as 


X[l  + 0^Af|  J 


The  subscripts  L,  TE,  and  °°  indicate  local,  trailing  edge,  and  free  stream 
conditions,  respectively.  This  formula  is  a compressible  flow  adaptation  of  the 
Nash  and  Bradshaw  method  (Reference  5)  for  prediction  of  magnification  of 
roughness  drag  by  pressure  gradients  for  incompressible  flow.  This  parameter 
can  be  easily  calculated  using  local  flow  parameters  from  existing  Euler  or 
Navier-Stokes  computational  solutions  for  the  geometry  under  consideration. 

TLNS3D  (multi-block  Navier-Stokes  CFD  code)  was  run  to  calculate  the 
turbulent  flow  magnification  factor-for  both  the  upper  and  lower  surface  of  the 
full  scale  Reference  H configuration  at  wind  tunnel  Reynolds  numbers.  Typical 
turbulent  boundary  layer  magnification  factors  at  the  location  of  the  trip  disks 
ranged  from  0.50  to  0.80. 
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The  generalized  equation  for  CDdisk  presented  above  is  based  on  the 
measurements  made  by  Gaudet  and  Winter  (Reference  6)  and  Palhster 
(Reference  7.)  Their  experiments  measured  the  drag  of  various  forms  of 
excrescence  over  a Mach  number  range  of  0.2  to  2.8.  The  data  is  known  to  be 
valid  for  the  following  conditions:  finite  length  cylinders,  a turbulent  boundary 
layer  ahead  of  the  excrescence,  zero  pressure  gradient,  two  dimensional  flow,  and 
h/8  ranging  from  0.5  to  2.0. 

For  circular  cylinders  which  extend  into  the  free  stream  (h>8)  the  drag  is 
calculated  by  pro-rating  the  submerged  and  exposed  areas  with  the  applicable 
drag  coefficient  and  then  summing  the  two  components.  Cos  and  represent 
the  drag  of  the  finite  and  infinite  length  cylinder,  respectively.  Both  Cm  and 
can  be  obtained  using  plots  created  by  Gaudet  and  Pallister,  given  Reynolds 
number  based  on  cylinder  diameter  (Rd)  and  Mach  number.  These  plots  are 
shown  later  is  this  paper.  The  above  equation  was  derived  for  cylinder  heights 
which  were  equal  to  or  greater  than  the  flat  plate  turbulent  boundary  layer  they 
were  immersed  in.  For  cylinders  immersed  in  the  boundary  layer 

CD„  = CD5  - FC, 

Where  F can  be  determined  from  the  plots  of  Reference  6. 
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Presented  above  is  CFD  solution  of  the  turbulent  boundary  layer  thickness  for  the 
Reference  H configuration  at  wind  tunnel  conditions.  TLNS3D  (multi-block 
Navier-Stokes  CFD  code)  was  run  to  calculate  the  boundary  layer  height  for  both 
the  upper  and  lower  surfaces. 


Typical  2.7%  model  scale  turbulent  boundary  layer  heights  at  the  location  of  the 
trip  disks  range  from  0.02  to  0.033  inches. 
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Laminar  Boundary  Layer  Height  Correction 



• Experimental  determination  of  disk  drag  coefficient  based 
on  turbulent  flow  upstream  of  cylinder 

• Model  disks  have  laminar  flow  upstream 

• Need  to  scale  estimated  5mr*  to  laminar  boundary  layer 


thickness 


For  HSCT  class  wind  tunnel  models,  it  has  been  observed  using  flow 
visualization  techniques  that  because  of  the  low  Reynolds  numbers  and  model 
scales  being  tested,  there  exists  a laminar  boundary  layer  on  the  portion  of  the 
wing  ahead  of  the  trip  disks.  In  order  to  use  the  equation  proposed  for  CM,t  in 
this  analysis  one  significant  assumption  must  be  made.  Calculation  of  CDdisk 
using  the  results  of  Gaudet  and  Pallister  is  applicable  as  long  as  the  laminar 
boundary  layer  height  on  the  leading  edge  of  the  wing  of  the  wind  tunnel  model 
is  approximately  simulated. 

For  this  study,  only  turbulent  boundary  layer  height  estimates  were  available. 
Following  the  above  assumption,  the  CFD  estimate  must  then  be  corrected  to  an 
approximate  laminar  height  by  using  the  ratio  of  the  incompressible,  flat  plate 
laminar  and  turbulent  boundary  layer  height  equations.  The  simplified  form  is 
presented  above. 

It  is  acknowledged  that  the  laminar  boundary  layer  assumption  is  not  entirely 
appropriate  for  the  methodology  outlined  here  because  the  original  data  used  to 
generate  the  drag  coefficient  of  a cylinder  was  taken  in  a turbulent  boundary 
layer.  However,  it  is  proposed  that  proper  modeling  of  the  trip  to  boundary  layer 
height  ratio  is  the  most  important  factor  in  the  CDdisk  equation  in  order  to  correctly 
estimate  the  disk  drag  coefficient.'  Also,  the  primary  source  of  drag  for  an 
excrescence  comes  from  the  pressure  drag  on  the  aft  face  of  the  disk.  Therefore 
using  a laminar  boundary  approximation  would  not  significantly  affect  this  flow 
phenomena  for  the  disk  Reynolds  numbers  tested.  Finally,  no  other  data  has  been 
found  which  directly  expresses  the  supersonic  drag  of  a cylindrical  excrescence 
in  a laminar  boundary  layer. 
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Laminar  Boundary  Layer  Correction  (cont.) 
— — 

• A drag  correction  is  required  to  account  for  the  velocity  profile 
differences  between  laminar  and  turbulent  boundary  layers 

• A 7%  correction  to  CD5  was  derived  from  experimental  data 


To  compensate  for  the  differences  in  velocity  profile  (hence,  dynamic  pressure) 
associated  with  the  use  of  a laminar  boundary  layer  in  place  of  a turbulent  one,  a 
small  correction  factor  should  be  applied  to  the  calculation  of  CDdisk  . This  q 
correction  is  assumed  to  be  proportional  to  the  area  ratio  of  the  two  velocity 
profiles.  Based  on  experimentally  determined  velocity  profile  data  contained  in 
Reference  8 at  Mach=2.40  and  similar  Reynolds  number,  a correction  of 
approximately  7%  was  determined. 

This  correction  may  be  different  for  the  transonic  Mach  numbers  because  both 
boundary  layer  thickness  and  velocity  profile  ratios  will  change.  As  applied 
currently,  this  correction  does  not  vary  with  trip  height  for  trips  inside  the 
boundary  layer  and  is  only  applicable  for  drag  predictions  where  h>8.  A 
formulation  would  have  to  be  derived  before  any  future  trip  drag  predictions  can 
be  made  for  trips  immersed  in  the  boundary  layer  (h<5). 

The  equation  for  is  rewritten  above  to  denote  the  laminar  boundary  layer 
height  and  velocity  profile  correction  factor. 


1678 


Knowing  Mach  number  and  the  Reynolds  number  based  on  disk  diameter  (Rd), 
CD5  and  can  be  obtained  from  the  Figures  above.  With  8^,  and  CD6  and 
now  available,  CDdiii.  can  then  be  calculated. 
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Trip  Drag  Solution  Process 
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To  simplify  the  analysis  about  the  entire  wing,  the  necessary  components  within 
the  excrescence  equation  were  calculated  only  at  the  midpoint  of  each  wing 
segment  (inboard,  mid-wing,  and  outboard)  for  both  the  upper  and  lower  surface. 
The  excrescence  drag  was  then  calculated  at  each  of  the  6 wing  segments  shown 
above  and  the  segment  incremental  trip  drags  were  then  summed  to  represent  the 
trip  drag  for  both  wings.  at  each  segment  represented  only  the  frontal  area 
associated  with  the  number  of  trips  present  in  that  segment. 
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The  comparison  between  excrescence  method  results  and  wind  tunnel  data 
utilizing  the  variable  roughness  height  method  (with  a linear  fit)  for  the 
Reference  H configuration  at  wind  tunnel  conditions  is  presented  above.  A 
significant  offset  (approximately  3 cts.,  average)  can  be  seen  between  the  two 
methods  with  the  excrescence  method  consistently  predicting  trip  drags 
approximately  double  that  of  the  experimental  method.  For  the  baseline  wind 
tunnel  trip  height  (0.01 1 in.)  the  variable  roughness  height  calculates  a trip  drag 
correction  of  2.5  cts.  while  the  excrescence  method  predicts  a correction  of  4.5 
cts.  The  two  methods  show  similarities  in  trip  drag  trend  and  slope,  with  the 
excrescence  method  again  predicting  a slightly  greater  slope. 


Although  there  appears  to  be  little  agreement  between  the  two  methods,  neither 
of  the  two  method  can  be  dismissed  as  “incorrect”  because  of  the  uncertainties 


associated  with  the  experimental  data  and  the  assumptions  used  in  the 
excrescence  method.  Additional  comparisons  are  made  in  following  slides  and  an 
assessment  of  each  of  the  significant  components  of  both  methods  is  also 
presented  later. 
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Presented  above  is  a comparison  of  excrescence  and  variable  roughness  height 
methods  for  a lift  coefficient  corresponding  to  the  minimum  drag  coefficient 
(Cl=0.015)  on  the  Ref.  H drag  polar.  As  observed  previously,  the  excrescence 
method  predicts  a higher  trip  drag  than  the  experimental  method  while  the  slope 
trend  with  h/8  (trip  disk  height/average  boundary  layer  height)  of  the  two  curves 
is  similar. 

As  seen  in  other  trip  drag  assessments,  the  variable  roughness  method  shows  a 
significant  trip  drag  variation  with  changing  lift  coefficient.  When  compared  to 
the  cruise  condition  (CL=0.12,  previous  slide)  at  the  baseline  trip  height,  the 
variable  roughness  height  method  shows  a 0.9  ct.  increase  in  trip  drag  (2.5  cts.  vs. 
3.4  cts.)  The  excrescence  method,  however,  shows  that  trip  drag  has  only  a slight 
dependence,  0. 1 ct,  on  lift  coefficient  (4.5  cts.  vs.  4.6  cts.) 
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Uncertainties  in  Trip  Drag  Methods 

. — . 


Experimental  Method 

• Data  uncertainty 

• Applicable  data  range 

• Extrapolation  order 

• Applicable  extrapolation  range 


Trip  Ht. 


Excrescence  Method 

• Turbulent  magnification  factor 

• Boundary  layer  height  estimate 

• Compressibility  assumptions 


# Model  Wing  Leading  Edge 
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To  fully  understand  the  problems  in  determining  trip  drag  using  either  the 
excrescence  or  variable  roughness  height  methods,  the  uncertainties  of  each 
method  must  be  carefully  reviewed.  The  primary  uncertainties  in  each  method 
are  listed  above  and  explained  below. 

Variable  Roughness  Height  Method 

• There  exists  scatter  in  each  data  point  taken  due  to  tunnel  repeatability  issues. 
Extrapolation  to  zero  trip  height  can  vary  on  the  order  of  ~±0.5  ct.  due  to  data. 

• No  firm  guidelines  exist  for  selecting  data  points  to  fit  Selection  of  various 
combinations  of  trip  data  can  vary  the  extrapolation  results  on  the  order  of  ~±1  ct. 

• There  is  litde  experience  in  selection  of  order  of  curve  fit  for  trip  disks  in 
supersonic  flow.  The  effect  on  extrapolation  is  on  the  order  of  ~±2  cts. 

• It  difficult  to  understand  how  a curve  fit  of  drag  data  for  disks  extending  into 
the  ffeestream  can  properly  approximate  disks  immersed  in  the  boundary  layer. 
Effect  on  extrapolation:  ~±?  cts. 

Excrescence  Method 

• Boundary  layer  transition  on  a wind  tunnel  model  does  not  immediately  begin 
at  the  trip  location  for  the  entire  model.  Perhaps  correlation  would  improve  with 
a magnification  factor  calculated  for^a  laminar  boundary  layer. 

• Method  fidelity  could  be  improved  by  calculating  directly  the  laminar  boundary 
layer  height  using  CFD  tools. 

• Compressible  formulations  for  the  turbulent  to  laminar  boundary  layer  height 
and  skin  friction  corrections  would  be  more  accurate. 
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Another  comparison  between  methods  is  presented  above  for  the  TCA 
configuration  based  on  preliminary  trip  drag  estimates  for  the  TCA  (UPWT 
1679).  The  comparison  between  the  two  methods  show  similar  trends  as  for  the 
Reference  H configuration.  An  offset  (an  average  of  5 cts.)  can  again  be  seen 
between  the  two  methods  with  the  excrescence  method  consistently  predicting 
trip  drags  approximately  double  that  of  the  experimental  method.  A slight  slope 
difference  can  also  be  seen  between  the  two  methods.  For  the  baseline  wind 


tunnel  trip  height  (0.012  in.)  on  the  TCA,  the  variable  roughness  height  calculates 
a trip  drag  correction  of  5.2  cts.  while  the  excrescence  method  predicts  8.9  cts. 


At  first  glance,  application  of  the  excrescence  method  appears  to  only  further 
cloud  the  trip  drag  issue,  but  this  “analogy”  allows  us  to  begin  to  understand  the 
flow  mechanisms  affecting  trip  drag  (represented  by  each  term  in  the  excrescence 
equation).  To  aid  in  understanding,  the  plot  of  the  2.7%  Ref.  H calculation  is  also 
included.  Comparisons  of  the  full  scale  boundary  layer  height  and  magnification 
factor  calculations  between  the  two  models  showed  no  significant  differences. 
Remaining  sources  for  trip  drag  differences  then  only  include  A^et  and  Cp^. 
Because  TCA  has  more  wing  area  than  the  Ref.  H,  it  is  known  that  there  is 
approximately  11%  more  disk  A^  on  the  TCA,  which  can  explain  some  of  the 
increased  drag  over  the  Ref.  H model.  Also,  because  the  TCA  is  a smaller  scale 
model  than  the  Ref.  H,  the  baseline  trip  h/S  ratio  is  greater,  leading  to  a larger 
value  of  CrmA.  and  a higher  trip  drag  located  further  up  the  curve.  If  the  trip  were 
reduced  to  a height  known  to  trip  (h/fel.2),  the  trip  drag  could  be  reduced  to  drag 
ranges  previously  predicted  by  both  methods  for  the  Ref.  H configuration. 


Currently  there  is  no  explanation  for  the  absolute  drag  level  difference  between 
the  two  methods. 
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Leading  Edge  Laminar  Flow  Correction 


• Application  of  wind  tunnel  data  for  prediction  of  full  scale 
flight  performance  polars  requires  a correction  for  model 
leading  edge  laminar  flow 

- Flow  visualization  techniques  have  shown  significant 
regions  of  laminar  flow  on  wind  tunnel  models 


• Approximately  equal  to  2 to  3 cts. 


Right  - Wind  Tunnel  + Laminar  Flow  Correction 

(oofT.  tor  trip  drag) 
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A final  correction  must  be  used  on  the  outlined  methodology  for  accurate 
application  of  wind  tunnel  data  to  full  scale  flight  conditions.  It  has  been  verified 
in  the  wind  tunnel  that  laminar  flow  occurs  ahead  of  the  trip  disks  on  the  model 
over  the  entire  span  of  the  wing  (hence,  the  need  for  trips.)  This  region  of 
laminar  flow  on  the  model  represents  approximately  6.5%  (Reference  H)  to  8.2% 
(TCA)  of  the  projected  wing  reference  area.  This  is  a significant  percentage  of 
the  wing  reference  area  that  should  not  be  neglected. 

Using  the  incompressible  Blausius  and  Prandtl  flat  plate  boundary  layer  skin 
friction  equations  as  an  approximation,  a drag  increment  can  be  calculated  to 
represent  the  turbulent  skin  friction  difference  between  a leading  edge  laminar 
flow  (model)  and  turbulent  flow  (full  scale.)  The  equation  representing  this 
correction  is  shown  above.  The  magnitude  of  this  correction  is  on  the  order  of  2 
to  3 cts. 
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Photograph  of  Laminar  Flow  on  Model 


Presented  above  are  sample  images  of  a sublimation  run  done  on  the  1.7%  TCA 
model  2b  (UPWT  1679)  to  illustrate  the  extent  of  laminar  flow  present  on  typical 
HSCT  wind  tunnel  models.  The  sublimation  technique  shows  areas  of  laminar 
and  turbulent  flow  by  using  the  relative  skin  friction  (scrubbing  power)  of  each 
flow  to  sublime  Flourine  powder  off  of  the  wing’s  surface.  In  this  image  lighter 
colored  areas  represent  laminar  flow  and  grayer  areas  represent  turbulent  flow. 

In  this  example,  the  outboard  wing  transitions  to  turbulent  flow  at  the  trip  disks. 
The  inboard  section  of  the  wing  trips  immediately  at  the  disks  but  takes 
approximately  0.75  inches  to  transition  to  fully  turbulent  flow.  The  two  large 
wedges  of  laminar  flow  seen  in  the  pictures  are  where  trip  disks  have  been 
intentionally  removed  to  observe  natural  transition.  The  black  streaks  behind  the 
inboard  trip  disks  represent  additional  scrubbing  due  to  the  shed  vortices  of  the 
disks  themselves. 
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Summary  of  Corrections  to  Benchmark  Data 




• Excrescence  method  trip  drag  corrections  were  applied  to 
benchmark  Ref.  H data  over  a range  of  Mach  numbers 

- Correction  for  model  leading  edge  laminar  flow  included 

• No  variation  with  CL  was  included 

- "Engineering  judgment"  used  to  fill  in  benchmark  corr.  table 

• Only  a select  number  of  CFD  solutions  were  available 

• Relative  magnitude  of  correction  values  between  Mach 
numbers  compared  against  wind  tunnel  estimates 


Trip  Drag  Correction  Table  - Ref.  H 


Mach 

0.90 

0.90 

(w/flap) 

1.20 

2.40 

Correction 

-0.00008 

-0.00011 

0.0 

-0.00025 

20 


The  methodology  previously  outlined  was  employed  to  generate  a table  of  values 
to  be  added  directly  to  the  benchmark  wind  tunnel  data  to  correct  the  data  for 
drag  due  to  model  transition  disks.  Because  boundary  layer  heights  and 
magnification  factors  had  to  be  taken  from  expensive  Navier-Stokes  CFD 
solutions,  only  a select  number  of  conditions  were  run.  The  table  above  shows 
the  conditions  for  which  the  trip  disk  drag  increments  ware  actually  calculated. 

Corrections  at  the  specified  Mach  numbers  were  calculated  for  the  cruise 
condition  and  all  intermediate  steps  matched  the  tunnel  conditions  of  the 
benchmark  data.  The  corrections  were  assumed  to  be  constant  throughout  the 
range  of  CL.S  covered  by  the  data  and  input  as  such  in  the  tables.  For  Mach 
numbers  or  configurations  within  the  benchmark  data  which  lied  between  the 
calculated  values,  the  data  were  linearly  interpolated  between  the  two  nearest 
applicable  conditions.  There  are  two  exceptions  to  this  statement  First  data  at 
M=0.95  and  0.98  were  input  as  twice  (2x)  the  data  calculated  at  M=0.90.  This 
assumption  was  made  based  on  the  results  of  a test  conducted  by  McDonnell 
Douglas  Aerospace  on  their  model  M2.4-7A  Opt  5 (Reference  9)  in  September, 
1994.  Second,  corrections  for  trip  disk  drag  at  M=1.07  and  1.10  were  assumed  to 
be  the  same  as  that  calculated  at  M=1.20.  Because  the  transonic  assumptions 
were  based  on  the  results  of  a model'M2.4-7A  test,  they  should  be  verified  for  the 
TCA  or  Reference  H configurations  when  data  becomes  available. 

The  application  the  trip  drag  corrections  to  the  benchmark  data  is  shown  below. 
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For  application  to  FY  96  Technology  Integration  Task  20,  the  results  of  the 
previously  outlined  excrescence  methodology  were  used  to  correct  Reference  H 
based  benchmark  wind  tunnel  data  for  trip  drag  effects.  Comparison  of 
excrescence  method  results  to  the  experimentally  based  variable  roughness  size 
method  show  significant  differences,  but  the  uncertainties  in  the  experimental 
methods  raise  many  other  questions.  Until  a better  understanding  of  the  trip  drag 
phenomena  can  be  acquired,  multiple  trip  drag  estimation  methods  should 
continue  to  be  explored. 

Application  of  the  excrescence  method  to  the  problem  of  model  trip  drag 
provides  a valuable  analogy,  based  on  boundary  layer  flow  phenomena,  to  help 
understand  the  mechanisms  involved  in  boundary  layer  transition.  Also,  it  has 
the  potential  to  reduce  the  need  for  expensive  and  lengthy  wind  tunnel  testing  for 
trip  drag  determination.  Finally,  because  the  method  is  flow  phenomena  based,  an 
accurate  solution  should  be  possible  once  all  component  uncertainties  are 
resolved. 

Almost  as  important  as  trip  drag  in  wind  tunnel  to  flight  or  CFD  correlations  is 
the  proper  accounting  for  the  laminar  flow  which  exists  in  wind  tunnel  data. 
Current  predictions/correlations  do  not  account  for  this  drag  component  and 
could  be  in  error  by  as  much  as  2 to  3 counts  at  wind  tunnel  conditions. 

It  is  apparent  that  currently  no  one  method  is  mature  enough  to  fully  quantify  the 
trip  drag  problem,  therefore  additional  work  on  the  topic  is  recommended. 
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Recommendations  for  Future  Work 

*5^" 

• Assess  laminar/turbulent  and  compressible/incompressible 
assumptions  used  in  current  excrescence  method 

- Can  affect  magnification  factor  and  CMisk 

• Detailed  trip  drag  study  in  the  wind  tunnel 

- Quantify  the  effects  of  height,  location,  grit  type, 

Reynolds  number,  planform,  model  scale,  etc. 

• Compare  excrescence  & variable  Reynolds  no.  methods 

- Little  supersonic,  high  Reynolds  number  data  available 

• Perform  more  theoretical/experimental  comparisons 

- Understand  the  drag  of  trips  immersed  in  the  B.L. 

- Use  the  method  to  understand  linear  vs.  quadratic  fits 

- Explore  trip  drag  dependence  on  lift  coefficient 


The  above  items  detail  some  topics  for  additional  work  which  should  be 
completed  to  better  understand  the  trip  drag  phenomena  at  transonic  and 
supersonic  conditions. 


1689 


References 


1.  Brastow.  A.  L,  Kicks.  R M.,  and  Harris  Jr.,  R V.,  Use  of  Grit-Type  Boundary-Layer-Transition  Trips 

on  Wind-Tunnel  Models.  NASA  TN  D-3579.  Sept.,  1966. 

2.  Lines,  T.  R.,  “Aerodynamic  Wind  Tunnel  Test  Summary  • Epoxy  Disk  Trip  Evaluation  (BTWT 

1837),’  Boeing  Internal  Document  Feb.,  1982. 

3.  ESDU,  An  introduction  to  aircraft  excrescence  drag.  Item  No.  90029,  Engineering  Sciences  Data 

Unit  London,  March  1993. 

4.  ESDU,  Calculation  of  excrescence  drag  magnification  due  to  pressure  gradients  at  high  subsonic 

speeds,  item  No.  87004,  Engineering  Sciences  Data  Unit  London,  November  1969. 

5.  Nash,  J.  F.,  and  Bradshaw,  P.,  The  Magnification  of  Roughness  Drag  by  Pressure  Gradients, 

Journal  of  the  Royal  Aeronautical  Society,  January,  1967. 

6.  Gaudet  L.  Winter,  K.  G.,  Measurements  of  the  Drag  of  Some  Characteristic  Aircraft  Excrescences 

Immersed  in  Turbulent  Boundary  Layers,  RAE  TM  Aero  1538,  Bedford,  October,  1973. 

7.  Paibster.  K.  C..  Wind  Tunnel  Measurements  of  the  Transonic  Drag  of  Excresoences  Immersed  in  a 

Turbulent  Boundary  Layer,  ARA  #37.  Bedford,  December,  1974. 

8.  ODonrteN,  R M.,  Experimental  Investigation  at  a Mach  Number  of  2.41  of  average  SWn-Frictton 

Coefficients  and  Velocity  Profiles  for  Laminar  and  Turbulent  BoSundary  Layers  and  an 
Assessment  of  Probe  Effects,  NASATN  3122,  Jan,  1954. 

9.  Cheung.  Samson,  Narducd,  Robert  and  Agrawat  Shreekarrt,  Post-Test  Data  Analysts  for  the 

1.675%  M2.4-7A  Opts  Model  In  the  NASA  Langley  16*  Transonic  Wind  Tunnel,  NASA  Contract 
NAS 1-20220,  Task  3,  January  21, 1995. 


23 


The  references  listed  above  were  used  in  the  development  of  the  trip  disk 
excrescence  methodology. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 
OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1 hour  per  response,  including  the  time  for  reviewing  instructions,  searchmg  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway,  Suife  1204.  Arlington,  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503 

1.  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 

December  1999 

3.  REPORT  TYPE  AND  DATES  COVERED 

Conference  Publication 

4.  TITLE  AND  SUBTITLE 

1997  NASA  High-Speed  Research  Program  Aerodynamic  Performance 
Workshop 

Volume  I — Configuration  Aerodynamics 

5.  FUNDING  NUMBERS 

WU  537-07 

6.  AUTHOR(S) 

Daniel  G.  Baize,  Editor 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  ~ 

NASA  Langley  Research  Center 
Hampton,  VA  23681-2199 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

L-17916C 

9.  SPONSORINGWIONITOR1NG  AGENCY  NAME(S)  AND  ADDRESS(ES) 

National  Aeronautics  and  Space  Administration 
Washington,  DC  20546-0001 

10.  SPONSORINGVMONITORING 
AGENCY  REPORT  NUMBER 

NASA/CP- 1999-209691/ 
VOL1/PT2 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT  1 

1 12b.  DISTRIBUTION  CODE 

Unclassified-Unlimited 

Subject  Category  02  Distribution:  Nonstandard 

Availability:  NASA  CASI  (301)  621-0390 

13.  ABSTRACT  (Maximum  200  words) 

The  High-Speed  Research  Program  and  NASA  Langley  Research  Center  sponsored  the  NASA  High-Speed 
Research  Program  Aerodynamic  Performance  Workshop  on  February  25-28,  1997.  The  workshop  was  designed  to 
bring  together  NASA  and  industry  High-Speed  Civil  Tranport  (HSCT)  Aerodynamic  Performance  technology 
development  participants  in  area  of  Configuration  Aerodynamics  (transonic  and  supersonic  cruise  drag  prediction 
and  minimization),  High-Lift,  Flight  Controls,  Supersonic  Laminar  Flow  Control,  and  Sonic  Boom  Prediction.  The 
workshop  objectives  were  to  (1)  report  the  progress  and  status  of  HSCT  aeroydamic  performance  technology 
development;  (2)  disseminate  this  technology  within  the  appropriate  technical  communities;  and  (3)  promote 
synergy  among  the  scientist  and  engineers  working  HSCT  aerodynamics.  In  particular,  single-  and  multi-point 
optimized  HSCT  configurations,  HSCT  high-lift  system  performance  predictions,  and  HSCT  Motion  Simulator 
results  were  presented  along  with  executive  summaries  for  all  the  Aerodynamic  Performance  technology  areas. 


14.  SUBJECT  TERMS 

High-Speed  Research;  High-Speed  Civil  Transport 

15.  NUMBER  OF  PAGES 

941 

16.  PRICE  CODE 

A99 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 
Unclassified 

19.  SECURFTY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 

20.  LIMITATION 
OF  ABSTRACT 

UL 

NSN  7540-01*280-5500 

Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39-18 
298-102 

